出错 Q2 (第 8 行) [b, a] = cheby2(n, Rs, Ws); % 根据阶数和规格化截止频率求解系数
时间: 2024-04-01 21:36:59 浏览: 16
这个错误提示是 MATLAB 的错误提示,意思是在使用 cheby2 函数时出错了。可能的原因是输入参数的值不正确,比如阶数 n、规格化截止频率 Ws、衰减 Rs 等参数的值设置不当。您可以检查一下输入参数的值是否符合要求,特别是检查一下阶数 n 是否设置过高,导致求解系数的过程出现了数值问题。如果输入参数的值都正确,但仍然出现错误,那可能是 cheby2 函数本身出现了问题,您可以尝试使用其他滤波函数替代。
相关问题
出错 Untitled (line 23) [b, a] = cheby1(n, Rp, [wp1, wp2], 'stop');
这个错误可能是由于切比雪夫I型数字带阻滤波器的阶数n计算错误导致的。请检查计算阶数的代码是否正确,以下是计算阶数的代码:
```
% 计算通带最大衰减Rp和阻带最小衰减As对应的通带和阻带截止频率
ep = sqrt(10^(Rp/10) - 1);
a = sqrt(10^(As/10) - 1);
B = ws2 - ws1;
n = ceil(acosh(a/ep) / acosh(B / (2 * pi))); % 计算阶数
```
请检查代码中的Rp、As、ws1、ws2、wp1、wp2和Fs的值是否正确,并检查计算阶数的公式是否正确。如果还有问题,请提供更多的错误信息,我会尽力帮助你解决。
请根据要求帮我修改下面的代码,要求如下:用 MATLAB 直接法设计椭圆型数字高通滤波器,要求:通带ωp=0.3π,Rp=1dB;阻带ωs=0.2π,As=20dB。请描绘滤波器的绝对和相对幅频特性、相频特性、零极点分布图,列出系统传递函数式。代码如下:ws1=0.15; ws2=0.85; %数字滤波器的阻带截止频率 ws=[ws1,ws2]; wp1=0.25;wp2=0.75; %数字滤波器的通带截止频率 wp=[wp1,wp2]; %求数字系统的频率特性[H,w]=freqz(b,a); Rp=1;As=20; %输人波器的通阻带衰减指标 [n,we]=cheb1ord(wp,ws,Rp,As); %计算阶数n和截止频率 [b,a]=cheby1(n,Rp,we); %直接求数字带通滤波器系数 [H,w]=freqz(b,a); %求数字系统的频率特性 dbH=20 * log10((abs(H)+eps)/max(abs(H)));%化为分贝值 subplot(2,2,1); plot(w/pi,abs(H)); title('幅频响应'); subplot(2,2,2); plot(w/pi,angle(H)); title('相频响应'); subplot(2,2,3); plot(w/pi,dbH); title('幅频响应 dB'); subplot(2,2,4); zplane(b,a); %根据H(z)绘制零极点
根据要求,修改后的 MATLAB 代码如下:
%数字高通滤波器设计
wp = 0.3*pi; %数字滤波器的通带截止频率
ws = 0.2*pi; %数字滤波器的阻带截止频率
Rp = 1; %通带最大衰减量
As = 20; %阻带最小衰减量
%使用cheby1函数进行设计
[n, wn] = cheb1ord(wp, ws, Rp, As); %计算阶数n和截止频率
[b, a] = cheby1(n, Rp, wn, 'high'); %直接求数字高通滤波器系数
%绘制幅频响应、相频响应、幅频响应dB和零极点分布图
[H, w] = freqz(b, a);
dbH = 20*log10(abs(H));
subplot(2,2,1); plot(w/pi, abs(H)); title('数字高通滤波器幅频响应');
subplot(2,2,2); plot(w/pi, angle(H)); title('数字高通滤波器相频响应');
subplot(2,2,3); plot(w/pi, dbH); title('数字高通滤波器幅频响应dB');
subplot(2,2,4); zplane(b, a); title('数字高通滤波器零极点分布图');
%打印系统传递函数式
syms z;
H = poly2sym(b,z) / poly2sym(a,z);
disp('数字高通滤波器系统传递函数式:');
disp(H);
相关推荐
![rar](https://img-home.csdnimg.cn/images/20210720083606.png)
![zip](https://img-home.csdnimg.cn/images/20210720083736.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)
![](https://csdnimg.cn/download_wenku/file_type_ask_c1.png)