matlab程序设计上机练习(2010研)c

《matlab程序设计上机练习(2010研)c》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、MATLAB程序设计上机练习C1、按顺序进行如下的操作（1）用Matlab命令产生一个4阶魔术方阵，并计算与；（2）求的特征值；（3）计算的各列的总和与平均值；（4）若，求方程组的解。2、绘制下列曲线的图形：。在命令窗口输入：ezplot(cos(x^2),[0,2*pi]);93、假如有一组实测数据：求其回归直线，画回归直线图形并计算最小误差平方和。1）回归直线求解：>>A=[0.11;0.111;0.121;0.131;0.141;0.151;0.161;0.171;0.181;0.21]>>b=[42;

2、43.5;45;45.5;45;47.5;49;53;50;55]>>a=Aba=124.675329.3474>>symsx>>y=124.6753*x+29.3474;2）回归直线9>>ezplot('y',[0.1,0.2])3）残差平方和（最小误差平方和）;求残差的M文件：functionS=fangcha(x,y)b=zeros(1,10);fori=1:109b(i)=124.6753*x(i)+29.3474;endbc=b-yS=0;fori=1:10S=S+c(i)^2;end命令窗口;>

3、>a=[0.10.110.120.130.140.150.160.170.180.2];>>b=[4243.54545.54547.549535055];>>S=fangcha(a,b)计算结果;b=Columns1through9941.814943.061744.308445.555246.801948.048749.295450.542251.7890Column1054.2825c=Columns1through9-0.1851-0.4383-0.69160.05521.80190.54870.295

4、4-2.45781.7890Column10-0.7175S=14.09904、求下列各结果：（1）用Matlab展开：；（2）用Matlab求极限：。（1）命令窗口：>>symsx>>y=(x^2+1)*(x^4-x^2+1);>>f=expand(y)9结果：f=x^6+1（2）命令窗口：>>symsxn>>f=(1-x/n)^nf=(1-x/n)^n>>limit(f,n,inf)结果：ans=exp(-x)5、用Matlab命令求方程的根。命令窗口;>>p=[1-301];>>roots(p)结果：a

5、ns=2.879490.6527-0.53216、求微分方程的通解。命令窗口：>>symsxy>>dsolve('(x^2-1)*Dy+2*x*y-cos(x)=0')结果：ans=1/2*cos(x)/x+exp(-2*x/(x^2-1)*t)*C17、求微分方程组在初始条件下的特解，并画出图形。命令窗口：>>symsxyt>>S=dsolve('Dx+5*x+y=exp(2*t)','Dy-x-3*y=0','x(0)=-1','y(0)=1','t');>>x=S.x,y=S.y作图：9>>ezplot

6、(S.x)>>ezplot(S.y)99