| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 1044 人关注过本帖
标题:(有问题) 广义切比雪夫的传输零点提取程序
只看楼主 加入收藏
caoaihua3
Rank: 1
等 级:新手上路
帖 子:38
专家分:0
注 册:2006-4-20
收藏
 问题点数:0 回复次数:0 
(有问题) 广义切比雪夫的传输零点提取程序

我最近在做一个滤波器的东西,要用matlab 计算.下面的东西就是在从书上找出的一段程序,但是好像不是行.


请大家帮帮忙.

         广义切比雪夫的传输零点提取:


clc
clear
RL=22;
c=sqrt(10^(RL/1O)-1);
ep=1/c;
ws=5
Ar=22;
Ass=70;
m=2000;
b=0;
for As=Ass:0.1:100;
    b=0;
    v1=acosh(10^((Ar+As)/20);
    for n=4:1:14
        y=v1-(n-1)*acosh(ws);
        if y<0
            disp('不合要求');
           break;
       end;
       v2=cosh(v1-(n-1)*acosh(ws));
       w1=(6-6.4);
       dw=x/m;
       if w1<ws
           k=0;
           for w=(4.4):dw:(6);
               Ls=10*log10((1+ep^2*cosh(acosh(w*w1-1)/(w-w1))+(n-1)*acos(w))*cosh(acosh(w*w1-1)*acosh(w))));
               if Ls<(As)
                   break;
               end;
               k=k+1;
           end
           if k==(m+1)
               disp('得到解');
               b=1;
               break;
           end;
       end
       if b==1
           break;
       else
           n=n+1;
       end
   end
   if b==1
       for w=-3;0.01;3;
           As1=10*log10(abs(1/(1+ep^2*cosh(acosh((w*w1-1)/(w-w1)+(n+1)*acosh(w))*cosh(acosh((w*w1-1)/(w-w1))+(n-1)*acosh(w))))));
           figure(1);
           hold on
          plot(w,As1,'-');
      end
      break;
  end
end
legend('传输特性');
xlabel('归一化频率');
ylabel('衰减(db)');
As
n
k
wl
      
      

运行一后提示
??? Error using ==> filter
Not enough input arguments.

请大家帮帮忙,改一下,谢谢了.

[此贴子已经被作者于2006-4-26 11:52:51编辑过]

搜索更多相关主题的帖子: 切比雪夫 传输 
2006-04-26 08:48
快速回复:(有问题) 广义切比雪夫的传输零点提取程序
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.034698 second(s), 8 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved