MATLAB符号运算（7）

实验目的

1）掌握定义符号对象和创建符号表达式的方法；
2）掌握符号运算基本命令和规则；
3）掌握符号表达式的运算法则以及符号矩阵运算；
4）掌握符号微积分和符号方程求解的基本方法。

实验要求

1）实验前，要认真预习实验任务，了解实验目的和实验内容；
2）实验时，要利用MATLAB语言编写程序代码形成独立的M文件，并调试程序使其能正确运行；
3）实验后，按要求编写实验报告，源程序要有适当的注释，以提高程序的可读性。

f1=sym('cos(x)+sqrt(1-sin(x)^2)') 
syms y t 
f2=y/exp(-2*t) 
%即sym可以在创建变量同时赋值，赋值后仍然为sym型；而syms不能在创建同时赋值
%syms是定义符号变量sym是将字符或者数字转换为字符

syms x f;
x=1.723;
f=(3*x^3+x^2-1)/(x^2+1);
t=vpa(subs(f,x,1.723),20) %vpa(S,d) S数值型符号对象，d是精度

syms a n x; %symsum求和
s1=a*x^a;
s2=a^2;
t1=symsum(s1,a,1,n)   %1到n这个区间，对s1符号表达式中的a
t2=symsum(s2,a,1,100)

在这里插入代码片syms x y
f1=(exp(x)-exp(y))/(x-y);
y1=limit(f1,x,y) %limit(f,x,a)求函数f中x到a的极限
f2=tan(x)^tan(2*x);
y2=limit(f2,x,pi/4,'right')  %right pi/4^+;left pi/4^- 

这里是引用

%compose(f,g,z)是f(g(z))
syms x v u
y=sqrt(1+u^2);
u=log(v);
v=exp(-x);
compose(compose(y,u,x),v,x) 

syms x y;
z=x^9+7*y^4-x^5*y^3; 
y1=diff(z,x,2)
y2=diff(z,y,2)  %diff(a,b,c) 对这a中含b的自变量求c阶导
y3=diff(diff(z,x),y)

syms x y z
s1=2*x-y+3*z-5;
s2=3*x+y-5*z-5;
s3=4*x-y+z-1;
[x,y,z]=solve(s1,s2,s3)
%调用solve函数：使用solve函数求解方程或方程组的符号解，例如：sol =solve(S1，S2，...Sn);

syms x
f1=x^4/(1+x^2);
f2=3*x^2-x+1;
t1=int(f1,x) %int(S,v)以v为自变量求f1的不定积分
t2=int(f2,x,0,2) 
%int(S,v，a,b)以v为自变量求f1的在（a,b)定积分

syms x y z;
[y,z]=dsolve('Dy-z=sin(x)','Dz+y=1+x','y(0)=2','z(0)=7','x')
%符号微分方程求解函数dsolve的调用格式：
% w=dsolve(|eq1,eq2,…′,cond1,cond2,…′,'v′);
% 对常微分方程(组)eq1，eq2，…中指定的符号自变量v，
% 求符号w，cond2为条件若没有指定变量v，则默认变量

实验结论

第1题知道sym与syms区别。第2题掌握求精确度，第3456题都是函数使用。第7题要掌握solve。第89题也是函数的使用。记得当时在做第7题的时候，我看资料没找到相关内容，然后就空在那里。我然后问小组的其他成员，他们有找到的，团队里写作业能更加更快，他在网上找到的。虽然有些实验我能单独完成，但有时候会遇到一两个不会的，恰恰这个时候就体现了团队的重要性。学习了本章，提高了程序的编写能力和加强了逻辑思维的培养。