（模拟数据从prob_data401.mat获得）目的l 强调：要非常慎用数值导数计算 。
l 练习mat数据文件中数据的获取 。
l 实验数据求导的后果l 把两条曲线绘制在同一图上的一种方法 。
解答（1）从数据文件获得数据的指令假如prob_data401.mat文件在当前目录或搜索路径上clearload prob_data401.mat （2）用diff求导的指令dt=t(2)-t(1);
yc=diff(y)/dt;
%注意yc的长度将比y短1plot( 。

32、t,y,b,t(2:end),yc,r)grid on （3）用gradent求导的指令（图形与上相似）dt=t(2)-t(1);
yc=gradient(y)/dt;
plot(t,y,b,t,yc,r)grid on 说明l 不到万不得已 ， 不要进行数值求导 。
l 假若一定要计算数值导数 ， 自变量增量dt 要取得比原有数据相对误差高1、2个量级以上 。
l 求导会使数据中原有的噪声放大 。
n 2.采用数值计算方法 ， 画出在区间曲线 ， 并计算 。
提示l 指定区间内的积分函数可用cumtrapz指令给出 。
l 在计算要求不太高的地方可用find指令算得 。
目的l 指定区间内的积分函数的数值计算法和cumtrapz指令 。
l find指令的应用 。
解答dt=1e-4;
t=0:dt:10;
t=t+(t=0)*eps;
f=sin(t)./t;
s=cumtrapz(f)*dt;
plot(t,s,LineWidth,3)ii=find(t=4.5);
s45=s(ii) s45 =1.6541n 3.求函数的数值积分 ， 并请采用符号计算尝试复算 。
提示l 数值积分均可尝试 。
l 符号积分的局限性 。
目的l 符号积分的局限性 。
解答dx=pi/2000;
x=0:dx:pi;
s=trapz(exp(sin(x 。

来源：(未知)

【学习资料】网址：/a/2021/0121/0021123973.html

标题：MATLAB|MATLAB教程2012a习题解答1-7章完整版-张志涌-北航( 四 )