- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
中央财经大学MATLAB第4讲课余作业
《第4讲 符号计算》课余作业
1、创建符号表达式f(x)=sin(x)+cos(x)-tan(x)
syms x
f= sin(x)+cos(x)-tan(x)
f =
cos(x) + sin(x) - tan(x)
2、计算习题1中的符号表达式在x=0, x=1, x=2*pi处的值
syms x
f=sym(sin(x)+cos(x)-tan(x));
subs(f,x,[0 1 2*pi])
ans =
1.0000 -0.1756 1.0000
3、计算符号表达式f(x)=x+cos(x)-sin(x)在x=pi-1处的值,并将结果设置为以下5中精度,即分别为小数点后1位、2位、10位、20位和50位有效数字。
syms x
f=x+cos(x)-sin(x);
a=subs(f,x,pi-1);
a1=vpa(a,1) % 在Matlab2009a版本中,此程序无效;在Matlab6.5中可用
fprintf(a1=\n%3.1f\n,a) % 输出一位有效数字
a2=vpa(a,2)
a10=vpa(a,10)
a20=vpa(a,20)
a50=vpa(a,50)
a1 =
0.75981936291375673508952104384662
a1=
0.8
a2 =
0.76
a10 =
0.7598193629
a20 =
0.75981936291375673509
a50 =
0.75981936291375673508952104384661652147769927978516
4、设x为符号变量,f(x)=x^4+2*x^2+1, g(x)=x^3+6*x^2+3*x+5,试进行如下运算:
(1)f(x)+g(x)
(2)f(x)*g(x)
(3)对f(x)进行因式分解
(4)求g(x)的反函数
syms x
f=x^4+2*x^2+1;
g=x^3+6*x^2+3*x+5;
f+g
f*g
horner(f)
finverse(g)
ans =
x^4 + x^3 + 8*x^2 + 3*x + 6
ans =
(x^4 + 2*x^2 + 1)*(x^3 + 6*x^2 + 3*x + 5)
ans =
x^2*(x^2 + 2) + 1
Warning: finverse(x^3 + 6*x^2 + 3*x + 5) is not unique.
In sym.finverse at 46
ans =
3/(x/2 + ((x/2 - 15/2)^2 - 27)^(1/2) - 15/2)^(1/3) + (x/2 + ((x/2 - 15/2)^2 - 27)^(1/2) - 15/2)^(1/3) - 2
5、设x为符号变量,令y=x^4+2*x^2+1,z=y^2,w=sin(z),试求w关于x的符号表达式,并求当x=5时,w的小数点之后具有1位、2位、10位、20位和50位有效数字的数值解。
syms x
y=x^4+2*x^2+1;
z=y^2;
w=sin(z)
w0=subs(w,x,5)
fprintf(w1=\n%3.1f\n,w0)
vpa(w0,2)
vpa(w0,10)
vpa(w0,20)
vpa(w0,50)
w =
sin((x^4 + 2*x^2 + 1)^2)
w0 =
-0.0673
w1=
-0.1
ans =
-0.067
ans =
-0.06734017262
ans =
-0.067340172619870686255
ans =
-0.067340172619870686254728298081317916512489318847656
6、使用sym函数生成如下符号矩阵, a=sym([1/x,1/(x+1);1/(x+2) 1/(x+3)]), b=sym([x, 1; x+2,0]),试对着两个符号矩阵进行如下操作:
(1)a-b
(2)a*b
(3)求a的倒置
(4)求a的行列式
(5)求b的逆
(6)求a的秩
(7)求a^3
a=sym([1/x,1/(x+1);1/(x+2) 1/(x+3)]);
b=sym([x, 1; x+2,0]);
a-b
a*b
a % 求a的倒置
det(a) % 求a的行列式
inv(b) % 求b的逆
rank(a) % 求a的秩
a^3
ans =
[ 1/x - x, 1/(x + 1) - 1]
[ 1/(x + 2) - x - 2, 1/(x + 3)]
ans =
[ (x + 2)/(x + 1) + 1, 1
您可能关注的文档
最近下载
- 横河涡街流量计第10版DY中文说明书.pdf VIP
- 混凝土挡土墙施工方案..doc VIP
- 六年级梦想之星英语初赛试卷含答案.doc VIP
- 智慧养老机构运营与管理:智慧养老机构健康服务管理PPT教学课件.pptx VIP
- 语文练习部分六年级第一学期(答案).pdf VIP
- 新22G02 钢筋混凝土结构构造 .docx VIP
- 《应急管理》ppt课件.pptx VIP
- 智慧养老机构运营与管理:智慧养老机构生活和快乐服务PPT教学课件.pptx VIP
- 商务部对外贸易司-汽车行业:2022中国二手车出口国别指南.docx VIP
- 智慧养老机构运营与管理:智慧养老机构出入院管理PPT教学课件.pptx VIP
文档评论(0)