機械機構(gòu)計算機輔助分析

1、 二分法原理 通過每次把f(x)的零點所在小區(qū)間收縮一半的方法，使區(qū)間的兩個端點逐步迫近函數(shù)的零點，以求得零點的近似值，這種方法叫做二分法。 二分法也稱對分法或區(qū)間套法，是逐次逼近一個實根的簡單方法。 02xx21.用二分法求解方程二分法程序框圖 matlab函數(shù)m文件如下：vfunction x=erfen(fx,xa,xb,n,eps)vx=xa;fa=eval(fx);vx=xb;fb=eval(fx);v disp( n xa xb xc fc );vfor i=1:n xc=(xa+xb)/2;x=xc;fc=eval(fx);v x=i,xa,xb,xc,fc;vdisp(x),i

2、f fc*fa0 xb=xc; else xa=xc; end if (xb-xa)/2eps,break,endendfx=x2-x-2;erfen(fx,-1.5,0.5,20,10(-5)ans = -1.0000fx=x2-x-2;erfen(fx,1.5,2.5,20,10(-5)ans = 2.0000牛頓法原理 牛頓法是解非線性方程（組）的一種重要迭代方法。牛頓法的一般公式：)( / )f(n1nnnxfxxx 牛頓法也常稱作切線法，其求單實根的收斂速度較快，但求重實根的收斂速度較快。2、用牛頓法求解方程0e-ex牛頓法程序框圖 matlab函數(shù)m文件如下：function x=

3、newton(fx,dfx,x0,eps,n)x=x0; f0=eval(fx); df0=eval(dfx);n=0; disp( n xn xn+1 fn+1 );while n=n x1=x0-f0/df0; x=x1;f1=eval(fx); x=n,x0,x1,f1; disp(x); if abs(x0-x1)eps fprintf(the procedure was successful.) return else n=n+1; x0=x1;f0=f1; endendif n=n+1 fprintf(the method failed after n iterations. )

4、endfx=exp(1)-exp(x);dfx=-exp(x);newton(fx,dfx,0.5,10(-5),100) the procedure was successful.ans = 1.0000 梯形法原理 給定步長，確定梯形的高之后，求出所有梯形的面積之和，再根據(jù)定積分的幾何意義可知曲頂高梯形的面積就是定積分的值。復(fù)化梯形公式為：niiixfxfhft112)()()(nabh/ )( 其中,n是區(qū)間的等分數(shù)。ba,3、用梯形法求解10sinxdx梯形法程序框圖 matlab函數(shù)m文件如下：vfunction y=t_quad(a,b,n)v h=(b-a)/n;v x1=a;v

5、 s=0;v for i=1:nv t=1/2*h*(sin(x1)+sin(x1+h);v s=s+t;v x1=x1+h;vendv svt_quad(0,1,100)vs = 0.4597歐拉法原理 歐拉（euler）法（亦稱切線法）是基于查分法的最簡單的顯式單步方法。盡管其精度較低，但算法簡單，起源最早，也是其它較復(fù)雜數(shù)值解法的一個很好的向?qū)А?歐拉法公式：),(1jjjjtyhfyy) 1,.,2 , 1 , 0(nj4、歐拉法和改進求歐拉法求：微分方程的解y(0)=1,h=0.2,區(qū)間0,1。yx2y y歐拉法程序框圖 matlab函數(shù)m文件如下：function x,y=eule

6、r(fun,x0,y0,xfine,h) n=(xfine-x0)/h; x(1)=x0;y(1)=y0; for n=1:n x(n+1)=x(n)+h; y(n+1)=y(n)+h*feval(fun,x(n),y(n);endfun=inline(y-2*x/y,x,y);x,y=euler(fun,0,1,1,0.2)x = 0 0.2000 0.4000 0.6000 0.8000 1.0000y =1.0000 1.2000 1.3733 1.5315 1.6811 1.8269改進歐拉法：改進歐拉法程序框圖 matlab函數(shù)m文件如下：function x,y=euler_g(f

7、un,x0,y0,xfine,h) n=(xfine-x0)/h; x(1)=x0;y(1)=y0; for n=1:n x(n+1)=x(n)+h; ym=y(n)+h*feval(fun,x(n),y(n); y(n+1)=y(n)+h/2*(feval(fun,x(n),y(n)+feval(fun,x(n+1),ym);endfun=inline(y-2*x/y,x,y);x,y=euler_g(fun,0,1,1,0.2)x = 0 0.2000 0.4000 0.6000 0.8000 1.0000y = 1.0000 1.1867 1.3483 1.4937 1.6279 1.7

8、542黃金分割法程序框圖輸入?yún)^(qū)間a,b，精度a1=b-0.618*(b-a); y1=f(a1);a2=a+0.618*(b-a); y2=f(a2);a2-a1y1停a=0.5(a1+a2)b=a2;a2=a1;y2=y1;a1=b-0.618*(b-a)y1=f(a1);a=a1;a1=a2;y1=y2;a2=a+0.618*(b-a)y2=f(a2);yyn 5、用黃金分割法求解方程f(x)=3x4-16x3+30 x2-24x+8在區(qū)間-10，10，=0.01 matlab函數(shù)m文件如下：function y=gold(f,a,b,eps) while abs(b-a)eps x1=a

9、+0.618*(b-a); x2=a+0.382*(b-a); f1=3*x14-16*x13+30*x12-24*x1+8 ; f2=3*x24-16*x23+30*x22-24*x2+8 ; if f1f2 b=x1;x1=x2; x2=a+0.382*(b-a); else a=x2;x2=x1; x1=a+0.618*(b-a); end endx=(x1+x2)/2運行結(jié)果如下：f=3*x4-16*x3+30*x2-24*x+8;gold(f,-10,10,0.01)x = 1.99846、斐波那契法求：f(x)=x+20/x在0.2，1 上的極值，=0.001 function y

10、=fibonacci(a,b,eps,n) f(1)=1;f(2)=1; for i=3:n f(i)=f(i-1)+f(i-2); x1=a+f(i-1)/f(i)*(b-a); x2=a+f(i-2)/f(i)*(b-a); while abs(x2-x1)eps f1=x1+20/x1; f2=x2+20/x2; if f1f2 b=x1;x1=x2; x2=a+f(i-2)/f(i)*(b-a); else a=x2;x2=x1; x1=a+f(i-1)/f(i)*(b-a); end end end x=(x1+x2)/2fibonacci(0.2,1,0.001,50)x = 0.

11、99897、用進退法求f(x)= 3x4-16x3+30 x2-24x+8 ，a0=0.5，h=0.5，=0.001 matlab函數(shù)m文件如下：function y=jintui(a0,h,eps) x1=a0; x2=x1+h; i=1;j=1; while(abs(x2-x1)eps) f1=3*x14-16*x13+30*x12-24*x1+8; f2=3*x24-16*x23+30*x22-24*x2+8; if(f1f2) x1=x2; x2=x1+2i*h; i=i+1; else x2=x1; x1=x2-h/4j; j=j+1; end endx=(x1+x2)/2f= 3*

12、x4-16*x3+30*x2-24*x+8運行結(jié)果如下：jintui(0.5,0.5,0.001)x = 1.9998f = 3.5751e-0078、質(zhì)點在振動分離篩上的分析(1)推導(dǎo)出k1,k2,k3。(2)分析=-25，-10，0，10，25，推導(dǎo)出=030的變化范圍。左移右移躍起的k的變化范圍。繪制=030時，k的變化規(guī)律。 （1）篩面上任一點位移： r為oa的長度。速度方程：加速度方程：物料在篩面上受到的牽連慣性力為： 當(dāng)曲柄在左半軸周（2,3象限)時，其慣性力的方向指向x的負方向，摩擦力的方向指向d。當(dāng)物料向下移動且沒有躍起時，則： 為物料相對篩面的加速度。 整理得：由于-1 co

13、swt 0,可得：0tansin)cos(cos0cos)sin(cosf22nnrfamgawtmrwmaamgawtmrwwtrwxcos2 wtrwxsinwtrw coss2wtrcosxra0)sin(g)cos(cos)sin(coscos22aawtrwawtmrwamgfn12)cos()sin(kaagrw當(dāng)曲柄在右半軸周（1,4象限)時，慣性力的方向指向x的正方向，摩擦力的方向指向c。當(dāng)物料向上移動且沒有躍起時，則：由于-1 coswt 0,可得：物料躍起條件：當(dāng) 時，由 式可以看出，只有coswt 時，只有coswt0時，才有可能 ，即物料在曲柄轉(zhuǎn)至右半周時躍起。（2 2

14、）clear; clear; clcclc; ;a=input(a=input(輸入角度：輸入角度：); );b b=a=a* *pi/180pi/180; ; c c=linspace(0.0001,30=linspace(0.0001,30* *pi/180,100); pi/180,100); d d=30=30* *(pi/180); (pi/180); k1=sin(k1=sin(d d- -b b)./cos()./cos(c c- -b b+ +d d); ); k2=sin(k2=sin(d d+ +b b)./cos()./cos(c c- -b b- -d d); );k3

15、=cos(k3=cos(b b)./abs(sin()./abs(sin(c c- -b b););0fna322)sin(cos0cos1cos)sin(cosrkaagrwwtagawtwc1=c*180/pi;plot(c1,k1,-b,c1,k2,-g,c1,k3,-r);xlabel(c);ylabel(k);title(k1,k2,k3曲線圖曲線圖);legend(k1 ,k2 ,k3 ); axis(0 30 0 8); 9（1）對于空間vx,vy,vz的表達式。（2）聯(lián)合收割谷物清選（小石頭，飽滿谷物，不成熟谷物，短秸桿）.其中v0=5,=10m/s,=10清涼室高度是0.3m

16、,u=0-30m/s，速度分別是25 m/s， 16 m/s ，14 m/s的位移曲線。（1）氣流的速度記為u，vl記為質(zhì)點的飄浮速度，則質(zhì)點相對于氣流的相對速度v的表達式為：轉(zhuǎn)換成三維微分方程：uvvfgmdtvdmgvuvvvvgdtdvvuvvvvgdtdvuvuvvvvgdtdvzxzylzyxzylyxxzylx222222222222)()()()(質(zhì)點初始速度v0在x、y、z方向上的分量：將時間分為若干小等分，每等分為t，求出t時間后的速度，再依據(jù)已知速度求出下一個t后的速度，以此類推，可以求所有點的速度。 令：sinsincoscoscos000000vvvvvvzyxgvuv

17、vvvgvvvfvuvvvvgvvvfuvuvvvvgvvvfzxzylzyxyxzylzyxxxzylzyx222232222222221)(),()(),()()(),(),(),(2),(),(2),(),(2t111111111111332211iiiiiiiiiiiiiiiiiiiiiiiizyxzyxzzzyxzyxyyzyxzyxxxvvvfvvvftvvvvvfvvvftvvvvvfvvvfvv在微小的t時間內(nèi),其速度可看作初始速度和終結(jié)速度的平均值,物料的位移量為：質(zhì)點的位移：（2）氣流分別為25m/s，16m/s，14m/s的位移曲線iiiiiiiiizzzyyyxxx11

18、12/)(2/)(2/)(111tvvztvvytvvxiiiiiizziyyixxifor j=1:1:4 u=0; if j=1 xl=24; elseif j=2 xl=10; elseif j= 3 xl=6.5; else xl=3; end clear; clc; w=30;h=0.3;l=1;v0=5;a=10*pi/180; b=a;g=9.8; t=0.001;h=w/2000; for i=1:1:2000 vx=v0*cos(a); vy=-v0*sin(a); u=u+h; x1=0; y1=0; while x1=0.3 & y1=1 f1=-g/xl2*sq

19、rt(vx+u*sin(b)2+(vy-u*cos(b)2)*(vx+u*sin(b)+g; f2=-g/xl2*sqrt(vx+u*sin(b)2+(vy-u*cos(b)2)*(vy-u*cos(b); vx1=vx+t*f1; vy1=vy+t*f2; f3=-g/xl2*sqrt(vx1+u*sin(b)2+(vy1-u*cos(b)2)*(vx1+u*sin(b)+g; f4=-g/xl2*sqrt(vx1+u*sin(b)2+(vy1-u*cos(b)2)*(vy1-u*cos(b); vx2=vx+t/2*(f1+f3); vy2=vy+t/2*(f2+f4); dx=(vx2+

20、vx)*t/2; dy=(vy2+vy)*t/2; vx=vx1; vy=vy1; x2=x1+dx; y2=y1+dy; x1=x2; y1=y2; end if j=1 v1(i)=y2; elseif j=2 v2(i)=y2; elseif j=3 v3(i)=y2; else v4(i)=y2; end endendwhile 1 fengsu=input(請輸入風(fēng)速請輸入風(fēng)速:); if fengsu0 break; else num=input(請輸入輸出參數(shù)個數(shù)請輸入輸出參數(shù)個數(shù):); n=round(fengsu/h); if n=0n=1; end out(1)=v1(n)

21、; out(2)=v2(n); out(3)=v3(n); out(4)=v4(n); if num=1 location=out(1) end if num=2 location=out(1:2) end if num=3 location=out(1:3) end if num=4 location=out(1:4) end end end運行結(jié)果請輸入風(fēng)速:25請輸入輸出參數(shù)個數(shù):4location = -0.0311 0.0786 0.3931 1.0029請輸入風(fēng)速:16請輸入輸出參數(shù)個數(shù):4location = -0.0423 0.0008 0.0884 1.0061請輸入輸出參數(shù)個數(shù):4location = -0.0442 -0.0111 0.0520 1.0023 10、曲柄滑塊機構(gòu)運動學(xué)分析clear;clc;w=input(角速度(轉(zhuǎn)數(shù)):r=);r=input(oa桿長度:r=);l=input(ab桿長度:l=);h=input