编号28试.docVIP

数学建模校级选修课Report Of course experiment Mathematic Modeling 《数学建模》课 程 数学建模校级选修课 Report Of course experiment Mathematic Modeling 论 文 试题编号 28 班 级 电气124班 姓 名 常亚飞 学 号 3120421124 教 师 胡 钢 总 成 绩 西安理工大学理学院应用数学系 二0一四年春季学期 数学建模作业1：MATLAB 作业 数学建模作业1：MATLAB 作业1 ?作业1题目（本题共30分） 1.1 用两种方法求。 1.2绘制连续调制波形y=sin(t)sin(9t), 。 1.3解常微分方程组。 1.4 用控制语句编写一个M文件，实现某种运算，例如计算有限项级数的和等。 说明：每一道题需将相关程序与结果（图或计算结果）附在每个题目后，且有M文件的需要将M文件存入上交的作业文件夹中。 ?作业1答题正文 ?试题1.1用两种方法求 解： 方法一：用int命令进行符号积分。 编写代码如下： clear; syms x; I=int(1./(x.^3-2*x-5),x,0,2); 运行结果为： I = log(12/(643*(717788808900^(1/2)/1063390828 - 1/1286)^(1/3)) - 69*(2/(643*(717788808900^(1/2)/1063390828 - 1/1286)^(1/3)) - … 由于得出的是符号结果，故还需要在代码后加入一行double(I),将符号结果转化为数值结果，最终结果为 ans = -0.4605 方法二：以定积分的梯形积分法为原理，利用trapz(x,I)求数值解，代码如下： clear; x=0:0.1:2; I=1./(x.^3-2*x-5); trapz(x,I) 运行结果为： ans = -0.4682 如果将步长改为0.01，则结果更为准确，I=-0.4606 方法三：利用quadl(‘fun’,a,b)求出高精度的自适应递推数值积分，先编写M文件 %M 函数 fun1.m function f=fun1(x) f=1./(x.^3-2*x-5); MATLAB代码为： clear; quad(fun1,0,2) vpa(quad(fun1,0,2),10) %以十位有效数字显示 最后得到结果为： ans = -0.4605017397 比较三种方法，如果要求结果精度不高的话，选择第一种方法结果更为简单准确，第二种方法结果的准确程度取决于所取步长，而第三种方法需要编写M文件及调用，较繁琐。 ?试题1.2绘制连续调制波形y=sin(t)sin(9t), 。 解：依题意编写代码如下： t=(0:pi/100:pi);%产生一组步长为pi/100的列向量 y1=sin(t)*[1,-1];%生成一组对称的正弦函数构成最终图像的包络 y2=sin(t).*sin(9*t); t3=pi*(0:9)/9;%生成一个序列 y3=sin(t3).*sin(9*t3); plot(t,y1,r:,t,y2,b,t3,y3,bo) %同一个平面上画出三个函数的图像 xlabel(X);ylabel(Y); legend(包络,包络,正弦函数,零点); axis([0,pi,-1,1])%设定横坐标，纵坐标的范围 运行得到图形： ?试题1.3解常微分方程组。 解：用dsolve法求解析解代码如下： clear; s=dsolve(Du=v,Dv=w,Dw=u,u(0)=0,v(0)=0,w(0)=1); u=simplify(s.u) v=simplify(s.u) w=simplify(s.u) 求得结果为：u = exp(t)/3-cos((3^(1/2)*t)/2)/(3*exp(t)^(1/2)) - (3^(1/2)*sin((3^(1/2)*t)/2))/(3*exp(t)^(1/2)) v=exp(t)/3-cos((3^(1/2)*t)/2)/(3*exp(t)^(1/2)) - (3^(1/2)*sin((3^(1/2)*t)/2))/(3*exp(t)^(1/2)) w=exp(t)/3-cos((3^(1/2)*t)/2)/(3*exp(t)^(1/2)) - (3^(1/2)*sin((3^(1/2)*t)/2