南京邮电大学数学实验练习题参考答案 下载本文

专业 姓名 学号 成绩

当?=3.6时,得:

-(18 x (x - 1))/510.90.80.70.60.50.40.30.20.1000.10.20.30.40.5x0.60.70.80.91

当?=3.84时,得:

-(96 x (x - 1))/2510.90.80.70.60.50.40.30.20.1000.10.20.30.40.5x0.60.70.80.91

数学实验实验报告

17

专业 姓名 学号 成绩

2.5 对于Martin迭代,取参数a,b,c为其它的值会得到什么图形?参考下表(取自63页练习13)

a m -m -m m/1000 m/1000 m/100 -m/10 b m -m m/1000 m/1000 m m/10 17 c m m -m 0.5 -m -10 4

解:取m=10000;迭代次数N=20000; 在M-文件里面输入代码:

function Martin(a,b,c,N)

f=@(x,y)(y-sign(x)*sqrt(abs(b*x-c))); g=@(x)(a-x); m=[0;0]; for n=1:N

m(:,n+1)=[f(m(1,n),m(2,n)),g(m(1,n))]; end

plot(m(1,:),m(2,:),'kx'); axis equal

在命令窗口中执行Martin(10000,10000,10000,20000),得:

数学实验实验报告

18

专业 姓名 学号 成绩

2500020000150001000050000-5000-2-1.5-1-0.500.511.52x 104

执行Martin(-10000,-10000,10000,20000),得:

50000-5000-10000-15000-20000-25000-2-1.5-1-0.500.511.52x 104

数学实验实验报告

19

专业 姓名 学号 成绩

执行Martin(-10000,10,-10000,20000),得:

0-2000-4000-6000-8000-10000-12000-10000-8000-6000-4000-200002000

执行Martin(10,10,0.5,20000),得:

302520151050-5-10-20-100102030

数学实验实验报告

20