数学建模Matlab练习

发布 2023-05-17 18:45:28 阅读 7093

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元,...