上传者: 38733333
|
上传时间: 2022-03-29 13:35:27
|
文件大小: 49KB
|
文件类型: -
% 示例 1 f1 = @(t) t.^2; Lapf1 = @(s) LapTrans(f1,s); sc = 1:0.1:10; sp = 1:10; yc = 2./sc.^3; yp = arrayfun(Lapf1, sp); 图1) 情节(sc,yc,'r',sp,yp,'*'); Legend('精确变换', '近似变换'); title('数值拉普拉斯变换 (t^2)')
% 示例 2
f1 = @(t) sin(t); Lapf1 = @(s) LapTrans(f1,s); sc = 1:0.1:10; sp = 1:10; yc = 1./(1 + sc.^2); yp = arrayfun(Lapf1, sp); 图(2) 情节(sc,yc,'r',sp,yp,'*'); Legend('精确变换', '近似变换'); title('数字拉普拉斯变换(sin(