matlab第七次作业
Matlab上机作业部分参考答案.ppt
[rank(A), rank([A B])]
ans =
34 由得出的结果看，A, [A;B] 两个矩阵的秩不同，故方程是
矛盾方程，没有解。
5. 试求下面齐次方程的基础解系
7. 建立如下一个元胞数组，现在要求计算第一个元胞第4行第 2列加上第二个元胞+第三个元胞里的第二个元素+最后一个元 胞的第二个元素。
a={pascal(4),'hello';17.3500,7:2:100}
解： >> a={pascal(4),'hello';17.3500,7:2:100} a=
[ 173/34, 151/34]
6. 求解方程组的通解
x1 2x2 4x3 6x4 3x5 2x6 4 2x1 4x2 4x3 5x4 x5 5x6 3
3x1 6x2 2x3 5x5 9x6 1 2x1 3x2 4x4 x6 8
4x2
5x3
2x4
x5
参考答案： （1） >> limit(sym('(tan(x) - sin(x))/(1cos(2*x))')) ans = 0 （2） >> y = sym('x^3 - 2*x^2 + sin(x)'); >> diff(y) ans = 3*x^2-4*x+cos(x) （3） >> f = x*y*log(x+y); >> fx = diff(f,x) fx = y*log(x+y)+x*y/(x+y)
MATLAB习题及参考答案经典.doc
习题：1， 计算⎥⎦⎤⎢⎣⎡=572396a 与⎥⎦⎤⎢⎣⎡=864142b 的数组乘积。
2， 对于B AX =，如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ，求解X 。
3， 已知：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。
4， 角度[]604530=x ，求x 的正弦、余弦、正切和余切。
(应用sin,cos,tan.cot)5， 将矩阵⎥⎦⎤⎢⎣⎡=7524a 、⎥⎦⎤⎢⎣⎡=3817b 和⎥⎦⎤⎢⎣⎡=2695c 组合成两个新矩阵： （1）组合成一个4⨯3的矩阵，第一列为按列顺序排列的a 矩阵元素，第二列为按列顺序排列的b 矩阵元素，第三列为按列顺序排列的c 矩阵元素，即 ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡237912685574（2）按照a 、b 、c 的列顺序组合成一个行矢量，即 []2965318772546， 将(x -6)(x -3)(x -8)展开为系数多项式的形式。
(应用poly,polyvalm)7， 求解多项式x 3-7x 2+2x +40的根。
(应用roots)8， 求解在x =8时多项式(x -1)(x -2) (x -3)(x -4)的值。
(应用poly,polyvalm)9， 计算多项式9514124234++--x x x x 的微分和积分。
(应用polyder,polyint ，poly2sym)10， 解方程组⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡66136221143092x 。
(应用x=a\b)11， 求欠定方程组⎥⎦⎤⎢⎣⎡=⎥⎦⎤⎢⎣⎡5865394742x 的最小范数解。
(应用pinv) 12， 矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=943457624a ，计算a 的行列式和逆矩阵。
(应用det,inv)13， y =sin(x )，x 从0到2π，∆x =0.02π，求y 的最大值、最小值、均值和标准差。
第七次实验作业1
实验报告（二）院（系）理学院课程名称：数学实验日期班级B1111 学号12 实验室文理楼209 专业数学教育姓名樊美林计算机号实验名称数学实验成绩评定所用软件Matlab 教师签名实验目的或要求1.实验步骤、心得体会1. 求下列函数的极小点：1) ()2123222118294xxxxxXf+-++=；、function f = fun1(x)f=x(1)^2+4*x(2)^2+9*x(3)^2-2*x(1)+18*x(2);x0=[1,1,2]x=fminunc('fun1',x0);y=fun1(x)y=-21.25002)()212122212223xxxxxxXf-+-+=；function f = fun1(x)f=x(1)^2+(3/2)*x(2)^2+2*x(1)*x(2)+x(1)-2*x(2);x0=[1,1]x=fminunc('fun1',x0)y=fun1(x)y =-4.75003)()()42121f X x x=-+.function f = fun2(x)f=(x(1)-1)^4+x(2)^2；x0=[0,1]x=fminunc('fun2',x0);y=fun2(x)y =4.0848e-010第1），2）题的初始点可任意选取，第3）题的初始点取为()TX1,00=.2. 梯子长度问题一楼房的后面是一个很大的花园. 在花园中紧靠着楼房有一个温室,温室伸入花园2m,高3m,温室正上方是楼房的窗台. 清洁工打扫窗台周围,他得用梯子越过温室,一头放在花园中,一头靠在楼房的墙上. 因为温室是不能承受梯子压力的,所以梯子太短是不行的.现清洁工只有一架7m 长的梯子,你认为它能达到要求吗? 能满足要求的梯子的最小长度为多少?问题分析：（对问题作出分析，如果问题过于简单，可以不作分析） 模型建立：（建立数学模型） 结果：（求解的结果）分析：（对求解结果进行分析）其他：（模型的缺点和优点以及改进方向，如果没有，可以不写）模型的求解：（求解所用的软件，求解该问题的源代码）3. 陈酒出售的最佳时机问题某酒厂有批新酿的好酒,如果现在就出售,可得总收入R 0=50万元(人民币),如果窖藏起来待来日(第n 年)按陈酒价格出售,第n 年末可得总收入60n eR R (万元),而银行利率为r=0.05,试分析这批好酒窖藏多少年后出售可使总收入的现值最大. (假设现有资金X 万元,将其存入银行,到第n 年时增值为R(n)万元,则称X 为R(n)的现值.)并填下表：第一种方案：将酒现在出售，所获50万元本金存入银行； 第二种方案：将酒窖藏起来，待第n 年出售。
MATLAB作业7参考答案
MATLAB作业7 参考答案1、试生成满足正态分布2N的30000 个伪随机数，对其均值和方差进行验(0.5,1.4)证，并用直方图的方式观察其分布与理论值是否吻合，若改变直方图区间的宽度会得出什么结论？解:用下面的语句可以生成随机数，并计算均值和方差，可见，其结果接近给定的数值。
>> x=normrnd(0.5,1.4,30000,1);>>m=mean(x), s=std(x)m =0.49974242123102s =1.40033494141044>> xx=-5:0.3:5; yy=hist(x,xx); bar(xx,yy/length(x)/0.3); hold onx0=-5:0.1:5; y0=normpdf(x0,0.5,1.5); plot(x0,y0)>> xx=-5:0.8:5; yy=hist(x,xx); bar(xx,yy/length(x)/0.8);hold on; plot(x0,y0)2、某研究者对随机抽取的一组保险丝进行了实验，测出使保险丝烧断的电流值为10.4, 10.2,12.0, 11.3, 10.7, 10.6, 10.9, 10.8, 10.2, 12.1，假设这些值α≤的条件下求出这些保险丝的溶断电流及其满足正态分布，试在置信水平0.05置信区间。
解:方法①由normfit() 函数可以直接求出置信区间，亦即溶断电流的均值为10.92，其置信区间为(10:43; 11:41)。
>> x=[10.4,10.2,12,11.3,10.7,10.6,10.9,10.8,10.2,12.1];>> [m1,s1,ma,sa]=normfit(x,0.05); m1, mam1 =10.92000000000000ma =10.43271643434768 11.40728356565233方法②采用T 检验函数即可判定是否接受均值为mean(x) 的检验，也能求出同样的均值与置信区间>> x=[10.4,10.2,12,11.3,10.7,10.6,10.9,10.8,10.2,12.1];mean(x)ans =10.92000000000000>> [H,p,ci]=ttest(x,mean(x),0.05)H =p =1ci =10.43271643434768 11.407283565652333、假设测出某随机变量的12 个样本为9.78, 9.17, 10.06, 10.14, 9.43, 10.60, 10.59, 9.98, 10.16,10.09, 9.91, 10.36，试求其方差及方差的置信区间。
数学模型数学建模第七次作业数理统计实验
数学模型第七次作业数理统计实验7.1实验目的与要求●学会对数据的参数进行估计和作相应的假设检验●学会对分布进行检验和数据的秩检验●建立相应的统计模型，并用R软件求解7.2 基本实验1. 区间估计已知某种灯泡寿命服从正态分布，在某星期所生产的该灯泡中随机抽取10只，测得其寿命（单位：小时）为1067 919 1196 785 1126 936 918 1156 920 948(1) 试问这批灯泡中大约95%的灯泡至少使用多少小时；(2) 求这批灯泡能够使用1000小时以上的概率。
解：(1)根据题意，使用R软件求解，编辑程序如下：> X<-c(1067,919,1196,785,1126,936,918,1156,920,948)> t.test(X,al="g")得到如下结果：One Sample t-testdata: Xt = 23.9693, df = 9, p-value = 9.148e-10alternative hypothesis: true mean is greater than 095 percent confidence interval:920.8443 Infsample estimates:mean of x997.1由此知道这批灯泡中大约95%的灯泡至少使用920.8443小时。
(2)> x<-c(1067,919,1196,785,1126,936,918,1156,920,948)> x[1] 1067 919 1196 785 1126 936 918 1156 920 948> pnorm(1000,mean(x),sd(x))[1] 0.5087941由此知道求这批灯泡能够使用1000小时以上的概率为50.87941%2. 假设检验I正常男子血小板计数均值为225×109/L，今测得20名男性油漆作业工人的血小板计数值(单位:109/L)220 188 162 230 145 160 238 188 247 113 126 245 164 231 256 183 190 158 224 175 问油漆工人的血小板计数与正常成年男子有无差异，并说明油漆作业对人体血小板计数是否有影响。
Matlab实验指导书(含答案)汇总
实验Matlab 操作环境熟悉、实验目的1．初步了解Matlab 操作环境。
2．学习使用图形函数计算器命令funtool 及其环境。
二、实验内容熟悉Matlab 操作环境，认识命令窗口、内存工作区窗口、历史命令窗口；学会使用format 命令调整命令窗口的数据显示格式；学会使用变量和矩阵的输入，并进行简单的计算；学会使用who 和whos 命令查看内存变量信息；学会使用图形函数计算器funtool ，并进行下列计算：1．单函数运算操作。
求下列函数的符号导数(1) y=sin(x);(2) y=(1+x)A3*(2-x); 求下列函数的符号积分(1) y=cos(x);(2) y=1/(1+x^2);(3) y=1/sqrt(1-x^2);(4) y=(x-1)/(x+1)/(x+2);求反函数(1) y=(x-1)/(2*x+3);(2) y=exp(x);(3) y=log(x+sqrt(1+x^2));代数式的化简(1) (x+1)*(x-1)*(x-2)/(x-3)/(x-4);(2) sin (x)A2+cos(x)A2;(3) x+sin(x)+2*x-3*cos(x)+4*x*sin(x);2．函数与参数的运算操作。
从y=xA2通过参数的选择去观察下列函数的图形变化(1) y1=(x+1)A2⑵ y2=(x+2)A2(3) y3=2*x^2⑷ y4=x^2+2(5) y5=x^4⑹ y6=x^2/23．两个函数之间的操作求和(1) sin(x)+cos(x)(2) 1+x+x^2+x^3+x^4+x^5乘积(1) exp(-x)*sin(x)(2) sin(x)*x商(1) sin(x)/cos(x);⑵ x/(1+x^2);(3) 1/(x-1)/(x-2);求复合函数(1) y=exp(u) u=sin(x)(2) y=sqrt(u) u=1+exp(xA2)(3) y=sin(u) u=asin(x)(4) y=sinh(u) u=-x三、设计提示1．初次接触Matlab 应该注意函数表达式的文本式描述。
《MATLAB7.x》程序设计语言（第二版）课后题答案
《MATLAB7.x》程序设计语言（第二版）部分课后习题答案（楼顺天，姚若玉，沈俊霞编著）说明：所有答案均是本人在备考过程中亲自整理的，收录了一部分题目的答案。
答案可能存在不足甚至谬误，很多算法也可能不是最优的，仅供参考。
本人尽力整理的一点点心血，希望对大家有所帮助。
第二章8、 a=rand(5,5);[i,j]=find(a>0.5);for u=1:length(i);b(u)=a(i(u),j(u));endb'10、aa=any(a');i=find(aa==0)a(i,:)=[];第四章1、function flag=isprime(m)%m为素数，flag=1%m不为素数，flag=0k=sqrt(m);flag=1;for i=2:k;if rem(m,i)==0;flag=0;breakendendm=input('input an integer\n');flag=isprime(m);if flag==1disp([num2str(m),'是素数'])elsedisp([num2str(m),'不是素数'])end4、function [x1,x2]=jfc(a,b,c)d=b^2-4*a*c;if d>0;x1=(-b-sqrt(d))/(2*a);x2=(-b+sqrt(d))/(2*a);elseif d==0;x1=-b/(2*a);a=input('a=');b=input('b=');c=input('c=');[x1,x2]=jfc(a,b,c);if x1~=x2;disp('原方程存在两个不同的根')disp(['x1=',num2str(x1)])disp(['x2=',num2str(x2)]);elseif x1==x2;disp('原方程存在两个相同的根')disp(['x1=x2=',num2str(x1)])elsedisp('原方程的根不存在')end%MATLAB上机作业，184页习题4，求二次方程的实根function [s1,s2]=solve(a,b,c)d=b^2-4*a*c;if d>0s1=(-b-sqrt(d))/(2*a)s1=(-b-sqrt(d))/(2*a)%disp(['原方程有两个不同的根'num2str(s1)'和'nums2tr(s2)]) elseif d==0s1=(-b-sqrt(d))/(2*a);s1=(-b-sqrt(d))/(2*a);%disp(['原方程有两个相同的根'num2str(s1))elsedisp('原方程无实根')endend5、x=-3:0.01:3;if x>=-3 & x<-1;y=(-x.^2-4*x-3)/2;elseif x>=-1 & x<1;y=-x.^2+1;elsey=(-x.^2+4*x-3)/2;endplot(x,y)6、%MATLAB上机作业，184页习题6，点不同按钮产生不同分布的数s=menu('请选择随机数类型','U[-10,10]','U[-5,5]','U[-1,1]','N(0,1)');switch scase 1,n=(rand(1)-0.5)*20;case 2,n=(rand(1)-0.5)*10;case 3,n=(rand(1)-0.5)*2;endn7、load a.txt[m,n]=size(a);w=zeros(m,1);for i=1:m;w(i)=(2*a(i,1)+3*a(i,2)+2*a(i,3)+4*a(i,4)+2.5*a(i,5)+a(i,6))/14.5;endw8、function k=jj(i)if i<=10k=0.1*i;elseif i<=20k=1+(i-10)*0.05;elseif i<=40k=1.5+(i-20)*0.02;elsek=1.9+(i-40)*0.01;ends=0;a=2;b=1;for i=1:15;s=s+a/b;c=b;b=a;a=a+c;ends9、a=ones(15,1);b=ones(15,1);c=ones(15,1);a(1)=2;c(1)=2;for i=2:15;a(i)=b(i-1)+a(i-1);b(i)=a(i-1);c(i)=a(i)./b(i)endsum(c)第五章3、t=1:1:10;t=t';y=[4.842,4.362,3.754,3.368,3.169,3.083,3.304,3.016,3.012,3.005]';x1=[ones(size(t)),exp(-t)];x2=[ones(size(t)),t.*exp(-t)];p1=x1\y;p2=x2\y;grid on,legend('给定数据','y1拟合','y2拟合')5、（1）function dy=cwf(t,y)dy=[5*y(1)-5*y(2)-6*y(3);3*y(1)-2*y(2)+5*y(3);2*y(1)-y(2)-4*y(3)];x0=[1,-4,5];tspan=[13,16];[t,y]=ode45('cwf',tspan,x0);plot(t,y)（2）function dy=cwf(t,y)dy=[y(1)+2*y(2)-3*y(3)+y(4);3*y(1)+y(3)-2*y(4);y(1)-2*y(2)+5*y(4);2*y(1)+3*y(2)+y(4)];x0=[1,-1,2,1];tspan=[15,16];[t,y]=ode45('cwf',tspan,x0);plot(t,y)7、t=1:10;tt=1:0.01:10;y=[15,39.5,66,85.5,89,67.5,12,-86.4,-236.9,-448.4];p1=polyfit(t,y,2);y1=polyval(p1,tt);p2=polyfit(t,y,3);y2=polyval(p2,tt);plot(t,y,'b*',tt,y1,'r-',tt,y2,'c-');legend('样本点','二次拟合','三次拟合'),grid on10&11、pa=[2 3 -4];pb=[4 -2 5];pc=[3 -2 5 6];d1=conv(pa,pb)[q2,r2]=deconv(pc,pa)[q3,r3]=deconv(pc,pb)dy1=polyder(pa,pb)[q2,r2]=polyder(pc,pa)[q3,r3]=polyder(pc,pb)12、x=-5:1:5;xx=-5:0.01:5;y=10.*exp(-abs(x));y1=interp1(x,y,xx,'nearest');y2=interp1(x,y,xx,'linear');y3=interp1(x,y,xx,'spline');y4=interp1(x,y,xx,'cubic');plot(x,y,'*',xx,y1,'r-',xx,y2,'c-',xx,y3,'k-',xx,y4,'g-')grid on,legend('样本点','最临近内插','线性内插','三次样条内插','三次曲线内插')13、a=rand(1,50);amax=max(a)vara=(std(a))^2b=randn(1,50);bmax=max(a)bmin=min(a)ub=mean(a)varb=(std(a))^214、t=[0,0.2,0.4,0.6,0.8,1,2,5]';tt=0:0.01:5;y=[1,1.51,1.88,2.13,2.29,2.4,2.6,-4]';p1=polyfit(t,y,3);y1=polyval(p1,tt);x=[ones(size(t)),exp(t)];p2=x\y;y2=p2(1)+p2(2)*exp(tt);plot(t,y,'k*',tt,y1,'r-',tt,y2,'b-')legend('样本点','三次多项式拟合','指数拟合'),grid on 17、function y=jifen1(x)y=exp(-2*x);function y=jifen2(x)y=exp(2*x);function y=jifen3(x)y=x.^2-3*x+0.5;z1=quadl('jifen1',0,2)z2=quadl('jifen2',0,2)z3=quadl('jifen3',-1,1)19、function z=jifen(x,y)z=exp(-x*y)-2*x*y;q=dblquad('jifen',0,1,-1,1)20、function dy=cwf(t,y)dy=[0.5-y(1);y(1)-4*y(2)];x0=[1,-0.5];tspan=[0,25];[t,y]=ode45('cwf',tspan,x0);plot(y(1,:),y(2,:));grid on亲爱的朋友，上文已完，为感谢你的阅读，特加送另一篇范文，如果下文你不需要，可以下载后编辑删除，谢谢！道路施工方案1、工程概况2、编制说明及编制依据3、主要施工方法及技术措施3．1施工程序3．2施工准备3．3定位放线3. 4土方开挖3．5卵石路基施工3．6天然砾基层施工3. 7高强聚酯土工格楞3．8水泥稳定砂砾基层施工3．9路缘石施工3. 10玻璃纤维土工格栅施工3．11沥青面层施工3. 12降水施工4、质量控制措施5、雨季施工安排6、安全技术措施1.工程概况本项目建设的厂址位于新疆石河子市。
实验七--离散系统分析的MATLAB实现讲解学习
实验七 离散系统分析的MATLAB 实现一、实验目的1、掌握利用MATLAB 绘制系统零极点图的方法；2、掌握离散时间系统的零极点分析方法；3、学习离散系统响应的MATLAB 求解方法；4、掌握用MATALB 实现离散系统频率特性分析的方法；5、深刻理解离散系统的系统函数零极点对系统频响的影响，可以根据 零极点知识设计简单的滤波器。
二、基本原理（一）离散系统零极点线性时不变离散系统可用线性常系数差分方程描述，即()()N Miji j a y n i b x n j ==-=-∑∑ （1）其中()y k 为系统的输出序列，()x k 为输入序列。
将式（1）两边进行Z 变换，00()()()()()Mjjj Nii i b zY z B z H z X z A z a z-=-====∑∑ （2） 将式（2）因式分解后有：11()()()Mjj Nii z q H z Cz p ==-=-∏∏ （3）其中C 为常数，(1,2,,)j q j M =L 为()H z 的M 个零点，(1,2,,)i p i N =L 为()H z 的N 个极点。
系统函数()H z 的零极点分布完全决定了系统的特性，若某系统函数的零极点已知，则系统函数便可确定下来。
因此，系统函数的零极点分布对离散系统特性的分析具有非常重要意义。
（二）离散系统零极点图及零极点分析 1、零极点图的绘制设离散系统的系统函数为()()()B z H z A z =则系统的零极点可用MATLAB 的多项式求根函数roots()来实现，调用格式为：p=roots(A) 其中A 为待求根多项式的系数构成的行矩阵，返回向量p 则是包含多项式所有根的列向量。
如多项式为231()48B z z z =++，则求该多项式根的MATLAB 命令为为：A=[1 3/4 1/8]; P=roots(A) 运行结果为： P =-0.5000 -0.2500需注意的是，在求系统函数零极点时，系统函数可能有两种形式：一种是分子、分母多项式均按z 的降幂次序排列；另一种是分子、分母多项式均按1z -的升幂次序排列。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
解：
（1）
假设：
每个年龄组的个体独立，且不受外界影响；
变量：
幼龄兔——a0(n)
中龄兔——a1(n)
老龄兔——a2(n)
按年龄分组的种群增长（Leslie矩阵）模型
得到的是：
v =
-0.7265 0.0804 - 0.4885i0.0804 + 0.4885i
-0.5484 -0.4344 + 0.3688i-0.4344 - 0.3688i
-0.4140 0.65590.6559
d =
1.3247 0 0
0 -0.6624 + 0.5623i 0
0 0 -0.6624 - 0.5623i
可知，a(n)=A*a(n-1)
A =
0 4.0000 3.0000
0.5000 0 0
0 0.2500 0
[v, d]=eig(A)
得到：
v =
-0.9474 0.9320 0.2259
-0.3158 -0.3560 -0.5914
-0.0526 0.0680 0.7741
d =
1.5000 0 0
兔子繁殖问题3
如果一对兔子每一个月可以生一对兔子，并且兔子在出生二个月以后就具有繁殖后代的能力，三个月后就离开群体。由一对兔子开始，一年可以繁殖成多少对兔子？求这个种群的稳定分布。
假设：
1、一个月生一对兔子；
2、幼兔经过两个月之后成为成兔；
3、成兔在生了兔子之后离开这个群体
变量：
一月兔——a1(n)
(1)开始每组各有1000只，求30年后各组分别有多少只;并确定种群的固有增长率和稳定分布。
(2)如果饲养者每5年出售一次动物，出售量为龄组i在这5年的增量，记出售量与该龄组存量之比为本时段收获系数H，即hi(n)xi (n)=xi (n)-xi (n-1)，H(n)=diag(h1(n), h2 (n), h3(n))。建立收获模型。
0 -1.3090 0
0 0 -0.1910
由b=A^30*a
可得：
b =
1.0e+008 *
5.3810
1.8429
0.2977
即在30年后，幼龄兔的数量是5.3810*10^8只，中龄兔的数量是1.8429*10^8只，老龄兔的数量是0.2977*10^8只。
由m=v(:,1)/sum(v(:,1))
可得：
m =
0.7200
0.2400
0.0400
幼龄兔、中龄兔、老龄兔所占的比率是：72%、24%、4%。
（2）如果饲养者每5年出售一次动物，出售量为龄组i在这5年的增量，记出售量与该龄组存量之比为本时段收获系数H，即hi(n)xi (n)=xi (n)-xi (n-1)，H(n)=diag(h1(n), h2 (n), h3(n))。建立收获模型。
因为每五年出售的仅仅是5年里的增量，所以，幼龄兔的数量保持不变
t(:,1)=v(:,1)/sum(v(:,1))
得到的是：
tHale Waihona Puke =0.43020.3247
0.2451
得出结论：
一月兔在年底占43.02%；
二月兔在年底占32.47%；
三月兔在年底占24.51%；
一群动物最高年龄为15岁(年)，繁殖周期为5年，因此每5岁一组分成3个年龄组，各组繁殖率为0, 4, 3，存活率为1/2,1/4。建立种群增长模型。
二月兔——a2(n)
三月兔——a3(n)
a1(n)=a2(n-1)+a3(n-1)
a2(n)=a1(n-1)
a3(n)=a2(n-1)
推知，a(n)=A*a(n-1)
A =
0 1 1
1 0 0
0 1 0
a=A^12*a
得到：
a =
12
9
7
结论：得到的一月兔是12对，二月兔是9对，三月兔是7对。
[v,d]=eig(A)