Inputting Power Spectral Density into System Transfer Function to obtain Output PSD
33 次查看(过去 30 天)
显示 更早的评论
I have come across the relation where and refer to the output and input spectral densities respectively, and is the system transfer function.
If I have already defined the system transfer function in the complex frequency domain using the following (for example)
s = tf('s');
H = (1+s*T2)/(1+s*(T2 + T1));
How would I go about inputting an arbitrary Power Spectral Density in this transfer function with the goal of observing the output spectral density? I am assuming the above relation is still valid in the complex frequency space but I may be wrong.
Any help would be greatly appreciated!
0 个评论
采纳的回答
Mathieu NOE
2020-12-2
hello
if your input PSD is a numerical vector defined for a frequency vector freq, you simply have to do dot products like this :
PSD_input = .... vector
freq = .... vector
s = j*2*pi*freq; % complex pulsation vector
% compute transfer function modulus
tf_mod = abs((1+s*T2)./(1+s*(T2 + T1));
% PSD out
PSD_out = (tf_mod.^2).*PSD_input; % that's all folks !
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Parametric Spectral Estimation 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!