专业 姓名 学号 成绩

当?=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