南京邮电大学数学实验MATLAB2023

慈云数据 2024-03-20 技术支持 43 0

声明

使用R2023b版本完成,参考教材为数学实验(第三版)金正猛 王正新等编 科学出版社。参考了不少前人智慧的结晶,在此感谢,不保证全对,仅供大家参考,请勿抄袭,有错漏、疑问之处,烦请通过QQ或者评论区留言交流指正,欢迎大家多交流。联系QQ:3468039783

南京邮电大学数学实验MATLAB2023
(图片来源网络,侵删)

模块一:基础练习

1.1

Clear
syms x
m=628;
f=(log(1+x-m*x^2)-x)/(1-cos(x));
limit(f)
f=((sqrt(m*x^2+2)-atan(m*x))/x);
limit(f,x,inf)
pretty(limit(f,x,inf))

常用函数见书P5页,微积分部分见书P25页。pretty(f)将f显示为数学书写形式,见课本P13页 

1.2

clear
syms x
m=628;
f=exp(m*x)*sin(x);
diff(f,2)
subs(f,x,0);
diff(f,6)

1.3

clear
syms x
m=628;
f=(x+sin(x))/(1+cos(x));
int(f,x)
pretty(int(f,x))
f=log(1+m*x)-m*x;
int(f,x,0,1)
pretty(int(f,x,0,1))

1.4

clear
syms x
m=628;
f=cos(x*(m/200+sin(x)));
taylor(f,x,'Order',5,'ExpansionPoint',0)
pretty(taylor(f,x,'Order',5,'ExpansionPoint',0))

在学校的机房电脑上可能不可以加'Order',答案可能会变得很长,请自行取舍'Order','ExpansionPoint',0等

南京邮电大学数学实验MATLAB2023
(图片来源网络,侵删)

1.5

clear
x=[];
x(1)=rand(1);
m=628;
for k=2:10
    x(k)=sqrt(m/100+x(k-1));
end
x

rand函数见书P5页

1.6 

clear
m=628;
A=[4,2,m-2;-3,0,5;1,5,2*m]
B=[3,4,0;2,0,-3;-2,1,1]
det(A)
inv(A)
eig(A)
[P,D]=eig(A);
X=P(1:3,1:1)
Y=P(1:3,2:2)
Z=P(1:3,3:3)
A\B  
A/B
C=[A B];
Q=rref(C)

向量与矩阵见书P11、12页,线性代数见书P27页

1.7

1.7.1

f.m

function y=f(x)
if 0
微信扫一扫加客服

微信扫一扫加客服

点击启动AI问答
Draggable Icon