用MATLAB编程完成一道数学题,要求如下

假设如下(1)由于自然繁殖,兔子数每年增长10%,狐狸数每年减少15%(2)由于狐狸吃兔子,兔子数每年减少狐狸数的0.15倍,狐狸数每年增加兔子数的0.1倍(3)第n年时,兔子数用Rn表示,狐狸数用Fn表示(4)初始时刻(既第0年)兔子数R0=100只,狐狸数有F0=30只,则
1.求Rn,Fn的关系式
2. 求Rn,Fn关于n的关系式
3.讨论当n越来越大时,兔子与狐狸的数量能否达到一个稳定的平衡状态说明理由。(画图表示Rn,Fn的变化趋势)

Rn=1.1Rn-1-0.15Fn-1

Fn=0.1Rn-1+0.85Fn-1

写成矩阵形式: Xn=A*Xn-1, Xn=[Rn Fn]', A=[1.1 -.15; 0.1 0.85]

根据差分方程组解法 Xn=A^nX0


A=[1.1 -0.15;0.1 0.85];

x(1,1)=100;

x(2,1)=30;

N=100;

for k=1:N

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

end

plot(x(1,:),'r-');

hold on;

plot(x(2,:),'b-');

legend('Rabbit','Fox');

grid on;



求Rn,Fn关于n的关系式:

求A的特征根 λ1=1,λ2=0.95,  Rn=a1λ1^n+b1λ2^n, Fn=a2λ1^n+b2λ2^n

再根据初始条件 可得 

Rn=210-110*0.95^n 

Fn=140-110*0.95^n

n->∞时, 0.95^n->0, Rn->210, Fn->140

追问

运行不了 具体的程序怎么写 可以截图给我吗 谢谢啦

追答

可以运行啊,出错信息是什么?

温馨提示:内容为网友见解,仅供参考
无其他回答

用MATLAB编程完成一道数学题,要求如下
Fn=0.1Rn-1+0.85Fn-1 写成矩阵形式: Xn=A*Xn-1, Xn=[Rn Fn]', A=[1.1 -.15; 0.1 0.85]根据差分方程组解法 Xn=A^nX0 A=[1.1 -0.15;0.1 0.85];x(1,1)=100;x(2,1)=30;N=100;for k=1:N x(:,k+1)=A*x(:,k);end plot(x(1,:),'r-');hold on;...

用MATLAB编程做一道数学题 需要具体的编写内容 求大神协助
R0=100;%兔子 F0=30;%狐狸 year=1;R(1) = (1+0.1)*R0-0.15*F0;F(1) = (1-0.15)*F0+0.1*R0;iter =1e-5;while 1 R(year+1) = (1+0.1)*R(year)-0.15*F(year);F(year+1) = (1-0.15)*F(year)+0.1*R(year);if (R(year+1)-R(year))<iter if (F(...

用matlab编写一个程序,要求如下,在一个excel表格中有两个子表格,sheet1...
按以下步骤做:ctrl+f3 在弹出的窗口的名称位置输入aa,引用位置输入=mid(get.workbook(1),find("[",get.workbook(1))+1,99)&t(now())添加,确定 在任意空白单元格输入=counta(aa)即可 如果还要更好的,hi我,给你一个统计表的个数并能输出各个表名称的文件 ...

一道Matlab编程题
1. 首先举一个简单的例子:求y=x^2 ,在x为[0,2]上的曲线长度。把下面的复制粘贴进MATLAB syms t x=t;y=t^2;df=@(t)(1+4*t.^2).^0.5; %%MATLAB早期版本不支持@功能 quad(df,0,1)答案ans=1.4789 2. 再回答你的问题:clc clear syms t x=sin (t);y=t^2;z=log(t...

一道需要matlab来实现的小问题,有关数值计算,can you help me?_百度知...
这个有点爱莫能助了...不好意思...数学方面的matlab不是很懂,等等别人答复吧~~

如何用matlab画一个微分方程组的图,具体要求如下
用matlab求解微分方程组可以用ode()函数求得其数值解。求解方法:1、首先,自定义微分方程组函数odefun(t,x),其主要内容 dy(1)=x(1)*(1-x(1))-x(2)+1\/x(1);dy(2)=x(2)*(-1+x(1));2、再用ode()函数求得其数值解。使用方法 [t,x]=ode45(@odefun,[0 20],y0);3、用...

用matlab编写一个程序
function subfun(a,b)x=a:0.01:b; %a,b为你所要的自变量取值范围 y=(x.^3+5*x-10)\/6.*(x+5<0)+exp(x-1).*(-5<=x<5)+log(7*x+4).*(x>=5);plot(x,y)保存后,在command window下输入 subfun(-10,10)就可以了……分段函数都可以引入一个向量,自变量在某个范围...

用matlab解一元二次方程(输入矩阵,输出矩阵)
SA阴影=(14-2y)3x=42x-6xy(1)求图1中12个小长方形xy面积Sxy=xy,根据题意,12Sxy+SA阴影+SB阴影=14x10 12xy+42x-6xy+SB阴影=140 SB阴影=140-6xy-42x(2)根据题意,SA阴影-SB阴影=112 (3)(1)(2)代入(3)得,42x-6xy-(140-6xy-42x)=112 42x-6xy-140+...

matlab数学实验的一道求极值题,怎么玩,我没下那个软件
第5题的极值,可以用fminbnd()极值函数,即可解决。clc fun=inline('(log(x))^2\/x','x'); %自定义函数 [x,fval,exitflag]=fminbnd(fun,0.5,1.5) %极值点

跪求数学实验几何物理插值问题的一道题,要求步骤详细,用MatLab...
以左边零点位坐标原点,建立坐标系。线性插值得到图形,再用数值积分可求面积 x=linspace(0,8.534,13)y=[0 0.914 5.060 7.772 8.717 9.083 9.144 9.083 8.992 8.687 7.376 2.073 0];x0=0:0.001:8.534;y1=interp1(x,y,x0);figure,plot(x,y,'k*',x0,y1,'-r')S=...

相似回答