y(2+2*i)=f(x(2+2*i)); end
plot(x,y,'r'); hold on; syms x;
ezplot(x,[0,1]); ezplot(f(x),[0,1]); axis([0,1,0,1]); hold off
练习12 对例2,对例2,取 a?4.5,25,55,120 观察图形有什么变化.试着提高迭代次数至26 000、28 000、100 000、500 000等观察图形
有什么变化.
>> Martin(45,2,-300,5000);
>> Martin(4.5,2,-300,5000);
>> Martin(25,2,-300,5000);
>> Martin(55,2,-300,5000);
>> Martin(120,2,-300,5000);
练习13 取参数a,b,c为其他的值会得到什么图形?参考表4.4.
表4.4 Martin迭代参数参考表
a b c -1000.1 -10 0 0.4 90 10 1 30 0 10 -10 100 -200 -4 -80 -137 17 10 4 100 -10
>> Martin(-1000,0.1,-10,5000);
>> Martin(-0.4,1,0,5000);
>> Martin(90,30,10,5000);