วันอังคารที่ 30 มกราคม พ.ศ. 2561
วันพฤหัสบดีที่ 25 มกราคม พ.ศ. 2561
วันพุธที่ 24 มกราคม พ.ศ. 2561
วันอังคารที่ 23 มกราคม พ.ศ. 2561
วันพุธที่ 17 มกราคม พ.ศ. 2561
วันอังคารที่ 16 มกราคม พ.ศ. 2561
วันจันทร์ที่ 15 มกราคม พ.ศ. 2561
สิ่งที่ อ.สอนในคาบ
การค้นหาไฟล์ใน Google
พิมพ์
Blasting ext pdf (หาไฟล์ pdf แบบระบุชื่อ)
* inurl:www...... ext:pdf (หาไฟล์ pdf แบบไม่ระบุชื่อจากเว็บ)
กูเกิ้ลวาดกราฟ
ไม่ต้องพิมพ์ Y=.....
3x^2 = 3*x^2
Scilab
Clear(a) ......เคลียเฉพาะ a
Clear(all) = clc .........เคลียทุกตัวแปร
u+v
u-v
u*v
u*v'
u'*v
u(3)
u(2)+v(1)
a=1;b=2;c=3;r=[a,b,c]
s={-1,2}
letters={'a','b','c','d'}
A=[1,2,3;4,5,6;1,-1,0]
B=[1,1,1,
2,3,-1,
5,7,-2]
u=[1,3,-5];v=[4,2,3];w=[-1,0,1];
C=[u;v;w]
D=[u'v'w']
r=[u,v,w]
A+B
C-D
A*B
B*A
C*u
D*v'
rank(A)
inv(A)
cond(B)
det(C)
A*inv(A)
inv(B)*B
spec(A)
trace(C)
-----------------------
เว้นวรรคแทน , ได้ แต่โอกาสพลาดง่ายไม่แนะนำ
x={.5 .7 .9 1.3 1.7 1.8]
y=[.1 .2 .75 1.5 2.1 2.4]
วาดกราฟ 2D เป็นจุด
(-1 เป็น +)
(-2 เป็น X)
(-3 เป็น วงกลมมีบวก)
-4 เป็นหน้าจั่วทึบ
-5 เป็นหน้าจั่วโปร่ง
-6 เป็นสามเหลี่ยมโปร่งชี้ขึ้น
-7 เป็นสามเหลี่ยมโปร่งชี้ลง
-8 เป็นหน้าจั่วโปร่งมีบวก
-9 เป็นวงกลมโปร่ง
-10 เป็นดอกจันมีบวก
-11 เป็นสี่เหลี่ยมโปร่ง
plot2d(x,y,style=-1)
ไม่ต้องปิดกราฟแรก คำสั่งต่อมาจะได้เส้นเพิ่ม ถ้าปิดไปแล้วจะได้แต่เส้น
plot2d(x,y)
พลอตกราฟรูปหัวใจ
x1=linspace(-1,0,50);
y1=-x1+sqrt(3-3*x1.^2);
y2=-x1-sqrt(3-3*x1.^2);
title('heart shape curve');
plot(x1,y1,'r');
plot(-x1,y1,'r');
plot(-x1,y2,'r');
plot(x1,y2,'r');
------------------------------
n=[3 5 9 17 33 65]; s=[.257 .0646 .0151 3.96e-3 9.78e-4 2.45e-4]; plot2d(n,s,style=-1) plot2d(n,s,style=5)
พิมพ์ใน Scinote
n=[3 5 9 17 33 65];
s=[.257 .0646 .0151 3.96e-3 9.78e-4 2.45e-4];
plot2d(n,s,style=4,logflag='ll')
plot2d(n,s,style=4,logflag='nl')
plot2d(n,s,style=4,logflag='ln')
plot2d(n,s,style=-1) ใส่จุด
ln คือเซมิล็อก y เป็น log แกน x เป็น normal
nl คือเซมิล็อก x เป็น log แกน y เป็น normal
xgrid(10); ใส่เส้นกริดในกราฟ
clear
n=10
p=100
i=5
for j=1:n
f(j)=p*(1+(i/100))^j;
printf('\nFuture value for each period.%3f\n',f(j))
j=j+1;
end
clear
n=input('input number of period...');
p=input('input Present value...');
i=input('input interest rate in percent...');
for j=1:n
f(j)=p*(1+(i/100))^j;
printf('\nFuture value for each period.%3f\n',f(j))
j=j+1;
end
clear
n=input("จำนวนปี...");
p=input("เงินฝาก...");
i=input("อัตราดอกเบี้ยต่อปี...");
for j=1:n
f(j)=p*(1+(i/100))^j;
printf('\nจำนวนเงินทั้งหมดที่ได้รับในแต่ละปี..%3f\n',f(j))
j=j+1;
end
s=[.257 .0646 .0151 3.96e-3 9.78e-4 2.45e-4];
plot2d(n,s,style=-1)
plot2d(n,s,style=5)
------------------------------------------------------------------------------------------
clear
a=input("แร่ป้อน=","string");
d=[1:10:1000]; ทำให้มันเป็นค่า 1 ถึง1000 โดย ห่างช่วงละ 10
u=input("ความหนืดของตัวกลาง(kg/m.s.)=");
g=input("แรงโน้มถ่วง(m/s^2)=");
s1=input("ความหนาแน่นของอนุภาคตัวกลาง1(kg/m^3)=");
s2=input("ความหนาแน่นของอนุภาคตัวกลาง2(kg/m^3)=");
S=(s1-s2);
Vt=(1/18)*(S/u)*g*((d*10^-6)^2);
xtitle("การตกตัวของอนุภาค"); ใส่ชื่อกราฟ
plot2d(d,Vt,style=9);
xgrid(10);
------------------------------------------------------------------------
D=input ("Blast Hole diameter in mm.=");
BH=D/15
printf('\nBench height >= %.2f metres\n',BH);
B=30*D
printf('\nBurden = %.2f\n',B);
S=1.5*B
printf('\nSpacing = %.2f\n',S);
Sd3=3*D
Sd15=15*D
printf('\nSubdrill Range= %.2f to %.2f\n',Sd3,Sd15);
SD=input ("Please select the Subdrill you want to use = ");
C=20*D
printf('\nRecommend Charge length >= %.2f\n',C);
SL=20*D*B
printf('\nStemming >= %.2f\n',SL);
BSR=BH/B
printf('\nBurden stiffness ratio = %.5f\n',BSR);
if(BSR>2&BSR<3.5)
BSR='good'
else BSR='very good'
end
SMS10=D/10
SMS20=D/20
printf('\nStemming material size= %.2f to %.2f\n',SMS10,SMS20);
L=BH+SD
printf('\nHole length = %.2f\n',L);
C2=L-SL
printf('\nCharge length = %.2f\n',C2);
N=input ("Number of Holes in a blast=");
DS=input ("Density of rock in t/m^3 =");
V=B*S*BH*N
printf('\nBlast Volume = %.2f\n',V);
T=V*DS
printf('\nBlasted tonnes = %.2f\n',T);
Vb=3.14*((D^2)/4000)*L
printf('\nVolume of blasthole = %.2f\n',Vb);
Ex=input ("Explosive density=");
Kg=Vb*Ex
printf('\nMass of explosive per hole = %.2f Kg\n',Kg);
PF=(Kg*N)/T
printf('\nPower Factors = %.2f Kg/m^3\n',PF);
-----------------------------------------------------------------------
ทำทิ้งไว้
clear;
clc
printf('Do you want shape 2D or 3D?');
T=input("If you want 2D Print 20 , 3D Print 30 =");
if T==20 then
printf('You want to find area of "Circle" or "Annulus"?');
T=input("If you want Circle Print 21 , Annulus Print 22 =");
if(T==21);
3D
if(T==30);
printf('You want to find area of "General cone or Pyramid" or "Frustum of cone"?');
T=input("If you want General cone or Pyramid Print 31 , Frustum of cone Print 32 =");
V=(3.14/3)*((r^2)+(R*r)+(R^2))*h
S=3.14*s*(R+r)+(3.14*(r^2))+(3.14*(R^2))
พิมพ์
Blasting ext pdf (หาไฟล์ pdf แบบระบุชื่อ)
* inurl:www...... ext:pdf (หาไฟล์ pdf แบบไม่ระบุชื่อจากเว็บ)
กูเกิ้ลวาดกราฟ
ไม่ต้องพิมพ์ Y=.....
3x^2 = 3*x^2
Scilab
Clear(a) ......เคลียเฉพาะ a
Clear(all) = clc .........เคลียทุกตัวแปร
exp a ...........exponantial a
save a........คือเก็บ a ไว้ แม้เคลียแล้ว แต่ยังสามาใช้คำสั่ง load a เรียกค่ากลับมาได้
ทำแนวตั้ง แนวนอน
x=-10.0:0.5:10.0
xf=x'
y=sin(x*%pi/10)
plot(x,y)
help plot
u=[-1,2,3]
v=[6,-10,0]u+v
u-v
u*v
u*v'
u'*v
u(3)
u(2)+v(1)
a=1;b=2;c=3;r=[a,b,c]
s={-1,2}
letters={'a','b','c','d'}
A=[1,2,3;4,5,6;1,-1,0]
B=[1,1,1,
2,3,-1,
5,7,-2]
u=[1,3,-5];v=[4,2,3];w=[-1,0,1];
C=[u;v;w]
D=[u'v'w']
r=[u,v,w]
A+B
C-D
A*B
B*A
C*u
D*v'
rank(A)
inv(A)
cond(B)
det(C)
A*inv(A)
inv(B)*B
spec(A)
trace(C)
-----------------------
เว้นวรรคแทน , ได้ แต่โอกาสพลาดง่ายไม่แนะนำ
x={.5 .7 .9 1.3 1.7 1.8]
y=[.1 .2 .75 1.5 2.1 2.4]
วาดกราฟ 2D เป็นจุด
(-1 เป็น +)
(-2 เป็น X)
(-3 เป็น วงกลมมีบวก)
-4 เป็นหน้าจั่วทึบ
-5 เป็นหน้าจั่วโปร่ง
-6 เป็นสามเหลี่ยมโปร่งชี้ขึ้น
-7 เป็นสามเหลี่ยมโปร่งชี้ลง
-8 เป็นหน้าจั่วโปร่งมีบวก
-9 เป็นวงกลมโปร่ง
-10 เป็นดอกจันมีบวก
-11 เป็นสี่เหลี่ยมโปร่ง
plot2d(x,y,style=-1)
ไม่ต้องปิดกราฟแรก คำสั่งต่อมาจะได้เส้นเพิ่ม ถ้าปิดไปแล้วจะได้แต่เส้น
plot2d(x,y)
พลอตกราฟรูปหัวใจ
x1=linspace(-1,0,50);
y1=-x1+sqrt(3-3*x1.^2);
y2=-x1-sqrt(3-3*x1.^2);
title('heart shape curve');
plot(x1,y1,'r');
plot(-x1,y1,'r');
plot(-x1,y2,'r');
plot(x1,y2,'r');
------------------------------
n=[3 5 9 17 33 65]; s=[.257 .0646 .0151 3.96e-3 9.78e-4 2.45e-4]; plot2d(n,s,style=-1) plot2d(n,s,style=5)
พิมพ์ใน Scinote
n=[3 5 9 17 33 65];
s=[.257 .0646 .0151 3.96e-3 9.78e-4 2.45e-4];
plot2d(n,s,style=4,logflag='ll')
plot2d(n,s,style=4,logflag='nl')
plot2d(n,s,style=4,logflag='ln')
plot2d(n,s,style=-1) ใส่จุด
ln คือเซมิล็อก y เป็น log แกน x เป็น normal
nl คือเซมิล็อก x เป็น log แกน y เป็น normal
xgrid(10); ใส่เส้นกริดในกราฟ
clear
n=10
p=100
i=5
for j=1:n
f(j)=p*(1+(i/100))^j;
printf('\nFuture value for each period.%3f\n',f(j))
j=j+1;
end
clear
n=input('input number of period...');
p=input('input Present value...');
i=input('input interest rate in percent...');
for j=1:n
f(j)=p*(1+(i/100))^j;
printf('\nFuture value for each period.%3f\n',f(j))
j=j+1;
end
clear
n=input("จำนวนปี...");
p=input("เงินฝาก...");
i=input("อัตราดอกเบี้ยต่อปี...");
for j=1:n
f(j)=p*(1+(i/100))^j;
printf('\nจำนวนเงินทั้งหมดที่ได้รับในแต่ละปี..%3f\n',f(j))
j=j+1;
end
กดรันใน Scinote แล้วพิมพ์ n=10,P=100,i=5
clear
a=input("เงินเดือน...");
b=input("จำนวนปี...");
i=input("อัตราดอกเบี้ยต่อปี...");
for j=1:b
p=a*12*b
f(j)=p*(1+(i/100))^j;
printf('\nจำนวนเงินทั้งหมดที่ได้รับในแต่ละปี.%d..%3f\n',j,f(j))
j=j+1;
end
การคำนวณแฟคทอเรียล
function fact=factorial(n)
fact=1
for k=1:n
fact =k*fact
end
endfunction
Run ใน Scinote แล้วพิมพ์
factorial(5) ใน Scilab ตามใจใส่เลขไร
เปลี่ยนชื่อฟังก์ชั่นได้นะ
function fact=Jinna(n)
fact=1
for k=1:n
fact =k*fact
end
endfunction
เปลี่ยนชื่อฟังก์ชั่นได้นะ
function fact=Jinna(n)
fact=1
for k=1:n
fact =k*fact
end
endfunction
Run ใน Scinote แล้วพิมพ์
Jinna(5)
-----------------------------
การใช้ฟังก์ชั่น IF ELSE
function s=signum(x)
if(x>0)
s='positive'
elseif(x<0)
s='negative'
else
s=0
end
endfunction
-----------------------------
การใช้ฟังก์ชั่น IF ELSE
function s=signum(x)
if(x>0)
s='positive'
elseif(x<0)
s='negative'
else
s=0
end
endfunction
Run ใน Scinote แล้วพิมพ์
signum(-123)
signum(123)
signum(0)
-------------------------------------------
clear
d=input("เส้นผ่านศูนย์กลางอนุภาค(m)=");
u=input("ความหนืดของตัวกลาง(kg/m.s.)=");
g=input("แรงโน้มถ่วง(m/s^2)=");
s1=input("ความหนาแน่นของอนุภาคตัวกลาง1(kg/m^3)=");
s2=input("ความหนาแน่นของอนุภาคตัวกลาง2(kg/m^3)=");
S=(s1-s2);
Vt=(1/18)*(S/u)*g*(d^2);
printf('\nความเร็วในการตกตัว=%.11f\n',Vt)
อันนี้คือเขียนสมการ
-------------------------------------------
clear
d=input("เส้นผ่านศูนย์กลางอนุภาค(m)=");
u=input("ความหนืดของตัวกลาง(kg/m.s.)=");
g=input("แรงโน้มถ่วง(m/s^2)=");
s1=input("ความหนาแน่นของอนุภาคตัวกลาง1(kg/m^3)=");
s2=input("ความหนาแน่นของอนุภาคตัวกลาง2(kg/m^3)=");
S=(s1-s2);
Vt=(1/18)*(S/u)*g*(d^2);
printf('\nความเร็วในการตกตัว=%.11f\n',Vt)
อันนี้คือเขียนสมการ
-----------------------------------------------------------------------------
x_vais=[.5 2 2.5 3 3.5 4 4.5 5]
y_vais=[0 2.7 4 5 6 7 8 9]
plot2d(x_vais,y_vais,leg='Line of Best Fit',style=10)
plot2d(x_vais,y_vais,style=-5)
xtitle('Data.Analysis')
------------------------------------
x_vais=[.5 2 2.5 3 3.5 4 4.5 5]
y_vais=[0 2.7 4 5 6 7 8 9]
plot2d(x_vais,y_vais,leg='Line of Best Fit',style=10')
plot2d(x_vais,y_vais,style=-5)
xtitle('Data Analysis')
a=[.5 2 3 4 5 6]
b=[0 2.5 3.5 4.5 5.5 6.5]
plot2d(a,b,style=5)
plot2d(a,b,style=-9)
xtitle('Data Analysis')
-------------------------------------
n=[3 5 9 17 33 65];s=[.257 .0646 .0151 3.96e-3 9.78e-4 2.45e-4];
plot2d(n,s,style=-1)
----------------------------------------------------------------------------------------------
clear
a=input("แร่ป้อน=","string");
d=input("เส้นผ่านศูนย์กลางอนุภาค(m)=");
u=input("ความหนืดของตัวกลาง(kg/m.s.)=");
g=input("แรงโน้มถ่วง(m/s^2)=");
s1=input("ความหนาแน่นของอนุภาคตัวกลาง1(kg/m^3)=");
s2=input("ความหนาแน่นของอนุภาคตัวกลาง2(kg/m^3)=");
S=(s1-s2);
Vt=(1/18)*(S/u)*g*(d^2);
printf('\nความเร็วในการตกตัวของ %s เท่ากับ%.11f\n',a,Vt)
การให้มันประกาศคำ
clear
a=input("แร่ป้อน=","string");
d=input("เส้นผ่านศูนย์กลางอนุภาค(m)=");
u=input("ความหนืดของตัวกลาง(kg/m.s.)=");
g=input("แรงโน้มถ่วง(m/s^2)=");
s1=input("ความหนาแน่นของอนุภาคตัวกลาง1(kg/m^3)=");
s2=input("ความหนาแน่นของอนุภาคตัวกลาง2(kg/m^3)=");
S=(s1-s2);
Vt=(1/18)*(S/u)*g*(d^2);
printf('\nความเร็วในการตกตัวของ %s เท่ากับ%.11f\n',a,Vt)
การให้มันประกาศคำ
clear
a=input("แร่ป้อน=","string");
d=[1:10:1000]; ทำให้มันเป็นค่า 1 ถึง1000 โดย ห่างช่วงละ 10
u=input("ความหนืดของตัวกลาง(kg/m.s.)=");
g=input("แรงโน้มถ่วง(m/s^2)=");
s1=input("ความหนาแน่นของอนุภาคตัวกลาง1(kg/m^3)=");
s2=input("ความหนาแน่นของอนุภาคตัวกลาง2(kg/m^3)=");
S=(s1-s2);
Vt=(1/18)*(S/u)*g*((d*10^-6)^2);
xtitle("การตกตัวของอนุภาค"); ใส่ชื่อกราฟ
plot2d(d,Vt,style=9);
xgrid(10);
------------------------------------------------------------------------
D=input ("Blast Hole diameter in mm.=");
BH=D/15
printf('\nBench height >= %.2f metres\n',BH);
B=30*D
printf('\nBurden = %.2f\n',B);
S=1.5*B
printf('\nSpacing = %.2f\n',S);
Sd3=3*D
Sd15=15*D
printf('\nSubdrill Range= %.2f to %.2f\n',Sd3,Sd15);
SD=input ("Please select the Subdrill you want to use = ");
C=20*D
printf('\nRecommend Charge length >= %.2f\n',C);
SL=20*D*B
printf('\nStemming >= %.2f\n',SL);
BSR=BH/B
printf('\nBurden stiffness ratio = %.5f\n',BSR);
if(BSR>2&BSR<3.5)
BSR='good'
else BSR='very good'
end
SMS10=D/10
SMS20=D/20
printf('\nStemming material size= %.2f to %.2f\n',SMS10,SMS20);
L=BH+SD
printf('\nHole length = %.2f\n',L);
C2=L-SL
printf('\nCharge length = %.2f\n',C2);
N=input ("Number of Holes in a blast=");
DS=input ("Density of rock in t/m^3 =");
V=B*S*BH*N
printf('\nBlast Volume = %.2f\n',V);
T=V*DS
printf('\nBlasted tonnes = %.2f\n',T);
Vb=3.14*((D^2)/4000)*L
printf('\nVolume of blasthole = %.2f\n',Vb);
Ex=input ("Explosive density=");
Kg=Vb*Ex
printf('\nMass of explosive per hole = %.2f Kg\n',Kg);
PF=(Kg*N)/T
printf('\nPower Factors = %.2f Kg/m^3\n',PF);
-----------------------------------------------------------------------
ทำทิ้งไว้
clear;
clc
printf('Do you want shape 2D or 3D?');
T=input("If you want 2D Print 20 , 3D Print 30 =");
if T==20 then
printf('You want to find area of "Circle" or "Annulus"?');
T=input("If you want Circle Print 21 , Annulus Print 22 =");
if(T==21);
r=input("Radius =");
d=2*r
printf('\nDiameter =%.2f \n',d);
d=2*r
printf('\nDiameter =%.2f \n',d);
A=3.14*(r^2);
C=2*3.14*r;
C=2*3.14*r;
printf('\nArea =%.2f \n',A);
printf('\nCircumference =%.2f \n',C);
elseif(T==22);
r=input("Inner radius =");
R=input("Outer radius =");
P=0.5*(r+R);
W=R-r;
A=3.14*((R^2)-(r^2));
printf('\nAverage Radius =%.2f \n',P);
printf('\nWidth =%.2f \n',W);
printf('\nArea =%.2f \n',A);
3D
if(T==30);
printf('You want to find area of "General cone or Pyramid" or "Frustum of cone"?');
T=input("If you want General cone or Pyramid Print 31 , Frustum of cone Print 32 =");
if(T==31);
if(T==32);
A=input("Area of base =");
h=input("Height =");
V=(1/3)*A*h;
printf('\nVolume =%.2f \n',V);
r=input("Top Radius =");
R=input("Base radius =");
h=input("Height =");
s=input("Slant height =");
S=3.14*s*(R+r)+(3.14*(r^2))+(3.14*(R^2))
printf('\nVolume =%.2f \n',V);
printf('\nSurface Area =%.2f \n',S);
!= คือเครื่องหมายไม่เท่ากับ
!= คือเครื่องหมายไม่เท่ากับ
----------------------------------------------------------
clc
x=1;
while x==1
D=input("How many Dimentions =?");
if D==2 then
ch1=input("Sector of circle press 1\nPararellogran press 2 ");
if ch1==1 then
r=input("radius= ?");
a=input("Angle in ra =?");
Ar=0.5*a*r*r;
s=a*r;
printf("Area = %.2f",Ar);
printf("\n Arc Length = %.2f",s);
elseif ch1==2 then
b=input("base = ?");
h=input("height = ?");
a1=input("Side = ")
Ar1=b*h;
P= (2*a1)+(2*b);
printf("Area = %.2f",Ar1);
printf("\n Perimeter = %.2f ",p);
else
disp("Error404");
end
elseif D==3 then
ch1=input("Right circular cone press 1\nQuare Pyramid press 2 ");
if ch1==1 then
r2=input("radius= ?");
h2=input("height = ?");
Vo1=0.33*3.14*r2*r2*h2;
s=(3.14*((r2^2)+(h2^2))^0.5)+(3.14*r2*r2);
printf("Volume = %.2f m3",Vo1);
printf("\n Surface Area = %.2f ",s);
elseif ch1==2 then
s=input("side = ?");
h=input("height = ?");
vo=0.33*(s*s)*h;
sf=s*(s+((s^2)+(4*h*h))^0.5);
printf("Vol = %.2f \n Surface Area = %.2f",vo,sf);
end
end
x=input("\n \n \n **** Repeat press 1 ****");
end
สมัครสมาชิก:
ความคิดเห็น (Atom)





