编号28试题.docVIP

《数学建模》课 程 论 文 试题编号 28 班 级 电气124班 姓 名 常亚飞 学 号 3120421124 教 师 胡 钢 总 成 绩 西安理工大学理学院应用数学系 二0一四年春季学期 (作业1题目（本题共30分） 1.1 用两种方法求。y=sin(t)sin(9t), 。 1.4 用控制语句编写一个M文件实现某种运算，例如计算有限项级数的和等。图或计算结果）附在每个题目后，且有M文件的需要M文件存入上交的作业文件夹中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)) - … 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))/(3*exp(t)^(1/2)) (试题1.4用控制语句编写一个M文件，实现某种运算，例如计算有限项级数的和等。clear; s(1)=1;N=100; for i=2:N s(i)=s(i-1)+1/i; end s(N) plot(1:100,s,*) 运行结果为：ans =5.1874 由图像可知：其部分和序列是单调递增形式趋于无穷大。 (作业1教师评分