matlab初步讲课.ppt

实验1、用克莱默法则求方程的根 A=[4 1 0;1 -1 5;2 2 3]; b=[6;14;-3]; x=zeros(1,3); if det(A)~=0 for i=1:3 C=A; C(:,i)=b; x(i)=det(C)/det(A); end x else disp(‘no answer’) end 实验2 昆虫繁殖问题 一、问题 一种昆虫，最长寿命为六周，将其分为三组：第一组0~2周龄；第二组2~4周龄；第三组4~6周龄。第一组为幼虫（不产卵），第二组每个成虫在两周内平均产卵100个，第三组每个成虫在两周内平均产卵150个。假设每个卵的成活率为0.09,第一组和第二组的昆虫能顺利进入下一个成虫组的存活率分别为0.1和0.2。设现有三个组的昆虫各100只，计算第2周、第4周、第6周后各个周龄的昆虫数目，并考虑下面问题： （1）以两周为一时间段，分析这种昆虫各周龄组数目演变趋势。昆虫数目是无限增长还是趋于灭亡？ （2）如果使用一种除虫剂可以控制昆虫的数目，使得各组昆虫的成活率减半，问这种除虫剂是否有效？ 了解建立数学模型的基本方法，运用线性代 数知识解决实际问题。 二、实验目的 1． 建立三个周龄组昆虫的种群动态分析的数学模型 二、实验内容与要求 记 为第k个时间段昆虫数分布向量。 当k= 0，1，2，3时，X(k)分别表示现在第0周后,2周后,4周后, 6周后的昆虫数分布向量。 如果有矩阵L 存在，使得 X(k+1) = LX(k) ，则称 L 为莱斯利矩 阵。 符号与假设： 以两周为一个时间段，某一时刻各周龄组的昆 虫数量用一个三维向量表示。由题设，初始时刻0~2周龄、 2~4周龄、4~6周龄的昆虫数量分别为： 建立数学模型并写出X(k)和X(k+1)的递推关系式，以及莱斯利矩阵L. 昆虫状态图 成长 产卵 L=[0 9 13.5;0.1 0 0;0 0.2 0]; X(:,1)=[100 100 100]; for i=1:3 X(:,i+1)=L^i*X(:,1); end X %求出各阶段的昆虫数目 2．根据数学模型计算并将数据结果填写到下表中空 白单元处. k (现在) k=1 (2周后) k=2(4周后) k=3(6周后) x1 100 ?2250 ? 360 ? 2052 x2 100 ? 10 ?225 ? 36 x3 100 ? 20 ?2 ?45 2．在除虫剂的作用下，情况将有所变化。重建三个 周龄组昆虫种群动态分析的数学模型 建立数学模型并写出X(k)和X(k+1)的递推关系式，以及莱斯利矩阵L. 除虫剂作用下昆虫状态图 成长 产卵 建立数学模型 format long L=[0,9,13.5;0.05,0,0;0,0.1,0]; X(:,1)=[100,100,100]; for i=1:30 X(:,i+1)=L^i*X(:,1); end X 所用命令 5．分析除虫剂是否有效 k 现在 k=1 (2周后) k=2 (4周后) k=3 (6周后) k=30 (60周后) x1 100 ? 2250 180 ?1019.25 ?0.15925050050414 x2 100 ?5 ?112.5 ?9 ?0.01084037648776 x3 100 ?10 ?0.5 ?11.25 ?0.00146906663282 考虑到环境保护，如果使用的除虫剂使昆虫存活率为零，将产生负面影响。你认为有效的除虫剂应该至少使昆虫存活率减少到多少为好？此时莱斯利矩阵的正特征值为多大？ 三、问题思考 例 画函数Z=(X+Y).^2的图形. x=-3:0.1:3; y=1:0.1:5; [X,Y]=meshgrid(x,y); Z=(X+Y).^2; surf(X,Y,Z) shading flat %将当前图形变得平滑 surf(x,y,z) ：画出数据点（x，y，z）表示的曲面 数据矩阵。分别表示数据点的横坐标、纵坐标、函数值 meshz(X,Y,Z)：在网格周围画一个curtain图(如,参考平面) [X,Y]=meshgrid(-3:.125:3); Z=peaks(X,Y);%产生一个凹凸曲面 meshz(X,Y,Z)%如果使用mesh则如右图所示 例 绘peaks