西安交通大学matlab数学实验课件7.pptVIP

* * * * * * * * * * * * * * * * 1. 方程求根的Matlab指令 2. “二分法” 3. “牛顿切线法” 4. 一般迭代法 5. 分叉与混沌 非线性方程求解 方程近似根的求法 实验问题: 分析：该问题可以转化为下列方程求根. 开普勒方程 近似解. f=inline(x-0.5*sin(x)-1); c=fzero(f,[0,2]) 方法１－－－－利用MATLAB指令 输出：c = 1.4987 显然该方程为一元非线性方程，没有现成的求解公式. c=fzero(f,[a,b]) 求f在区间[a,b]的零点 c=fzero(f,x0)　　求f在x0附近的零点 方法２－－ “二分法” f=inline(x-0.5*sin(x)-1); fplot(f,[0,2]); grid 首先画图，观察出解的范围 然后利用二分法求解 二分法简介 分析： 条件： 得到第n个区间 终止条件： clear; f=inline(x-0.5*sin(x)-1); a=0;b=1; r=1.0e-5; N=log((b-a)/r)/log(2); k=0; while n =fix(N)+1 k=k+1;c=(a+b)/2; if abs(f(c))==0 break; elseif f(c)*f(b)0 a=c; else b=c; end fprintf(k=%d, x=%.5f\n,k,c); end 程序： 方法３－－－－迭代法　　　 （１）牛顿切线迭代法 分析： 条件： 关于牛顿切线法有下面收敛性定理 定理 clear; f=inline( x-0.5*sin(x)-1); df=inline(1-0.5*cos(x)); d2f=inline(0.5*sin(x)); a=1;b=2;dlt=1.0e-5; if f(a)*d2f(a)0 x0=a; else x0=b; end m=min(abs(df(a)),abs(df(b))); k=0; while abs(f(x0))m*dlt k=k+1; x1=x0-f(x0)/df(x0); x0=x1; fprintf(k=%d x=%.5f\n,k,x0); end 程序： （２） 一般迭代法 ——迭代格式 ——迭代序列 ——迭代初值 ——迭代函数 如何构造迭代函数 使得产生的迭代序列收敛? 例题 构造迭代格式 取初值 迭代序列 构造迭代格式 取初值 迭代序列 实验考察 如何构造迭代函数 定理 如果 在[a,b]上连续且满足 则迭代格式 对任意初值 迭代序列 收敛. 如果 (2) 存在正的常数h，对 用一般迭代法求解开普勒方程 令迭代函数为 显然 所以，迭代格式收敛 Logistic方程与混沌 在生物学中，有一个刻画生物种群个体总量增长情况的著名的方程——逻辑斯谛（Logistic）方程： 其中xn为某生物群体的第n代的个体总数与该群体所能达到的最大保有量时的个体数之比。 选定初值和比例系数r的值后，由方程就能生成一个数列： 考察迭代格式(Logistic方程 ) 初值 1. 当参数r取值分别为1.2，2.5，3.2，3.5，3.8 考察其迭代序列的收敛情况 clc;clf; x=0.1; y=[ ]; r=1.2; %改变取值得到相应的图形 hold on axis([0 100 0 1]) for i=1:100 x=r*x*(1-x);y=[y,x]; plot(i,x,k.,markersize,10) fprintf(x(%d)=%.10f\n,i,x); end t=1:100; plot(t,y,k-); grid 程序 2. 将参数r取0，0.3，0.6，0.9，1.2，…，3.9的迭代序列 收敛情况放置到同一坐标系中观察其变化 clear;clf; hold on axis([0,4,0,1]);grid for r=0:0.3:3.9 x=[0.1]; for i=2:150 x(i)=r*x(i-1)*(1-x(i-1)); end pause for i=101:150 plot(r,x(i),k.); end text(r-0.1,max(x(101:150))+0.05,[\it{r}=,num2str(r)])