南京邮电大学数学实验MATLAB2023

声明

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

模块一：基础练习

1.1

clearsyms xm=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

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

1.3

clearsyms xm=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

clearsyms xm=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等

1.5

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

rand函数见书P5页

1.6 

clearm=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/BC=[A B];Q=rref(C)

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

1.7

1.7.1

f.m

function y=f(x)if 0