matlab使用方法讲述.ppt

[x,y]=dsolve(Dx=x+y,Dy=-x+y,x(0)=1,y(0)=2,t) x = exp(t)*(2*sin(t)+cos(t)) y = exp(t)*(2*cos(t)-sin(t)) 1、常微分方程的符号求解法 建模应用： 1.人口问题 已知美国1790-1900年的人口数据 年 1790 1800 1810 1820 1830 1840 1850 1860 人口(百万) 3.9 5.9 7.2 9.6 12.9 17.1 23.2 31.4 年 1870 1880 1890 1900 人口(百万) 38.6 50.2 62.9 76.0 估计指数模型 中的参数 N0 与 r. N0 与 r不是线性关系 对指数模型两边取对数，得 记y = lnN(t), a=lnN0 ,则有 a 与 r 是线性关系 程序如下： x=0:11; %离散的时间点 y=[3.9,5.9,7.2,9.6,12.9,17.1,23.2,31.4,38.6,50.2,62.9,76]; %对应时间点的数据 z=log(y); a=polyfit(x,z,1); %线性多项式插值 思考：如果你拟合的数据一直取到2000年，得到还是这个拟合方程么？ 2. 划艇比赛 通过课上的比例分析我们推得： 题中所给的表格 赛艇 2000米成绩 t (分) 种类 1 2 3 4 平均 单人 7.16 7.25 7.28 7.17 7.21 双人 6.87 6.92 6.95 6.77 6.88 四人 6.33 6.42 6.48 6.13 6.32 八人 5.87 5.92 5.82 5.73 5.84 n t 1 7.21 2 6.88 4 6.32 8 5.84 最小二乘法 利用4次国际大赛冠军的平均成绩对模型 t n – 1/ 9 进行检验 t n 1 2 4 8 7.21 6.88 6.32 5.84 ? ? ? ? 与模型巧合！ 程序： y=[7.21 6.88 6.32 5.84]; x=[1 2 4 8]; z=log(x); t=log(y); a=poly(z,t,1); 例1.23 画出函数 在-5 ? x ? 5的图形。 解: x=-5:.1:5;↙ %取绘图横坐标向量点x y=sin(x.^2); ↙ plot(x,y),grid on↙ %图1.9曲线y=sinx^2 图1.9曲线 1.5.3 plot(x1,y1,x2,y2,x3,y3…) 功能:在同一图形窗口画出多条不同颜色曲线，曲线 关系为 。 例1.25 在同一图形窗口画出三个函数 , ,y = x的图形,-2? x ? 2 。 解:Matlab命令： x=-2:.1:2; ↙ plot(x,cos(2*x),x,x.^2,x,x) ↙ legend(cos(2x),x^2,x) ↙ 图1.11 例1.25的绘图结果 1.5.4 ezplot(F,[xmin,xmax]) 功能:画出符号函数F在区间[xmin,xmax]内的图像 说明: F是只含有一个变量的函数。如果区间 [xmin,xmax]缺省,默认区间为 . 解: syms t　↙ Ezplot(‘2/3*exp(-t/2)*cos(3/2*t)’,[0,4*pi])↙ 例1.26 绘制 在 间的图形. 图1.12 符号函数的图形 1.5.5 二维特殊图形 例1.27练习指令bar,stairs,pie,pie3,stem,area. 解: x=1:5; subplot(2,3,1),bar(x),title(直方图)↙ subplot(2,3,2),stairs(x),title(阶梯图)↙ subplot(2,3,3),stem(x,rp),title(火柴棍图)↙ subplot(2,3,4),pie(x),title(饼图)↙ subplot(2,3,5),pie3(x),title(三