数学建模培训作业(MATLAB编程部分)
数学建模Matlab练习4
MATLAB练习四1. 绘制两个一元函数的曲线图（例如y=sin(x)、y=cos(x)），加上标注说明，包括坐标、抬头、曲线说明等。
%lianxi41.mx=0:pi/30:2*pi;y=sin(x);plot(x,y,'-ob','markersize',10);hold ony1=cos(x);plot(x,y1,'-.dm','markersize',10);legend('sinx','cosx');title('sin和cos函数的曲线图')xlabel('x');ylabel('y');2. 寻找一组数据拟合多项式，写出M文件，包括拟合的次数，数值的计算，数学表达式的呈现，图形的呈现等。
%lianxi42x=0:pi/20:pi;y=sin(x);n=input('n=');p=polyfit(x,y,n);x1=0:pi/20:pi*2;y1=sin(x1);y2=polyval(p,x1)plot(x1,y1,'-ob',x1,y2,':dm');legend('原曲线','拟合曲线');xlabel('x');ylabel('y');title('sin的拟合曲线和原曲线');axis([-0.2,6.5,-1.2,1.2])3. 求解一个代数方程，写出M文件，应该包括函数式M文件和命令式M文件。
%fc3.mfunction f=fc3(x)f='3*x-1=0';%lianxi43clc;clear;close;sprintf('原方程和其解析解的结果为：')s=solve(fc3)sprintf('保留10位有效数字的结果为：')digits(10)sprintf('数值解的结果为：')x=vpa(s)4. 求解一个代数方程组，写出M文件，应该包括函数式M文件和命令式M文件。
数学建模matlab例题参考及练习
数学实验与数学建模实验报告学院：专业班级：姓名：学号：完成时间：年月日承 诺 书本人承诺所呈交的数学实验与数学建模作业都是本人通过学习自行进行编程独立完成，所有结果都通过上机验证，无转载或抄袭他人，也未经他人转载或抄袭。
若承诺不实，本人愿意承担一切责任。
承诺人：年 月 日数学实验学习体会（每个人必须要写字数1200字以上，占总成绩的20%）练习1 一元函数的图形1． 画出x y arcsin =的图象.2． 画出x y sec =在],0[π之间的图象.3． 在同一坐标系中画出x y =，2x y =，3x y =，3x y =，x y =的图象. 4． 画出3232)1()1()(x x x f ++-=的图象，并根据图象特点指出函数)(x f 的奇偶性. 5． 画出)2ln(1++=x y 及其反函数的图象. 6． 画出321+=x y 及其反函数的图象.练习2 函数极限1．计算下列函数的极限.(1)xxx4cos12sin1lim4-+π→.程序：sym x;f=(1+sin(2*x))/(1-cos(4*x));limit(f,x,pi/4)运行结果：lx21ans =1(2).程序：sym x;f=(1+cos(x))^(3*sec(x));limit(f,x,pi/2)运行结果：lx22ans =exp(3)(3)22)2(sinlnlimxxx-ππ→.程序：sym x;f=log(sin(x))/(pi-2*x)^2;limit(f,x,pi/2)运行结果：lx23ans =-1/8(4)212lim xxex→.程序：x xx sec32)cos1(lim+π→sym x ;f=x^2*exp(1/x);limit(f,x,0)limit(f,x,0,'right')limit(f,x,0,'left')运行结果：lx24ans =NaNans =Infans =%左极限为零，存在，右极限为无穷大，在x 趋近于零时函数没有极限(5))215(lim 122x x x x +-∞→.程序：sym x ;f=5*x^2/(1-x^2)+2^(1/x);limit(f,x,inf)运行结果：>> lx25ans =-4(6)x x x x x -+-→32112lim .程序：sym x ;f=(x^2-2*x+1)/(x^3-x);limit(f,x,1)运行结果：>> lx26ans =0 (7)x x x 11lim 20-+→.程序：sym x ;f=(sqrt(1+x^2)-1)/x;limit(f,x,0)运行结果：>> lx27ans =0 (8))3sin(cos 21lim 3π--π→x xx .程序：sym x ;f=(1-2*cos(x))/sin(x-pi/3);limit(f,x,pi/3)运行结果：>> lx28ans =3^(1/2) (9)tgxx x )1(lim 0+→.程序：sym x ;f=(1/x)^tan(x);limit(f,x,0,'right')运行结果：>> lx29ans =(10)xx arctgx )2(lim π+∞→.程序：sym x ;f=(2/pi*atan(x))^x;limit(f,x,inf,'left')运行结果：>> lx210ans =Inf2． 解方程012=-⋅x x .程序：sym x ;X=solve(x*2^x-1)运行结果：>> lx202X =lambertw(0, log(2))/log(2)%方程有两个解3． 解方程1sin 3+=x x .程序：sym x ;X=solve(3*sin(x)+1-x)运行结果：>> lx203X =-0.538470451711254993610615326557454． 解方程03=++q px x .(p 、q 为实数)程序：X=solve('x^3+p*x+q=0','x')运行结果：X =((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3) - p/(3*((p^3/27 + q^2/4)^(1/2)- q/2)^(1/3))p/(6*((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)) - ((p^3/27 + q^2/4)^(1/2)- q/2)^(1/3)/2 - (3^(1/2)*i*(p/(3*((p^3/27 + q^2/4)^(1/2) -q/2)^(1/3)) + ((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)))/2p/(6*((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)) - ((p^3/27 + q^2/4)^(1/2)- q/2)^(1/3)/2 + (3^(1/2)*i*(p/(3*((p^3/27 + q^2/4)^(1/2) -q/2)^(1/3)) + ((p^3/27 + q^2/4)^(1/2) - q/2)^(1/3)))/2练习 3 导数及偏导数计算1.求下列函数的导数. (1))11)(1(-+=x x y程序：sym x ;f=(sqrt(x)+1)*(1/sqrt(x)-1);diff(f)运行结果：>> lx31ans =(1/x^(1/2) - 1)/(2*x^(1/2)) - (x^(1/2) + 1)/(2*x^(3/2))(2)x x x y ln sin =程序：sym x ;f=x*sin(x)*log(x);diff(f)运行结果：>> lx32ans =sin(x) + log(x)*sin(x) + x*cos(x)*log(x)2.求下列参数方程所确定的函数的导数.(1)⎩⎨⎧==t y t x 44程序：sym t ;f1=t^4;f2=4*t;diff(f2)/diff(f1)运行结果：>> lx321ans =1/t^3(2)⎩⎨⎧-=+=arctgt t y t x )1ln(2程序：sym t ;f1=log(1+t^2);f2=t-atan(t);diff(f2)/diff(f1)运行结果：>> lx322ans =-((t^2 + 1)*(1/(t^2 + 1) - 1))/(2*t)3.求下列隐函数的导数. (1)22ln y x x y arctg +=程序：syms x y ;f=atan(y/x)-log(sqrt(x^2+y^2));yx=-diff(f,x)/diff(f,y)运行结果；>> lx331yx =(x/(x^2 + y^2) + y/(x^2*(y^2/x^2 + 1)))/(1/(x*(y^2/x^2 + 1)) - y/(x^2 + y^2))(2)xy y x =程序：syms x y ;f=x^y-y^xyx=-diff(f,x)/diff(f,y)运行结果：>> lx332f =x^y - y^xyx =(x^(y - 1)*y - y^x*log(y))/(x*y^(x - 1) - x^y*log(x))4.设x e y x cos =，求)4(y .程序：sym x ;f=exp(x)*sin(x);diff(f,x,4)运行结果：>> lx34ans =(-4)*exp(x)*sin(x)5.验证x e y x sin =满足关系式:22=+'-''y y y程序： sym x ;f=exp(x)*sin(x);y2=diff(f,x,2);y1=diff(f,x,1);y=f;y2-y1*2+2*y=='0'运行结果：>> lx35ans =1%运行结果为1表示y2-y1*2+2*y=='0'成立6.设)ln(y x x u +=，求22x u ∂∂，22y u∂∂，y x u ∂∂∂2. 程序：syms x y ;f=x*log(x+y);uxx=diff(f,x,2)uyy=diff(f,y,2)f1=diff(f,x);uxy=diff(f1,y)运行结果：>> lx36uxx =2/(x + y) - x/(x + y)^2uyy =-x/(x + y)^2uxy =1/(x + y) - x/(x + y)^27.求下列多元隐函数的偏导数y zx z ∂∂∂∂,.(1)1cos cos cos 222=++z y x程序：syms x y z ;f=(cos(x))^2+(cos(y))^2+(cos(z))^2-1;zx=-diff(f,x)/diff(f,z)zy=-diff(f,y)/diff(f,z)运行结果：>> lx371zx =-(cos(x)*sin(x))/(cos(z)*sin(z))zy =-(cos(y)*sin(y))/(cos(z)*sin(z))(2)xyz e z =程序：syms x y z ;f=exp(z)-x*y*zzx=-diff(f,x)/diff(f,z)zy=-diff(f,y)/diff(f,z)运行结果：>> lx372f =exp(z) - x*y*zzx =(y*z)/(exp(z) - x*y)zy =(x*z)/(exp(z) - x*y)8.证明函数22)()(ln b y a x u -+-=(b a ,为常数)满足拉普拉斯方程: 02222=∂∂+∂∂y u x u（提示：对结果用simplify 化简）练习4 积分计算1.计算下列不定积分. (1)⎰+dx x x 12 (2)⎰+x xdx 2sin 12sin2.计算下列定积分.(1)⎰e xdx x 1ln (2)⎰ππ342sin dx x x 3.求⎰+t dx x x x 12)ln (ln 1并用diff 对结果求导.4.求摆线)cos 1(),sin (t a y t t a x -=-=的一拱(π≤≤20t )与x 轴所围成的图形的面积.5.计算二重积分(1)⎰⎰≤++122)(y x dxdy y x (2)⎰⎰≤++x y x dxdy y x 22)(22 6.计算⎰+L ds y x 22 L 为圆周)0(22>=+a ax y x7.计算⎰++-L dy y x dx y x )()(2222，其中L 为抛物线2x y =上从点(0，0)到点(2，4)的一段弧.练习5 matlab 自定义函数与导数应用1.建立函数x x a a x f 3sin 31sin ),(+=，当a 为何值时，该函数在3π=x 处取得极值，它是极大值还是极小值，并求此极值.2.确定下列函数的单调区间.（1）7186223---=x x x y （2）)0(82>+=x x x y3.求下列函数的最大值、最小值.（1）2332xx y -=41≤≤-x （2）312824≤≤-+-=x x x y练习6 matab 矩阵运算与数组运算1． 计算（1）⎪⎪⎪⎭⎫ ⎝⎛--521111204321+⎪⎪⎪⎭⎫ ⎝⎛--232002101041221 （2）⎪⎪⎭⎫ ⎝⎛-01301213⎪⎪⎪⎪⎪⎭⎫ ⎝⎛030101020501⎪⎪⎪⎭⎫ ⎝⎛-205101 (3)52422⎪⎪⎭⎫ ⎝⎛- 2．设⎪⎪⎪⎭⎫ ⎝⎛-=243121013A ,⎪⎪⎪⎭⎫ ⎝⎛-=112111201B ,求满足关系B X A =-23的X ．练习7 矩阵与线性方程组１．求下列矩阵的秩．(1)⎪⎪⎪⎭⎫ ⎝⎛-321110021 (2)⎪⎪⎪⎪⎪⎭⎫ ⎝⎛4820322513454947513253947543173125 ２．求下列矩阵的行列式，如可逆，试用不同的方法求其逆矩阵．（1）⎪⎪⎪⎭⎫ ⎝⎛--285421122 （2）⎪⎪⎪⎪⎪⎭⎫ ⎝⎛---6201111121324321 ３．设X ⎪⎪⎪⎭⎫ ⎝⎛-111012111=⎪⎪⎪⎭⎫ ⎝⎛521234311求X ． ４．解下列线性方程组．(1)⎪⎪⎩⎪⎪⎨⎧=--+=+-+=+-+=+-+6223312433862344224221432143214321x x x x x x x x x x x x x x x x (2)⎪⎪⎩⎪⎪⎨⎧-=+--=+--=-+-212201432143214321x x x x x x x x x x x x练习8 常微分方程与级数 求1-6题微分方程的通解1．1222+='y y y x 2．x y x y dx dy -+= 3．x x x y y +='cos 4．1)2sin cos (='+y y y x 5．x e y y y x 2cos 3=-'+'' 6．x x y y sin 14++=+'' 求7、8题初值问题的解7．⎪⎩⎪⎨⎧==-++-+=10)2(212222x y dx dy x xy y y xy x8．⎪⎪⎩⎪⎪⎨⎧===++==0000222,02V dt dxx x x a dt dx n dt x d t t9．给出函数x x e x f x x cos 2sin )(+=在点0=x 的7阶taylor 展开式以及在x=1处的 5阶taylor 展开式．10．判别下列级数的敛散性，若收敛求其和．(1)++++7151311(2)∑∞=+112n nntgπ11．求幂级数∑∞=--22)1(nnnnnx的和函数．12．求函数项级数∑∞=-1)2sin)1(nnnn xπ的和函数．。
matlab文件数学建模作业.docx
1.1(1)代码:syms x;m=433;limit((m*x-sin(m*x))/(x A3),^0) 结果：ans=81182737/6(2 )代码:syms x;m=433;limit((m*x-sin(m*x))/(x A3)/x/Inf) 结果：ans=01.2代码：syms x;m=433;y=exp(x)*cos(m*x/1000);diff(y,x,2)结果:(812511*cos((433*x)/1000)*exp(x))/1000000 - (433*sin((433*x)/1000)*exp(x))/5001.3代码：dblquad(,exp(x.*x+y.*y)'/0/l/0/l)结果：ans = 2.13941.4代码：syms x;m=433;i nt(x A4/(m A 2+4*x*x),x)结東：(81182737*atan((2*x)/433))/32 - (187489*x)/16 + x A3/121.5代码：syms x;m=433;diff(exp(x)*cos(m*x)/10)结果：-C231619285105554448813874688*cos(433*x)*exp(x) - 5350116883610368850333216*sin(433*x)*exp(x) 1.6代码：syms x;m=433;simple(taylor(sqrt(m/1000.0+x)/x/0/,Order,/5))结果：4330A(1/2)*(- (390625000*x A4)/35152125121 + (625000*x A3)/81182737 ・(1250*x A2)/187489 + (5*x)/433 + 1/100)1.7代码：n二2;A二g];while n<20;A(n+1)=A(n)+A(n-l);n=n+l；end;fprintf(,A=(,);for i = l:size(A/2);fprintf( %u,/A(i));if(i~=size(A/2)) fprintf(7);end;end;fprintfC)1);结果：0.24250.97011.8f.mfunction y 二f(x)if(x>=0 && x<=0.5);y=2*x; elseif(x>0.5 && x< = l);y=2*(l-x); else;y=NaN;e nd;1.10 代码：t=0:0.01*pi:6*pi;plot3(cos(t)z sinm/jhold on;grid on;plot3(2*cos(t)/2*sin(t)/t/k ,);1.11 代码：m=433;A=[4r 2/2;-3A5;l /5*m /3]； B=[l /3A-2/0/-3;2r l /l];det(A)/2*A-B /A*B /A.*B /A*B A (-l)/ inv(A)*B /A A 2/A ,⑴x=-3:0.1:3;m=433;a = m/600;f=@(x,a,b) l/sqrt(2*pi)/a*exp(-(x-b).A2./(2*a*a)); plot(x ，f(xaO»;hold on; plot&,f&a-i));hold on; plot (x,f (xai));⑵x 二二 433;f=@(x /a /b) l/sqrt(2*pi)/a*exp(-(x-b).A2./(2*a*a)); plot(x,f(x,l,0));hold on; plot(x,f&20»;hold on; plot(x,f(xAO));hold on;p I ot(x,f (x, m/100r 0));1.13m=433;ezmesh('433*x.*x+y A 4,,30); coIorm 叩([OQlj);4 9 108361.14绘图m=433;x=-2:0.01:2;y=x.A 5-m/200*x-0.1;plot(x,y);grid on;求根vpa(solve(y),3)得根为:1.22,-0.0462,-1.2求单调区间：vpa(solve(diff(y))z3)得:0.811,-0.811根据图象:单调增区间为(・8 , -0.811)4 0.811, +8 )单调减区间为:(-0.811A811)推理：当x<-0.811是"的导数大于0 ,故递增；-0.811<x<0.811时，y的导数小于0 ,故递减；x>0.811时，y的导数大于0 ,故递增。
数学建模Matlab作业
x=[-3:.1:3]; %确定x的值
y1=2.*x+5;
y2=x.^2-3.*x+1;
subplot(1,2,1); %建立二个窗口,并在第一个编辑
plot(x,y1);
axis([-4 4 -5 20]); %划定坐标范围
subplot(1,2,2); %建立二个窗口,并在第一个编辑
plot(x,y2);
axis([-4 4 -5 20]); %可省略,因为坐标范围就是以第二幅图为准的程序效果图:
x=[190 33 45 42 45];
explode = [0 1 0 0 0]; %输入分离切面向量
pie(x,explode,{'生活费','资料费','电话费','购买衣服费','其它费用'}); %绘图程序效果图:
第三题代码:(这一题我是按照点乘来做的,而且精度调的比较大,如果是按矩阵做,只需要把z的表达式一改就行了)
[x,y]=meshgrid(-5:.01:5); %x的精度为0.01
z=x.^2+y.^2+sin(x.*y);
meshc(x,y,z) %绘制等高图
执行效果图:
第四题代码:
x=0:pi/5:4*pi;
y=sin(x).*exp(-x/10);
xx=0:.1:4*pi; %取间隔为0.1
yy=spline(x,y,xx) %三次样条插值
plot(x,y,'o',xx,yy) %绘图,并用’o’标记执行效果图:。
北航5系—MATLAB编程部分练习大作业
MATLAB编程部分练习大作业班级_xxxx_ 姓名:_xxx_ 学号：_xxxxx_一、判断题：1）由矩阵a、b组成的方程组x*b=a的解，可以利用x=a/b进行求解。
（对）2）@符号作为子函数数柄，记录了子函数的详细信息，可作为主函数的输出参数。
（对）3）num2str函数表示将字符串或字符数组转换为数字或矩阵。
（错）4）n=nargin(‘fun’)用于返回函数fun定义的输入参数的个数。
（对）5）某m函数文件定义行语句为：Function varargout=myfun(x,varargin)，如果调用方式为myfun(a,b,c)则varargin是长度为2的单元数组。
（对）6）Roots用于求多项式等于0的根，根用行向量表示。
（错）7）最小二乘法拟合函数ployfit(x,y,n)中输入参数n为样本数据个数。
（错）8）RTW（Real-Time Workshop）是MatLab提供的实时仿真工具，只要硬件性能允许，所有SimuLink建立的模型都可以通过RTW实现实时运行。
（错）9）在MATLAB中函数tf2zp()可以用来求传递函数的零极点和增益。
（对）10）MATLAB语言是由C语言开发的，执行效率高，一般能够快速实时运行。
（错）二、选择题：1）下列（）函数可以创建一个m行、m列单位矩阵（对角线上为数值1，其它均为0）：A)eye(m); B) zeros(m); C)ones(m); D)rand(m);2）下列变量名中（）是合法的：A) Char_1; B) x*y; C) x\y; D) end2）定义a=['abc';'cde'] ，运行reshape(a,[3,2])语句产生矩阵为（）A) ['ab','cd','de']; B) ['ab';'cd';'de']; C) ['ad';'cc';'be']; D) ['ab';'cc';'de'] 3）定义A=[1,2,3;4,5,6;7,8,9]; 运行A(1,:)*A(:,3)语句计算结果为（）：A) [ 30,36,42]; B) 36; C) 30; D) 424）下列计算表达式中（）是错误的A）y=eval(‘sin(t)’); B) y=feval(‘sin(t)’,t); C) y=feval(@sin,t); D) y=feval(‘sin’,t);5）如果n阶方阵A可逆，则求A逆矩阵的命令是（）:A) inv(A); B) det(A); C) rank(A); D)以上都不对6）多项式x4-12x3+25x+116的系数向量为（）：A) [1,-12,25,116]; B) [116,25,-12,1]; C) [1,-12,0,25,116]; D) [116,25,0,-12,1];7）Varagin作为函数声明语句中最后一个输入参数，用于传送不定数目的参数，其数据类型为：A) 数值数组；B)单元数组；C)结构数组；D)函数数柄8）令y=[-1,6,15,-7,31,2,-4,-5]，运行[ynew,indx]=sort(y)，请问所得indx(2)的具体数值为：A) 1；B) 4；C) 7；D) 89）下列函数中，（）用于绘制x,y,z三维空间曲线：A) plot；B) plot3；C) mesh；D) peak三、编程计算题：1）利用结构数据类型定义变量A，用于储存某个学生以下三方面信息：姓名——Dafei、排名——第3以及三门课的成绩：88，99和98。
数学建模matlab作业
数统应数20121323003 王妍数统应数20121323022 胡可旺1 生成5阶矩阵，使其元素满足均值为1，方差为4的正态分布；代码：y=1+sqrt(4)*randn(5)结果:2，生成一个20行5列矩阵A，其元素分别以概率0.7、0.1、0.2取自0、1、4这三个整数，然后将“删除A的全零行”之后的结果赋给矩阵B；for i = 1:20for j = 1:5p = rand();if p<=0.7A(i,j) = 0;elseif p>0.7&&p<=0.8A(i,j) = 1;elseA(i,j) = 4;endendendB = [];for i = 1:20flag = true;for j = 1:5if A(i,j)~=0flag=false;break;endendif ~flagB = [B;A(i,:)];endendB3, （函数文件）对于输入的正整数，判断其是否为素数function [ output_args ] = u3( n )%U3 Summary of this function goes here% Detailed explanation goes herey=1;for i=2:1:sqrt(n)if rem(n,i)==0y=0;breakendendif(y==1)fprintf('是素数')elsefprintf('不是素数')end4，（函数文件）对于输入的正整数，求出其中0的个数；function [ i ] = u4( n )%U4 Summary of this function goes here% Detailed explanation goes herei=0;while n/10~=0if rem(n,10)==0i=i+1;endn=fix(n/10);end5，（函数文件）有n个人围成一圈，按序列编号。
MATLAB编程作业
《Matlab编程训练》作业专业学生姓名班级学号指导教师完成日期实训一 MATLAB 语言介绍和数值计算1．先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存变量。
122sin 851z e =+.2. 已知 1234413134787,2033657327A B --⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦，求下列表达式的值：(1) A+6*B 和A-B+I （其中I 为单位矩阵） A+6*B:A-B+I:(2)A*B和A.*BA*B程序：A=[12 34 -4;34 7 87;3 65 7]B=[1 3 -1;2 0 3;3 -2 7]c=A*B结果：A.*B程序：A=[12 34 -4;34 7 87;3 65 7]B=[1 3 -1;2 0 3;3 -2 7]D=A.*B结果：(3)A^3和A.^3A^3程序：A=[12 34 -4;34 7 87;3 65 7]E=A^3结果：A．^3程序：A=[12 34 -4;34 7 87;3 65 7]C=A.^3(4)A/B及B\AA/B程序：A=[12 34 -4;34 7 87;3 65 7]B=[1 3 -1;2 0 3;3 -2 7]C=A/B结果：B\A程序：A=[12 34 -4;34 7 87;3 65 7]B=[1 3 -1;2 0 3;3 -2 7]D=B\A结果：(5)将矩阵C=B\A的右下角2*2子矩阵赋给D, 并（3）保存变量（mat文件）程序：A=[12 34 -4;34 7 87;3 65 7];B=[1 3 -1;2 0 3;3 -2 7];C=B*inv(A);D=C(2:3,2:3)结果：3. 求得矩阵⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=34157864653434533145A 的每行最大元素所在的位置？(至少两种方法) 第一种：A=[5 14 33;45 43 3;65 4 6;78 15 34][RowMax Order]=max(A')结果：第二种：A=[5 14 33;45 43 3;65 4 6;78 15 34][Max_num,index]=max(A,[],2)结果：实训二 MATLAB 编程基础1. 求[25,1258]之间能被15整除的数的个数。
数学建模MATLAB培训2
模再型编设写在一t=个0时脚刻本缉M私文艇发件现走 Y
私c船lea，r 走私船以速度a平行于Y轴
正ts向=0行:0驶.05，:0缉.5私; 艇按速度b以指
向x走0=私[0船0]的'; 方向行驶(b>a)，在任
Q(c, at)
意[时t x刻]=otd缉e45私(@艇ji位si,于ts,点x0P);(x,y),则 M文件以ex4t1作为文件
0.3500 12.8170 4.5552 15.0000 7.0000
0.4000 13.9806 6.1773 15.0000 8.0000
0.4500 14.7451 8.0273 15.0000 9.0000
0.5000 15.0046 9.9979 15.0000 10.0000
3、数据拟合、polyfit、人口预报.
v=[x(:,2);y(2:end,2)];
%速度函数
L=length(h);dt=(T+12)/(L-1);
%用三点公式求加速度
ac(1)=(-3*v(1)+4*v(2)-v(3))/2/dt;
%起始点加速度
for k=2:L-1
ac(k)=(v(k+1)-v(k-1))/2/dt;
%中间点加速度
设置的搜索路径，然后在解方程的脚本文件中调用：
[t, x]=ode23(@f,ts,x0,options) 或
ts=[a, b]表示自变量的 初值到终值的区间
[t, x]=ode45('f',tspan,x0,options) 或ts=[t0,t1,t2,…tf]
options用于设置误差限（缺省时相对误差为10-3，绝对误差 为10-6），用语句实现为：
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
数学建模培训作业 （MATLAB 编程部分）
1． 请使用switch 语句将百分制的学生成绩转换为五级制的成绩输出。
2． 猜数游戏：首先由计算机随机产生一个 [1,100] 之间的一个整数，然后由用户猜测所产生的这个数。
根据用户猜测的情况给出不同的提示，如果猜测的数大于产生的数，则显示 “High” ，小于则显示 “ Low ” ，等于则显示 “You won ！”，同时退出游戏。
用户最多有 7 次机会。
3． 使用for 结构计算1+2+3+…+100。
4． 设计一个九九乘法表。
5． 使用while 结构计算1+2+3+…+100。
6． 求1!+2!+ …+10!的值。
7． 编程生成三对角矩阵。
1
10000011100000001110000000111000000011100
000001110000000111000000011100
1
1
轾犏犏犏犏犏犏犏
犏犏
犏犏犏犏犏犏犏犏臌
8． 计算分段函数的值，要求根据不同的x 输入，给出相应的结果。
223135
x x y x x ìï-ïïï+ïï=íï-ïïïï+ïî 110011x x x x ?-< ?>
9． 已知1
1111
1(1)435721
n n p -?+-++-- ，编程求 的近似值。
10．
输入下面的矩阵
12345678
910111213141516A 轾犏犏犏=犏犏
犏臌
编程求该矩阵的对角线元素之和，并找出最大和最小元素的值以及其所在的行、列号。
11．
求水仙花数。
如果一个三位数的个位数、十位数和百位数的立方和等于该数自身，
则称该数为水仙花数。
编一程序求出所有的水仙花数。
12． 给定两个实数a 、b 和一个正整数n ，计算()k a b +和()k a b -，其中n k ,,2,1 。
13． 编写函数，生成一个1!，2!，…，n!的阶乘表。
14． 编一个函数统计字符串中单词的个数。
15．
求n 阶勒让德多项式的值，其递归公式为：
(,)((2*
1)**(1,)
(p n x n x p n x n p n
n
=----- (0,)
1;(1,p x p x x
== 16． 编写一个判断任意输入正整数是否为素数的函数文件，并在命令窗口调用。
17．
编写一个万年历计算程序，当输入年月日后，能够计算出该日是星期几。