资源描述：
《数学建模matlab练习4》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、MATLAB练习四1.绘制两个一元函数的曲线图(例如y二sin(x)、y二cos(x)),加上标注说明，包括坐标、抬头、曲线说明等。%lianxi41.mx=0:pi/30:2*pi；y=sin(x)；plot(xzy,*-ob1z'markersize1,10)；holdonyl=cos(x)；plot(x,yl,'-.dm','markersize1,10)；legend(1sinx1,'cosx1)；title(1sin和cos函数的曲线图')xlabel(1x*)；ylabel(1y1)；2.寻找一组数据拟合多项式，写出M文件，包括拟
2、合的次数，数值的计算，数学表达式的呈现，图形的呈现等。%lianxi42x=0:pi/20:pi;y=sin(x)；n=input(*n=')；p=polyfit(x,y,n);xl=0:pi/20:pi*2;yl=sin(xl)；y2=polyval(p,xl)plot(xlzyl,'-ob1,xl,y2z':dm*)；legend(1原曲线J】拟合曲线「);xlabel(1x1)；ylabel(*y1)；title(1sin的拟合曲线和原曲线');axis([-0.2,6.5,-1.2,1.2])1.求解一个代数方程，写出M文件，应该包括
3、函数式M文件和命令式M文件。%fc3.mfunctionf=fc3(x)f=13*x-l=01；%lianxi43clc；clear；close；sprintf('原方程和Jt解析解的结果为：，)s=solve(fc3)sprintf(•保留10位有效数字的结果为：Jdigits(10)sprintf('数值解的结果为：Tx=vpa(s)2.求解一个代数方程组，写出M文件，应该包括函数式M文件和命令式M文件。%fcl.mfunctiony=fcl(x)y(1)=x(l)-0•7*sin(x(1))-0.2*cos(x(2))；y(2)=x(2
4、)-0.7*cos(x(l))+0.2*sin(x(2));y=[y(1)y(2)]；%lianxi44clcclear；close；sprintf('输入初始值:')x0=input('x0=')%[0.50.5];sprintf('原方程的解为：，)x=fsolve(*fcl',xO)1.求解一个常微分方程，写出M文件，应该包括函数式M文件和命令式M文件。functiondy=fc(x)dy=1Dy=-5*x1；closeclearclcs=dsolve(fc,1x1)2.求解一个常微分方程组，写出M文件，应该包括函数式M文件和命令式M文
5、件。%fc.mfunctiondy=fc(t,y)dy=[y(2);3*(l-y(l)x2)*y(2)-y(1)]；%lianxi.mclc;clear；close；[T/Y]=ode45('fc,,[020],[2;3]);plot(T,Y(:,1),心*,T,Y(:,2)宀」)title(9DE方程的解，)xlabel('T1)ylabel(1Y1)legend('Y1'z'Y21)3.输入一个矩阵，求出其行列式的值、逆矩阵、特征值、特征向量。%lianxi47clearclca=input「矩阵A—)%矩阵的输入必须带[]%a=[111
6、1;1234；13610;141020];det=det(a)%对应行列式的值inv=inv(a)%矩阵的逆[vzd]=eig(a)%矩阵的特征值和特征向量trace=trace(a)%矢巨卩车白勺迹rank=rank(a)%矢冃卩车的秩poly=poly(a)%特征多项式cond=cond(a)%条件数1.绘制一个二元函数的曲面图(例如y=sin(x+y)),加上标注说明，包扌舌坐标、抬头、曲面说明、网格、着色等。%lianxi48clear；close；clc；x=0:pi/30:pi*2；y=0:pi/30:pi*2；[X,Y]=mes
7、hgrid(x,y)；%将向量x,y指定的区域转化为矩阵X,Yz=sin(X+Y)；colormap(cool)；mesh(X,Y,z)%生成网格图xlabel(1x轴1)；ylabel('y轴*);zlabel('z轴')；title('sin(x+y)的曲线图T;legend('sin(x+y)曲面')；1.给出一个优化问题，写出M文件，进行求解。%lianxi491clear；close；clc；sprintf(，输入目标函数系数：T;f=input(1f=')；%f=[-5,4,2];sprintf(•输入系数矩阵：»);A=inpu
8、t(1A=')；%[6-11;124];sprintf(，输入右侧常数：
9、;b=input(1b=1)；%[8z10];vlb=[-1,0,0];vub=[3,2]
