MATLAB计算方法与实现.docVIP

（1）：恢复窗口：在Desktop中下拉式菜单中的Desktop Layout,选择Default来恢复。 （2）：在同一坐标系中，画出函数y=x^3-x-1和y=abs(x)*sin5x的图像。 x=-1:0.1:2;y1=x.^3-x-1; y2=abs(x).*sin(5*x); plot(x,y1,k,x,y2,:ro) legend(y1=x.^3-x-1,y2=abs(x).*sin(5*x)), xlabel(x),ylabel(y),title(y1,y2画在同一坐标系中) （3）：根据数据建立一个人口增长模型。（百万） 年 1850 1860 1870 1880 1890 1990 1910 1920 人口 23.2 31.4 38. 50.2 62.9 75.995 91.972 105.711 年 1930 1940 1950 1960 1970 1980 1990 2000 人口 123.203 131.699 150.697 179.329 203.212 226.505 249.633 281.422 解题思路：首先将表格中的数据转变为MATLAB能处理的矩阵，然后把人口数量看成年份的函数并绘制出这一函数图形。根据数学相关理论，用3,4阶多项式拟合这一函数，拟合时不计2000年的数据对，而是将这对数据用来检验并确定模型。最后用确定的模型预测2010年美国人口。 在Command window 中输入： t=1850:10:1990; p=[23.2,31.4,38.6,50.2,62.9,75.995,91.972,105.711,123.203,131.699,150.697,179.323,203.212,226.505,249.633]; %读取数据 plot(t,p,’o’);axis([1850 2020 0 400]); title(‘Population of the U.s.1850-1990’); ylabel(‘Millions’);%绘制出数据的函数图形并加以修饰 f1=polyfit(t,p,3);f2=polyfit(t,p,4);%对数据做3,4阶多项式拟合，结果分别为f1和f2 v=[polyval(f1,2000),polyval(f2,2000)];%计算当t=2000时多项式f1,f2的值 abs(v-251.422) %计算两个模型与2000年人口数的绝对误差。 ans = 30.8937 由计算结果可以确定3阶多项式可作为此问题的数学模型，因此进一步进行预测： V1=polyval(f1,2010) V1 = 311.6136 即由3阶多项式拟合可预测到2010年美国人口将达到311.6148百万人。 （4）：研究捕食者与被捕食者（Lotka-Voltrra）模型的相互作用系数?和?的影响。 X’1=x1- x1x2 ? X’2=x2+x1x2? X1(0)=x2(0)=1,0=t=10 建立M文件，再调用微分方程求解器ode45求出数值解。 function dx=lotka(t,x) global alpha beita dx=[x(1)-alpha*x(1)*x(2);-x(2)+beita*x(1)*x(2)]; 第二步：调用ode45命令求数值解。在Command Window中输入下面的命令。 global alpha beita alpha=0.1;beita=0.2; [t,x]=ode45(lotka,[0,10],[1,1]);%求模型的数值解 plot(t,x(:,1),k,t,x(:,2),:k)%画模型的图像解 title(alpha=0.1;beita=0.2,Lotka-Voltrra), xlabel(t),ylabel(x1,x2), legend(x1’,x2), 第二章 MATLAB的变量是以字母开头，由字母，数字，下划线组成。同一字母的大小写表示不同的变量。 常见的变量有四种：数值变量（double）字符变量（char）符号变量（sym）结构变量（struct）。MATLAB存储变量时用图标和文字来注明变量的类型。例如输入： x1=1/3 x2=’1/3’ x3=sym(‘1/3’) x4.x=1/3 可以在workspace中看到它们各自不同的图标和变量类型的文字解释。它们分别是double型,char型,sym型,struct型变量。 下面主要讨论数值型变量的有关问题。 在MATLAB中，不用事先定义变量的特性（实数或复数），也不用事先定义变量的维数。一旦某个变量被赋值，则该变量连同其值将被存放于Workspace中，直到该变量被清除