数值分析几十个程序.pdfVIP

  1. 1、本文档共25页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多

第一题关于舍入误差累计的效果模拟。x*是在(0,1)上服从均匀

分布的随机数,对x*取5位有效数字得到x,将产生的k(比如取为

10000)个x*(x)相加得到X*(X),研究X*-X的分布情况以及X*-X

和k的关系。将得到的结果用图形表示出来。

x10-5

12

10

8

6

4

2

0

-2

-4

0102030405060708090100

x10-12

2.5

2

1.5

1

0.5

差0

-0.5

-1

-1.5

-2

-2.5

01002003004005006007008009001000

K值

k=1;

kk=ones(1,1000);

xc=ones(1,1000);

whilek=1000

r=rand(1,k);

v=vpa(r,5);

x0=sum(r);

x1=sum(v);

cha=x0-x1;

xc(k)=cha;

kk(k)=k;

k=k+1;

end;

plot(kk,xc,rh)

可以看出随着处理的数据的增大误差也越来越大,但是还分布在横轴的两侧

第二题研究产生各种特定矩阵的方法(阶数在10-100),比如对称

阵,三对角阵,正定矩阵,正交阵,对角占优矩阵,说明如何生成。

(1)生成对称矩阵

k=10;

a=rand(10,10);

fori=1:k%先生成一个下三角矩阵;

forj=1:k

ifji

a(i,j)=0;

end;

end;

end;

aa=a;%加上他本身的转置;

g=a+aa

g=

1.97600.03770.88520.91330.79620.0987

0.26190.33540.67970.1366

0.03770.21350.65380.49420.77910.7150

0.90370.89090.33420.6987

0.88520.65381.48810.50000.47990.9047

0.60990.61770.85940.8055

0.91330.49420.50001.77300.02870.4899

0.16790.97870.71270.5005

0.79620.77910.47990.02870.14290.5216

0.09670.81810.81750.7224

0.09870.71500.90470.48990.52161

文档评论(0)

186****1451 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档