Matlab 讲义chap4 习题.docVIP

chap4 习题 根据题给的模拟实际测量数据的一组和 试用数值差分diff或数值梯度gradient指令计算，然后把和曲线绘制在同一张图上，观察数值求导的后果。（模拟数据从prob_data401.mat获得） 〖答案〗 首先将prob_data401.mat文件复制到当前工作目录，或者将搜索路径扩展到该文件所在目录，然后导入该文件 clear load(D:\教学\本科教学\2008-2009上学期\科学计算实践\prob_data401.mat) plot(t,y) axis([0,7,-2,2]) hold on yd=diff(y); yg=gradient(y); d=pi/1000; %注意：此处需要按照t的取值设置 dyd=yd/d; dyg=yg/d; plot(t(1:end-1),dyd,r) %plot(t,dyg,r) hold off 采用数值计算方法，画出在区间曲线，并计算。 符号计算方法 syms t x ft=sin(t)/t yx=int(ft,t,0,x) ft = sin(t)/t yx = sinint(x) sinint(4.5) ans = 1.6541 ezplot(yx) 数值计算方法 dt=0.05; t=0:dt:10; t(1)=eps; %这里用来修正t=0时的除零错误 ft=sin(t)./t; plot(t,ft) yx=dt*cumtrapz(t,ft); plot(t,yx) 〖答案〗 s45 = 1.6541 求函数的数值积分，并请采用符号计算尝试复算。 dx=pi/100; x=0:dx:pi; fx=exp(sin(x).^3); s=dx*trapz(fx) s = 5.1370 syms x fx=exp(sin(x)^3) s=vpa(int(fx,x,0,pi)) s =5.1370135567176141893469089750876 用quad求取的数值积分，并保证积分的绝对精度为。 fx=exp(abs(x)).*abs(sin(x)); quad(fx,-5*pi,1.7*pi,1e-9) ans = 3.6177e+006 求函数在区间中的最小值点。 ft=@(t)(sin(5*t)^2*exp(0.06*(t^2))-1.5*t*cos(2*t)+1.8*abs(t+0.5)) [xn0,fval,exitflag,output]=fminbnd(ft,-5,5) ft = @(t)(sin(5*t)^2*exp(0.06*(t^2))-1.5*t*cos(2*t)+1.8*abs(t+0.5)) xn0 = -1.2850 fval = -0.1860 exitflag = 1 output = iterations: 12 funcCount: 13 algorithm: golden section search, parabolic interpolation message: [1x112 char] 〖答案〗 最小值点是 -1.28498111480531 相应目标值是 -0.18604801006545 设，用数值法和符号法求。 数值计算方法 首先定义函数 function ydot=DyDt_ex4_6(t,y) %y1=y(t); %y2=dy(t)/dt; %Dy1=dy(t)/dt=y2 %Dy2=d2y(t)/t2 ydot=[y(2);1+3*y(2)-2*y(1)]; 然后进行如下步骤 tspan=[0,1]; y0=[1;0]; [tt,yy]=ode45(@DyDt_ex4_6,tspan,y0); plot(tt,yy(:,1)); 符号计算方法 dsolve(D2y-3*Dy+2*y=1,Dy(0)=0,t) ans = 1/2+exp(2*t)*C1-2*exp(t)*C1 〖答案〗 数值解 y_05 = 0.78958020790127 符号解 ys = 1/2-1/2*exp(2*t)+exp(t) ys_05 = .78958035647060552916850705213780 已知矩阵A=magic(8)，（1）求该矩阵的“值空间基阵”B ；（2）写出“A的任何列可用基向量线性表出”的验证程序（提示：利用rref检验）。 〖答案〗 三组不同的基 B1 = 64 2 3 9 55 54 17 47 46 40 26 27