北科大matlab 第四次作业.docVIP

下载本文档

13
0
约3.53千字
约 10页
2018-08-13 发布于贵州
举报
版权申诉

北科大matlab 第四次作业.doc

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 第四次作业

PAGE PAGE 2 《数学实验》报告实验名称?? 线性代数相关运算及数值方法计算定积分学院专业班级姓??? 名???????? ??? 学??? 号??????? ??? ?? 2012年5月【实验目的】掌握矩阵的基本运算、特征值、特征向量和线性方程组的求解；能熟练运用数值方法求定积分【实验任务】 P114 12,14 P115 21. (1),(2) 第七章习题P167 17. (2) 18 【实验程序】 1..A=[1 0 2 0 0 5;3 5 6 4 3 4;2 3 1 4 6 7;0 7 6 4 5 3;1 4 5 7 4 5;1 6 7 3 4 2] A det(A) rank(A) rref(A) 2. A=[2 1 1;1 2 1;1 1 2] poly(A) [V D]=eig(A) 3. A=[1 1 2 -1;-1 1 3 0;2 -3 4 -1] rref(A) 4. A=[1 -1 -1 1 0;1 -1 1 -3 1;1 -1 -2 3 -1/2] rref(A) 5.function y=jifen(x) y=x.*sin(x)./(1+(cos(x)).^2) x=0:0.01:pi; y=x.*sin(x)./(1+(cos(x)).^2); t=length(x); s1=sum(y(1:(t-1)))*0.01 s2=sum(y(2:t))*0.01 s3=trapz(x,y) s4=quad(jifen,0,pi) 6. function y = jifen(x) y=1./(1-sin(x)); x=0:0.01:pi/4; y=1./(1-sin(x)); t=length(x); s1=sum(y(1:(t-1)))*0.01 s2=sum(y(2:t))*0.01 s3=trapz(x,y) s4=quad(jifen,0,pi/4) 【实验结果】 1. A=[1 0 2 0 0 5;3 5 6 4 3 4;2 3 1 4 6 7;0 7 6 4 5 3;1 4 5 7 4 5;1 6 7 3 4 2] A = 1 0 2 0 0 5 3 5 6 4 3 4 2 3 1 4 6 7 0 7 6 4 5 3 1 4 5 7 4 5 1 6 7 3 4 2 A ans = 1 3 2 0 1 1 0 5 3 7 4 6 2 6 1 6 5 7 0 4 4 4 7 3 0 3 6 5 4 4 5 4 7 3 5 2 det(A) ans = -2878 rank(A) ans = 6 rref(A) ans = 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 0 0 0 0 0 0 1 2. A=[2 1 1;1 2 1;1 1 2] A = 2 1 1 1 2 1 1 1 2 poly(A) ans = 1.0000 -6.0000 9.0000 -4.0000 [V D]=eig(A) V = 0.4082 0.7071 0.5774 0.4082 -0.7071 0.5774 -0.8165 0 0.5774 D = 1.0000 0 0