第二次练习题。
1、 设,数列是否收敛?若收敛,其值为多少?精确到6位有效数字。
> f=inline('(x+7/x)/2');
syms x;
x0=3;for i=1:1:20
x0=f(x0);
fprintf('%g,%g',i,x0);end
本次计算运行到第三次结果稳定,可得:
数列收敛,收敛到2.64575
2、 设是否收敛?若收敛,其值为多少?精确到17位有效数字。
学号为单号,取。
> s=0;
for i=1:1:200
s=s+1/i^7;
fprintf('%g,%20.17f',i,s);end
运行至第190次后稳定,值为1.00834927738192070
书上习题:(实验四)
1,2,4,7(1),8,12(改为:对例2,取观察图形有什么变化.),13,14 。
练习1 编程判断函数的迭代序列是否收敛.
> f=inline('(x-1)/(x+1)')
x0=4;
for i=1:20
x0=f(x0);
fprintf('%g,%g',i,x0);end
由此可以发现迭代数列不一定收敛,迭代**现循环。
练习2 先分别求出分式线性函数、的不动点,再编程判断它们的迭代序列是否收敛.
运用上节的收敛定理可以证明:如果迭代函数在某不动点处具有连续导数且导数值介于-1与1之间,那末取该不动点附近的点为初值所得到的迭代序列一定收敛到该不动点.
1)解方程,得到x=-1,是函数f1(x)的不动点。
x=(x-1)/(x+3)
x =-1f1=inline('(x-1)/(x+3)')x0=-0.5;
for i=1:2000
x0=f1(x0);
fprintf('%g,%g',i,x0);end
2)解方程,得到x=-5和3,是函数f2(x)的不动点。
x=(-x+15)/(x+1)
x=-5,3;
format long;
f2=inline('(x-15)/(x+1)')x0=6;
for i=1:2000
x0=f2(x0);
fprintf('%g,%g',i,x0);end
由此可见由于迭代序列虽有不动点x=-1,但在此处导数不在-1与1之间,所以迭代数序列不收敛。
练习4 能否找到一个分式线性函数,使它产生的迭代序列收敛到给定的数?用这种办法近似计算.
x0=1;stopc=1;eps=10^(-17); a=1;b=2; c=1; d=1; k=0;
f=inline('(a*x+b)/(c*x+d)')
kmax=10;
while stopc>eps&k k=k+1;x0=f(a,b,c,d, x0);
stopc=abs(x0^3-2);
fprintf('%i,%g',k,x0)end
练习7 下列函数的迭代是否会产生混沌?
> f=inline('1-2*abs(x-1/2)')x=y=
x(1)=rand;
y(1)=0;x(2)=x(1);y(2)=f(x(1));
for i=1:10000
x(1+2*i)=y(2*i);
x(2+2*i)=x(1+2*i);
y(1+2*i)=x(1+2*i);
y(2+2*i)=f(x(2+2*i));
endplot(x,y,'r');
hold on;
syms x;
ezplot(x,[0,1]);
ezplot(f(x),[0,10]);
axis([0,1,0,1]);
hold off
练习8 函数=(01)称为logistic映射,试从“蜘蛛网”图观察它取初值为= 0.5产生的迭代序列的收敛性,将观察记录填人下表,若出现循环,请指出它的周期.
表4.3 logistic迭代的收敛性。
> f=inline('3.3*x*(1-x)')
> x=
> y=
> x(1)=0.5;
> y(1)=0;x(2)=x(1);y(2)=f(x(1));
> for i=1:10000
x(1+2*i)=y(2*i);
x(2+2*i)=x(1+2*i);
y(1+2*i)=x(1+2*i);
y(2+2*i)=f(x(2+2*i));
end> plot(x,y,'r');
> hold on;
> syms x;
> ezplot(x,[0,1]);
> f=inline('3.5*x*(1-x)')
> x=
y=x(1)=0.5;
y(1)=0;x(2)=x(1);y(2)=f(x(1));
for i=1:10000
x(1+2*i)=y(2*i);
x(2+2*i)=x(1+2*i);
y(1+2*i)=x(1+2*i);
y(2+2*i)=f(x(2+2*i));
endplot(x,y,'r');
hold on;
syms x;
ezplot(x,[0,1]);n
ezplot(f(x),[0,1]);
axis([0,1,0,1]);
hold off
> f=inline('3.56*x*(1-x)')
> x=
y=x(1)=0.5;
y(1)=0;x(2)=x(1);y(2)=f(x(1));
for i=1:10000
x(1+2*i)=y(2*i);
x(2+2*i)=x(1+2*i);
y(1+2*i)=x(1+2*i);
y(2+2*i)=f(x(2+2*i));
endplot(x,y,'r');
hold on;
syms x;
ezplot(x,[0,1]);
ezplot(f(x),[0,1]);
axis([0,1,0,1]);
hold off
> f=inline('3.568*x*(1-x)')
> x=
y=x(1)=0.5;
南邮2019数学建模校赛A题
a题。计划生育政策调整对人口数量 结构及其影响的研究。人口的数量和结构是影响经济社会发展的重要因素。从20世纪70年代后期以来,我国鼓励晚婚晚育,提倡一对夫妻生育一个孩子。该政策实施30多年来,有效地控制了我国人口的过快增长,对经济发展和人民生活的改善做出了积极的贡献。但另一方面,其负面影响也开始显...
南邮概率统计模拟试卷
南京邮电大学 2010 2011学年第 2 学期。概率论和数理统计 模拟试卷。院 系班级学号姓名。特别提醒 备用数据。一 填空题 共48分,每格3分 1 事件a,b,c中恰好有一个发生可表示为。2.已知,3.袋中有20个乒乓球,其中12个是黄球,8个是白球,今有两人依次随机地从袋中各取一球,取后不放...
南邮形势与政策作业世博会
上海世博会。姓名学号 世博会是人类科学的盛会,在本质上,是人文精神的科技化 物化和科学技术的人文化。世博会的创新促进了科技文化知识的传播和普及。现代世博会更加注重文明的对话,文化的交流。可以说文化贯串于世博会的始终,而世博会的文化影响力也是人类文明进步的助推器。世博精神是人文精神的延续。世博会是技术...