MATLAB的答案

如何做一个传递函数最小相位?

39视图(30天)
Govind Narayan Sahu
Govind Narayan Sahu 2023年5月16日在20分
评论道: Govind Narayan Sahu2023年5月21日16:41
亲爱的MATLAB社区,
我有一个工厂非最小相位传递函数。我想让它稳定的最小相位系统,这样我可以没有instaability逆。
%定义传递函数
H =特遣部队([-4.8 16000 0 0],[4.8 16080 286800 51160000]);
isminphase ([-4.8 16000 0 0], [4.8 16080 286800 51160000])
谢谢!
3评论
显示2年长的评论
Govind Narayan Sahu
Govind Narayan Sahu 2023年5月21日16:41
嗨 @Sam翟 ,
我想设计一个传递函数(Hinv)行为反向H, H的最终大小* Hinv(大约)等于1。这是因为当我通过一个输入信号,它在大约10 Hz放大输出。
H是可逆的,如果它有一个最低的阶段。
我读了这篇文章: https://doi.org/10.1080/10910344.2012.698970
我能创建一个最小相位的频率响应,但不能得到正确的传递函数使用“invfreqz”和“freqz”。
%定义传递函数
num = (-4.8 16000 0 0);
穴= (4.8 - 1.608 5.116 2.868 e04 e05 e07);
H =特遣部队(num穴);%逆传递函数窝/ Num
w = 2π* * (0.1:0.01:100);
%计算复杂的倒频谱
c_hat_n =传输线(日志(abs(挤压(freqresp (H, w)))));
N =长度(c_hat_n);
望远镜=圆(N / 2) + 1;
c_hat_n(望远镜:N) = 0;
c_hat_n = c_hat_n * 2;
h_hat_min = c_hat_n;
机构= exp (fft (h_hat_min));
w11 =角(机构);
[b] = invfreqz(机构)、w11 3, 3, [], 1000);
[h, w1] = freqz (b, a, N);
次要情节(211)
持有在;情节(w / 2π,abs(挤压(freqresp (H, w))));
情节(w / 2π,abs(机构),“r”)
情节(w / 2π,abs (h)、‘g’)
包含(的频率(赫兹))
ylabel (“大小(N / N)、“)
次要情节(212)
持有在;情节(w / 2π,阶段(挤压(freqresp (H, w))) + 2 *π);
情节(w / 2π,阶段(机构),“r”)
情节(w / 2 /π2 * pi-w1,‘g’)
包含(的频率(赫兹))
ylabel (的阶段,(程度))

登录置评。

在回答这个问题。

答案(0)

类别

找到更多的在转换在帮助中心和文件交换

标签

s manbetx 845


释放

R2022b

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!