数学建模Matlab练习4.docVIP

MATLAB练习四 1. 绘制两个一元函数的曲线图（例如y=sin(x)、y=cos(x)），加上标注说明，包括坐标、抬头、曲线说明等。 %lianxi41.m x=0:pi/30:2*pi; y=sin(x); plot(x,y,-ob,markersize,10); hold on y1=cos(x); plot(x,y1,-.dm,markersize,10); legend(sinx,cosx); title(sin和cos函数的曲线图) xlabel(x); ylabel(y); 2. 寻找一组数据拟合多项式，写出M文件，包括拟合的次数，数值的计算，数学表达式的呈现，图形的呈现等。 %lianxi42 x=0:pi/20:pi; y=sin(x); n=input(n=); p=polyfit(x,y,n); x1=0:pi/20:pi*2; y1=sin(x1); y2=polyval(p,x1) plot(x1,y1,-ob,x1,y2,:dm); legend(原曲线,拟合曲线); xlabel(x); ylabel(y); title(sin的拟合曲线和原曲线); axis([-0.2,6.5,-1.2,1.2]) 3. 求解一个代数方程，写出M文件，应该包括函数式M文件和命令式M文件。 %fc3.m function f=fc3(x) f=3*x-1=0; %lianxi43 clc; clear; close; sprintf(原方程和其解析解的结果为：) s=solve(fc3) sprintf(保留10位有效数字的结果为：) digits(10) sprintf(数值解的结果为：) x=vpa(s) 4. 求解一个代数方程组，写出M文件，应该包括函数式M文件和命令式M文件。 %fc1.m function y=fc1(x) y(1)=x(1)-0.7*sin(x(1))-0.2*cos(x(2)); y(2)=x(2)-0.7*cos(x(1))+0.2*sin(x(2)); y=[y(1) y(2)]; %lianxi44 clc clear; close; sprintf(输入初始值:) x0=input(x0=) %[0.5 0.5]; sprintf(原方程的解为：) x=fsolve(fc1,x0) 5. 求解一个常微分方程，写出M文件，应该包括函数式M文件和命令式M文件。 function dy=fc(x) dy=Dy=-5*x; close clear clc s=dsolve(fc,x) 6. 求解一个常微分方程组，写出M文件，应该包括函数式M文件和命令式M文件。 %fc.m function dy=fc(t,y) dy=[y(2);3*(1-y(1)^2)*y(2)-y(1)]; %lianxi.m clc; clear; close; [T,Y]=ode45(fc,[0 20],[2;3]); plot(T,Y(:,1),-,T,Y(:,2),--) title(ODE方程的解) xlabel(T) ylabel(Y) legend(Y1,Y2) 7. 输入一个矩阵，求出其行列式的值、逆矩阵、特征值、特征向量。 %lianxi47 clear clc a=input(矩阵A=)%矩阵的输入必须带[] %a=[1 1 1 1;1 2 3 4;1 3 6 10;1 4 10 20]; det=det(a)%对应行列式的值 inv=inv(a)%矩阵的逆 [v,d]=eig(a)%矩阵的特征值和特征向量 trace=trace(a)%矩阵的迹 rank=rank(a)%矩阵的秩 poly=poly(a)%特征多项式 cond=cond(a)%条件数 8. 绘制一个二元函数的曲面图（例如y=sin(x+y)），加上标注说明，包括坐标、抬头、曲面说明、网格、着色等。 %lianxi48 clear; close; clc; x=0:pi/30:pi*2; y=0:pi/30:pi*2; [X,Y]=meshgrid(x,y);%将向量x,y指定的区域转化为矩阵X,Y z=sin(X+Y); colormap(cool); mesh(X,Y,z)%生成网格图 xlabel(x轴); ylabel(y轴); zlabel(z轴); title(sin(x+y)的曲线图); legend(sin(x+y)曲面); 9. 给出一个优化问题，写出M文件，进行求解。 %lianxi491 clear; close; clc; sprintf(输入目标函数系数：); f=input(f=); %f=[-5,4,2]; sprintf(输入系数矩阵：); A=input(A=); %[6 -1 1;