matlab练习四。
1. 绘制两个一元函数的曲线图(例如y=sin(x)、y=cos(x)),加上标注说明,包括坐标、抬头、曲线说明等。
x=0:pi/30:2*pi;
y=sin(x);
plot(x,y,'-ob','markersize',10);
hold on
y1=cos(x);
plot(x,y1,''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文件。
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文件。
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
clcclear;
close;
sprintf('输入初始值:')
x0=input('x0=')
sprintf('原方程的解为:')
x=fsolve('fc1',x0)
5. 求解一个常微分方程,写出m文件,应该包括函数式m文件和命令式m文件。
function dy=fc(x)
dy='dy=-5*x';
closeclear
clcs=dsolve(fc,'x')
6. 求解一个常微分方程组,写出m文件,应该包括函数式m文件和命令式m文件。
function dy=fc(t,y)
dy=[y(2);3*(1-y(1)^2)*y(2)-y(1)];
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
clearclc
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=')
sprintf('输入右侧常数:')
b=input('b=')
vlb=[-1,0,0];
vub=[3,2];
x=linprog(f,a,b,vlb,vub)
10. 随意编写一个m文件,解决一个具体问题。
同第六题。
Matlab与数学建模综合练习
1 求的所有根。先画图后求解 2 求下列方程的根。3 求解下列各题 4 1 求矩阵的逆矩阵及特征值和特征向量。2 求点 1,1,4 到直线l x 3 1 y 0 z 1 2的距离。5 已知分别在下列条件下画出的图形 6 画下列函数的图形 1 7 设,数列是否收敛?若收敛,其值为多少?精确到6位有效数...
Matlab与数学建模综合练习
1 求的所有根。先画图后求解 2 求下列方程的根。3 求解下列各题 4 1 求矩阵的逆矩阵及特征值和特征向量。2 求点 1,1,4 到直线l x 3 1 y 0 z 1 2的距离。5 已知分别在下列条件下画出的图形 6 画下列函数的图形 1 7 设,数列是否收敛?若收敛,其值为多少?精确到6位有效数...
Matlab与数学建模综合l练习
1 某厂有一台制杯机,可生产两种型号的杯子,a型杯每6小时可生产100箱,b型杯每5小时可生产100箱,这台机器每周生产时间为60小时,生产出的产品堆放在仓库里,库容量为15000立方米,a型杯每箱占有空间10立方米,b型杯每箱占有空间20立方米,生产a型杯每箱可获利5元,b型杯每箱可获利4.5元,...