信噪分离120151123.docVIP

freq_zone=freq1-freq2; 频率子集区间，freq1是频率子集上限,freq2是频率子集下限 for n=0:r+1 dw1(n+1)=freq1+freq_zone*n; 在频率子集区间等间隔取r+2个频率点 for i1=0:r+1 ai1_temp=1; for k=0:r+1 ai_temp=ai1_temp/(cos(dw1(k+1))- cos(dw1(i1+1)); end ai(i1)=(-1)^i1*ai_temp; end for i1=0:r+1 delta_temp1=0; delta_temp2=0; for n=0:r+1 delta_temp1= delta_temp1+ai(n)*Hd(n); delta_temp2= delta_temp2+(-1)^n*ai(n)/W(n); end delta1(i1)= delta_temp1/delta_temp2; 用remez设计FIR低通滤波器问题。matlab 那位高手会设计这个滤波器，滤波器阶数P=50.通带[0,0.4pi].阻带[0.5pi,pi].通带最大偏差delta1=0.01.阻带最大偏差0.001.在[0,0.4pi]U[0.5pi,pi]间隔中等间隔采样N=200个点。谢谢！ 我来回答 窗体底端 窗体顶端 窗体底端 回答 (1) 网友回答?2014-04-22 b=remez(50,[0,0.4,0.5,1],[1,1,0,0],[0.01,0.001]);freqz(b,1) 采用频率fs?： 8000Hz 通带截至频率：1500Hz(归一化后为0.375) 阻带截至频率：2000Hz(归一化后为0.5) 通带波纹：0.01(即0.1737dB) 阻带衰减：0.1(即20dB) [n, fo, mo, w]? = remezord([1500 2000], [1 0], [0.01 0.1], 8000); b = remez(n, fo, mo, w); (1)[H,w]=freqz(B,A,N)(2)[H,w]=freqz(B,A,N,’whole’) 用matlab语言计算等纹波最佳逼近法设计线性相位FIR滤波器: 采频率fs?通带截至频率阻带截至频率通带波纹阻带衰减 [n, fo, mo, w]? = remezord([w1 w2], [1 0], [delta1 delta2], fs); b = remez(n, fo, mo, w); [H,w1]=freqz(B,1,N);求出滤波器频率响应  (1)中B和A分别为离散系统的系统函数分子、分母多项式的系数向量，返回量H则包含了离散系统频响在 0~pi范围内N个频率等分点的值（其中N为正整数），w则包含了范围内N个频率等分点。调用默认的N时，其值是512。 采频率fs?： 8000Hz 通带截至频率：1500Hz(归一化后为0.375) 阻带截至频率：2000Hz(归一化后为0.5) 通带波纹：0.01(即0.1737dB) 阻带衰减：0.1(即20dB) [n, fo, mo, w]? = remezord([1500 2000], [1 0], [0.01 0.1], 8000); b = remez(n, fo, mo, w); %[h w] = freqz(b, 1, 128); %plot(w/pi, abs(h)); freqz(b, 1, 128); 用Matlab语言计算相关函数的频谱： n=0:N-1; k =0:N-1; %序列范围 WN=exp(-j*2*pi/N); %设定wn因子 nk=n*k; %赋值nk转置为n行1列矩阵与k相乘 WNnk=WN.^nk; %矩阵点幂 Rk=r4m*WNnk; %相关函数的频谱 Xk=xn*WNnk; %输入信号x(n)的频谱 (5)信号频谱还原为时域信号： n=0，1,2……N-1 用Matlab语言将信号频谱还原为时域信号： n=0:N-1; k =0:N-1; %序列范围 WN1=exp(j*2*pi/N); %设定wn因子 kn=k*n; %赋值nk转置为n行1列矩阵与k相乘 WNnk1=WN1.^kn; %矩阵点幂 x2n=Xmtk*WNnk1; %求出时域信号 （6）小波变换 h=0.7071; [C,L]=wavedec(s,4,’db9’);四层小波变换 ca4=appcoef(C,L,’db9’,4); 第四层近似系数 cd1=detcoef(C,L,’db9’,1); 第一层细节系数 cd2=detcoef(C,L