发布时间:2022-08-17 13:26
序号 | 内容 |
---|---|
1 | 建立函数文件 |
2 | 实现简单三角波的积分 |
3 | 实现简单三角波的微分 |
键入以下代码,然后保存,保存的文件名字一定要和函数名字一样
function yt = fun1(t)%函数名为fun1
yt = tripuls(t,4,0.5);
end
h=0.001;t=-3:h:3;
y=fun1(t);
y1=diff(y)/h;
subplot(2,1,1);
plot(t,y);
title('f(t)');
subplot(2,1,2);
plot(t(1:length(t)-1),y1);
title('df(t)/dt')
t=-3:0.1:3;
y=tripuls(t,4,0.5);
for x=1:length(t)
y2(x)=quad('fun1',-3,t(x));
end
subplot(2,1,1);
plot(t,y);
title('f(t)');
subplot(2,1,2);
plot(t,y2);
title('integral of f(t)');