?? q603.m
字號:
%《MATLAB在電子信息課程中的應用》第六章例6.3程序q603
% 求任意線性系統的沖擊函數
% 電子工業出版社出版 陳懷琛 吳大正 高西全合著 2001年10月
a=input('多項式分母系數向量a=(書上取 poly([0,-1+2j,-1-2j,-2,-5]) ) ');
b=input('多項式分子系數向量b=(書上取 [8,3,1] ) ');
[r,p] = residue(b, a), % 求留數
disp('解析式h(t)=Σr(i)*exp(p(i)*t)')
disp('給出時間數組t=[0:dt:tf]')
dt=input('dt= '); % 給定時間數組
tf=input('tf= ');
t=0:dt:tf;h=zeros(1,length(t)); % h的初始化
for i=1:length(a)-1 % 根數為a的長度減1
h = h+ r(i)*exp(p(i)*t); % 疊加各根分量
end
plot(t,h),grid
?? 快捷鍵說明
復制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號
Ctrl + =
減小字號
Ctrl + -