MATLAB程序設(shè)計(jì)

1、一 多 項(xiàng) 式 擬 合多項(xiàng)式擬合是多項(xiàng)式運(yùn)算的一個(gè)重要組成部分，在工程及科研工作中都等到了廣泛的應(yīng)用。其實(shí)現(xiàn)一方面可以由矩陣的除法求解超定方程來(lái)進(jìn)行；另一方面在matlab中還提供了專門的擬合函數(shù)polyfit,其常用調(diào)用格式如下：polyfit(x,y,n):其中x, y為擬合數(shù)據(jù)，n為擬合多項(xiàng)式的階數(shù)。p,s=polyfit (x, y, n):其中p為擬合多項(xiàng)式系數(shù)常量，s為擬 合多項(xiàng)式系數(shù)常量的結(jié)構(gòu)信息。1 .用5階多項(xiàng)式對(duì)0, pi/2上的正弦函數(shù)進(jìn)行最小二乘擬合。x=0:pi/20:pi/2;y=sin(x);a=polyfit(x,y,5);x1=0:pi/30:pi*2;y1=

2、sin(x1);y2=a(1)*x1a5+a(2)*x1a4+a(3)*x1a3+a(4)*x1a2+a(5)*x1+a(6);plot(x1,y1,b-,x1,y2,r*)legend(原曲線,擬合曲線)axis(0,7,-1.2,4)2 .設(shè)原始數(shù)據(jù)為x,在11個(gè)點(diǎn)上測(cè)得的y值如下:x=-2.0 -1.6 -1.2 -0.8 -0.4 0 0.4 0.8 1.2 1.6 2.0;y=2.8 2.96 2.54 3.44 3.56 5.4 6.0 8.4 9.5 13.3 15;采用2階和10階多項(xiàng)式擬合，弁畫圖比較x=-2.0 -1.6 -1.2 -0.8 -0.4 0 0.4 0.8 1

3、.2 1.6 2.0;y=2.8 2.96 2.54 3.44 3.56 5.4 6.0 8.4 9.5 13.3 15;a,polyfit(x,y,2);x1=-2:0.01:2;y2=a(1)*x1.a2+a(2)*x1+a(3); plot(x,y, b- ,x1,y2, r-)1614 -12 -/ 10、.8 .6 lr4 -2 d-2-1.5-1-0.500.511.52x=-2.0 -1.6 -1.2 -0.8 -0.4 0 0.4 0.8 1.2 1.6 2.0;y=2.8 2.96 2.54 3.44 3.56 5.4 6.0 8.4 9.5 13.3 15;a=polyfi

4、t(x,y,10);x1=-2:0.01:2;y2=a(1)*x1.a10+a(2)*x1.a9+a(3)*x1.a8+a(4)*x1.a7+a(5)*x1.a6+a(6)*x1.a5 +a*x1.a4+a(8)*x1.a3+a(9)*x1.a2+a(10)*x1+a(11);plot(x,y, b- ,x1,y2, r-)300252015105-2-1.5-1-0.500.511.52二圖形的保持與子圖1.在4個(gè)子圖中繪制不同的三角函數(shù)圖。（二維子圖的繪制）x=0:0.1*pi:2*pi;subplot(2,2,1);plot(x,sin(x),-*);title( sin(x) );su

5、bplot(2,2,2);plot(x,cos(x), -o);title( cos(x) );subplot(2,2,3);plot(x,sin(x).*cos(x),-x);title( sin(x)*cos(x);subplot(2,2,4);plot(x,sin(x)+cos(x),-h);title( sin(x)+cos(x) )0.510.5 0-0.5-102468sin(x)cos(x)0-0.5 02468sin(x)*cos(x)sin(x)+cos(x)10-1-2 02對(duì)于函數(shù)z (x 1)2 (y 1)2, x的取值范圍02, y的取值范圍20,畫出四個(gè)三維子圖，其

6、中 1,1；2,1;1,2;0.8,0.5；要求每一個(gè)子圖有標(biāo)題。x,y=meshgrid(0:0.1:2,-2:0.1:0);subplot(2,2,1);z=(x-1)a2+(y+1)a2;mesh(z);title(子圖 1);subplot(2,2,2);z=2*(x-1)a2+(y+1)a2； mesh(z);title(子圖 2);subplot(2,2,3);z=(x-1)a2+2*(y+1)a2;mesh(z);title(子圖 3);subplot(2,2,4);z=0.8*(x-1)a2+0.5*(y+1)a2;mesh(z);title(子圖 4);子圖1子圖20 0子圖

7、3子圖4三.龍格-庫(kù)塔方法1 .用經(jīng)典的 r-k 方法求解 v =-2y+2x2+2x。其中(0=x=0.5 ) ,y(0)=1.編制函數(shù)文件(fun.m)fun.mfunction f=fun(x,y)f=-2*y+2*xa2+2*x在命令窗口輸入x,y=ode23(fun,0,0.5,1); x0.04000.09000.14000.39000.4400 0.49000.50000.92470.84340.77540.61050.6084 0.61540.61790 0.2900 0.3400 y1.00000.6440 0.62220.19000.71990.24000.6764再輸入：

8、plot(x,y,-o)可等到圖形如下:10.950.90.850.80.750.70.6500.050.10.150.20.250.30.350.40.450.52 .解剛性方程解如下剛性方程：(見課本251頁(yè))編制f.m文件function f=f(x,y)f=-2 1;998 -999*y+2*sin(x);999*(cos(x)-sin(x)在程序窗口輸入：tic;ode45(f,0 10,2,3);t1=toc t1 =3.4540tic;ode23(f,0 10,2,3);t1=toc t1 =2.9839tic;ode23s(f,0 10,2,3);t1=toc t1 =0.69

9、5732.521.510.50-0.5-1012345678910-1.523.用龍格庫(kù)塔方法求解常微分方程，初始條件為jx 4,s|t0 dt2 dt x t o 1編制函數(shù)文件(f.m)function ydot=f(t,y)ydot=y(2); 4 ;在命令窗口輸入：t,y=ode45(f,0:0.1:10,1,2);plot(t,y(:,1),t,y(:,2);25020015010050123456789104.設(shè)有一微分方程組x1x28s(t),已知t 0當(dāng)時(shí), x2(t) sin(2t)x2(0)0.5,求微分方程在t 0,50上的解，弁畫出xi編制函數(shù)文件(liao.m)fun

10、ction f=liao(t,x)f=x(2)+cos(t);sin(2*t);在命令窗口輸入：t,x=ode45(liao,0 50,0.5,-0.5);plot(x(:,1),x(:,2);xk0) 0.5,x2的軌跡。四,while/for 循環(huán)1.用for循環(huán)求1!+2!+3!+？。!的值sum=0;for i=1:20;prd=1;for k=1:i;prd=prd*k; endsum=sum+prd;endsumsum =2.5613e+018五.simulink電路仿真simulink是matlab軟件的擴(kuò)展，是實(shí)現(xiàn)動(dòng)態(tài)系統(tǒng)建模和仿真的一個(gè)軟件包。運(yùn)行一個(gè)仿真完整的步驟:設(shè)置仿真參數(shù)(2).啟動(dòng)仿真(3)仿真結(jié)果分析1、電阻電路的計(jì)算，電路如圖+ % -ra r4 rs i i-t+建立 simulink 模型，r1=2, r2=4, r3=12, r4=4, r5=12, r6=4,r7=2, us=10