matlab fft freqz,【急】请教 fft、freqz、bode 求相频响应的区别及原因

该楼层疑似违规已被系统折叠隐藏此楼查看此楼

已知冲击响应，要画传递函数的相频响应，先后分别使用 fft、freqz、bode函数，发现fft和freqz的趋势是一样的，但是角度数值不同。而bode的趋势则和fft、freqz都不同。代码如下，SP1是函数的冲击响应，采样率4800。

多谢各位啦！ ^___^

%% 用 fft 函数

H_complex_two = fft(SP1,4800);

H_complex = H_complex_two(1:end/2);

F = 1:2400;

H_angle = unwrap(angle(H_complex))/pi*180; % Degree

figure;plot(F,H_angle)

axis([20,150,-900,-100]);

grid on

ylabel('Phase (Degree)');title('FFT');

xlabel('Frequency (Hz)');

%% 用 freqz

[H,W] = freqz(SP1,1,4800);

W = W/pi*2400;

Ha = unwrap(angle(H))/pi*180;

figure;plot(W,Ha)

axis([20,150,-1600,-800]);grid on

xlabel('Frequency (Hz)');ylabel('Phase (Degree)');title('freqz')

%% 用 bode

htf = tf(SP1.',1,1/4800);

Wtf = 2*pi*(1:2400); % 幅频响应和相频响应分析频率(rad/sec)

[Htf_Mag,Htf_Phase_temp] = bode(htf,Wtf); % 原次级通道频响

Htf_Phase = zeros(1,2400);

for j = 1:2400

Htf_Phase(j) = Htf_Phase_temp(1,1,j);

end

Htf_Phase = unwrap(Htf_Phase);

figure;plot(Wtf/2/pi,Htf_Phase,'k');

axis([20,150,-30,40]);grid on;

xlabel('Frquency (Hz)');ylabel('Phase (Degree)');title('bode');

matlab fft freqz,【急】请教 fft、freqz、bode 求相频响应的区别及原因相关推荐

【 MATLAB 】信号处理工具箱之 fft 案例分析
上篇博文:[ MATLAB ]信号处理工具箱之fft简介及案例分析介绍了MATLAB信号处理工具箱中的信号变换 fft 并分析了一个案例,就是被噪声污染了的信号的频谱分析. 这篇博文继续分析几个小案例 ...
【 MATLAB 】信号处理工具箱之fft简介及案例分析
目录 Syntax Description Y = fft(X) Y = fft(X,n) Y = fft(X,n,dim) Examples Noisy Signal Syntax Y = fft( ...
MATLAB中的快速傅里叶变换FFT与IFFT
背景 FFT (Fast Fourier Transform)是离散傅立叶变换的快速算法,可以将一个信号从时域变换到频域.同时与之对应的是IFFT(Inverse Fast Fourier Trans ...
【STM32F407的DSP教程】第28章 FFT和IFFT的Matlab实现（幅频响应和相频响应）
完整版教程下载地址:http://www.armbbs.cn/forum.php?mod=viewthread&tid=94547 第28章 FFT和IFFT的Matlab实现(幅 ...
[Matlab科学计算] 频谱分析和FFT算法总结
频谱分析是一种非常重要的信号处理方法,在机械设备故障诊断.振动系统分析.电力系统.无线电通信.信息图像处理和自动控制等学科中都有重要应用.频谱分析的核心是1965年Cooely-Tukey发表的快速傅 ...
Matlab快速傅里叶变换程序（FFT）编写
Matlab快速傅里叶变换程序(FFT)编写 (不利用Matlab内置fft) 最近在学习数字信号处理,正好到快速傅里叶变换,写一个基于基2FFT算法的傅里叶变换.利用原理:DIT-FFT(时域抽取法 ...
深入浅出解释FFT（四）——fft分析信号频率和相位
很基础的问题往往很重要,做仿真时候有一个点的差错都会导致结果的错误.在网上找了前人写的东东,总结下希望对大家有帮助,让大家少走一些弯路. 1. 信号的时域采样点N和频域采样点数相同 %####### ...
matlab伯德图稳定性,11系统的稳定性Bode稳定判据介绍.ppt
7.6 由伯德图判断系统的稳定性;幅相曲线(-1,j0)点左侧的负实轴对数幅频特性L(ω)>0(即零分贝线以上的区域)对数相频特性-180°线 ;可见,同样可以利用伯德图来判别系统的稳定性.这 ...
matlab抛物偏微分方程,抛物型方程差分求解跪求MATLAB解抛物型偏微分方程的程序...
为什么抛物线方程与圆方程联立不能使用韦达定理很容易了解到,抛物线和圆的交点均在X轴上方\"其实这时你应该注意到一点就是,这两个交点的纵坐标是相等的,所以其实对应的是一个y值,也就是你列的一 ...
matlab 求复数的模长,matlab计算带有复数的函数，最后求复数函数的模，结果里面却有...
该楼层疑似违规已被系统折叠隐藏此楼查看此楼 matlab计算带有复数的函数,最后求复数函数的模,结果里面却有复数标志i,输入自变量,得到的结果确实是实数,不知道怎么回事; 代码如下: syms n2 ...

matlab fft freqz,【急】请教 fft、freqz、bode 求相频响应的区别及原因

matlab fft freqz,【急】请教 fft、freqz、bode 求相频响应的区别及原因相关推荐

最新文章

热门文章