How I can plot the magnitude and phase response oh the function

Question

Helda 2013-10-19

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/90780-how-i-can-plot-the-magnitude-and-phase-response-oh-the-function

回答： Abood Abusafea 2022-2-2

i want do that for this function y=(4*sin(50*t)/(6*t)

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

sixwwwwww 2013-10-19

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/90780-how-i-can-plot-the-magnitude-and-phase-response-oh-the-function#answer_100247

在 MATLAB Online 中打开

Dear Helda, here is an example showing amplitude and phase plots of your defined function:

t = 1:100;
y = 4 * sin(50 * t) ./ (6 * t);
figure, plot(t, abs(y)), title('Amplitude plot')
figure, plot(t, angle(y)), title('Phase plot')

I hope it helps. Good luck!

8 个评论
显示 6更早的评论隐藏 6更早的评论

Helda 2013-11-28

plz help me :/ for filter with a transfer function of H=(0.1667*s^3-0.5*s^2+0.5*s+0.1667)/(s^3-0.3333*s) what is the filter order ? what the filter type ? (FIR,IIR) why? plot the frequency response of the filter ? find the filter gain by matlab ?

Engr Muhammad Amir Shahzad 2020-4-11

@sixwwww. Why you used *fftshift* in plotting magnitude of signal?

请先登录，再进行评论。

Answer 2

David 2013-10-20

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/90780-how-i-can-plot-the-magnitude-and-phase-response-oh-the-function#answer_100317

在 MATLAB Online 中打开

"How I can plot the magnitude and phase response of the function

y=(4*sin(50*t)/(6*t)"

From what I've read, it seems you want the amplitude and phase of this function in the frequency domain. If this is the correct assumption to make, then you will need to make a lot more specifications. You will need to know your sampling rate, Fs, and either your time of observation or the number of points you have sampled. After you have figured these out, look into 'fft' function MATLAB provides.

An alternate route would be to use MATLAB's symbolic toolbox. You will want to look into how to create symbolic variables and symbolic equations as well as how to use the 'laplace', 'subs', and 'ezplot' functions. Good luck.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 3

Afshin Aghayan 2017-7-24

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/90780-how-i-can-plot-the-magnitude-and-phase-response-oh-the-function#answer_275291

编辑：Afshin Aghayan 2017-8-2

look at the following Matlab function, it can calculate phase spectrum as well as amplitude spectrum with a perfect accuracy:

https://www.mathworks.com/matlabcentral/fileexchange/63965-amplitude-and-phase-spectra-of-a-signal--fourier-transform-

This program calculates amplitude and phase spectra of an input signal with acceptable accuracy especially in the calculation of phase spectrum.The code does three main jobs for calculation amplitude and phase spectra. First of all, it extends the input signal to infinity; because for calculation Fourier transform(FT) (fft function in Matlab), we consider our signal is periodic with an infinite wavelength, the code creates a super_signal by putting original signal next to itself until the length of super_signal is around 1000000 samples, why did I choose 1000000 samples? Actually, it is just based on try and error!! For most signals that I have tried, a supper signal with 1000000 samples has the best output.

Second, for calculating fft in Matlab you can choose different resolutions, the Mathwork document and help use NFFT=2^nextpow2(length(signal)), it definitely isn't enough for one that wants high accuracy output. Here, I choose the resolution of NFFT=100000 that works for most signals.

Third, the code filters result of FT by thresholding, it is very important step! For calculating phase spectrum, its result is very noisy because of floating rounding off error, it causes during calculation "arctan" even small rounding off error produces significant noise in the result of phase spectrum, for suppressing this kind of noise you can define a threshold value. It means if amplitude of specific frequency is less than predefined threshold value (you must define it) it put zero instead of it.

These three steps help to improve the result of amplitude and phase spectra significantly.

IF YOU USE THIS PROGRAM IN YOUR RESEARCH, PLEASE CITE THE FOLLOWING PAPER:

Afshin Aghayan, Priyank Jaiswal, and Hamid Reza Siahkoohi (2016). "Seismic denoising using the redundant lifting scheme." GEOPHYSICS, 81(3), V249-V260. https://doi.org/10.1190/geo2015-0601.1