用Matlab求解差分方程问题市公开课一等奖百校联赛特等奖课件.pptxVIP

用Matlab求解差分方程问题;一、一阶线性常系数差分方程;模型建立;Matlab实现;;利用plot绘图观察数量改变趋势;plot(k,y1,k,y2,k,y3)在同一坐标系下画图

;;;k=(0:20);%一个行向量

y1=(20,-0.0324,5);%也是一个行向量

round([k’,y1’])%对k,y1四舍五入，但

%是不改变变量值

plot(k,y1)%ky1是行向量列向量都能够

也能够观察发展趋势，以及在较差条件下发展趋势，也能够考查每年孵化数量改变影响。;高阶线性常系数差分方程;一年生植物繁殖;模型及其求解;Xk=a1bcXk-1+a2b(1-a1)bcXk-2;Xk=a1bcXk-1+a2b(1-a1)bcXk-2;;结果分析：Xk=pXk-1+qXk-2(1)

x1+px0=0(2)

;本例中，用待定系数方法能够求出

b=0.18时,

c1=95.64,c2=4.36，

这么

实际上，

植物能一直繁殖下去条件是b0.191;线性常系数差分方程组;;模型及其求解;用矩阵表示

用matlab编程，计算x(k),观察n年以后3个城市汽车数量改变情况

;functionx=czqc(n)

A=[0.6,0.2,0.1;0.3,0.7,0.3;0.1,0.1,0.6];

x(:,1)=[200,200,200];

fork=1:n

x(:,k+1)=A*x(:,k);

end

假如直接看或者发展趋势，能够直接在命令窗

口（commondwindow）作，而不是必须编一个函数;A=[0.6,0.2,0.1;0.3,0.7,0.3;0.1,0.1,0.6];

n=10;

fork=1:n

x(:,1)=[200,200,200];

x(:,k+1)=A*x(:,k);

end

round(x);作图观察数量改变趋势;;直接输入x（:,1）值即可;6.6按年纪分组人口模型;假设

与

建模;按年纪分组种群增加;模型及其求解;记在时段k种群各年纪组数量为

X(k)=[x1(k),x2(k),····,xn(k)]’;这么，有x(k+1)=Lx(k),k=0,1,····

给定在0时段，各年纪组初始数量x(0)

就能够预测任意时段k,各年纪组数量

设一个群分成5个年纪组，

繁殖率b1=0,b2=0.2,,b3=1.8,b4=0.8,b5=0.2

存活率s1=0.5,s2=0.8,s3=0.8,s4=0.1

各年纪组现有数量都是100只，

用matlab计算x(k);b=[0,0.2,1.8,0.8,0.2];

s=diag([0.5,0.8,0.8,0.1]);

L=[b;s,zeros(4,1)];

x(:,1)=100*ones(5,1);

n=30;

fork=1:n

x(:,k+1)=L*x(:,k);

end

round(x)

k=0:30;

subplot(1,2,1),plot(k,x),grid

;将x(k)归一化后向量记做x’(k),称为种群按年纪组分布向量，即各年纪组在k时段在数量上占总数百分比。

y=diag(1./sum(x));%sum(x)对列求和

Z=x*y

Subplot(1,2,2),plot(k,z),grid

结果分析：时间充分长以后，种群按年纪组分布x’(k)趋向稳定。