求设计一个带通滤波器到matlab程序

%

用切比雪夫最佳一致逼近设计线性相位FIR带通滤波器；

%信号为0.5hz,

0.9hz,

1.1hz和1.5hz的正统信号叠加组如谈成

%通带为[0.9,1.1]

%频谱分辨率与信号实际长度N成正比

clear

all

f1=0.5f2=0.9f3=1.1f4=1.5t=0:1203N=length(t)fs=10M=512

x1=sin(2*pi*(f1/fs)*t)+sin(2*pi*(f2/fs)*t)+sin(2*pi*(f3/fs)*t)+sin(2*pi*(f4/fs)*t)

figure(1)

subplot(211)plot(t,x1)title('原信号')

y=fft(x1)

f=(0:1/N:1/2-1/N)*fs

subplot(212)plot(f,abs(y(1:N/2)))gridxlabel('hz')%处理前频谱

wc1=2*f2/fswc2=2*f3/fswc3=2*f4/fs%归一化角频率，用于下面的f1

f1=[0

wc1-0.05

wc1

wc2

wc2+0.05

1]

A=[0

0

1

1

0

0]%设置带裤橡余通或带阻，1为带通，0为带阻

weigh=[1

1

1

]%设置通带和阻带的权重

b=remez(60,f1,A,weigh)%传函分子

h1=freqz(b,1,M)%幅频特性

figure(2)

f=(0:1/M:1-1/M)*fs/2

subplot(211)plot(f,abs(h1))gridtitle('带通')

x2=filter(b,1,x1)

S1=fft(x2)

f=(0:1/N:1/2-1/N)*fs

subplot(212)plot(f,abs(S1(1:N/2)))gridxlabel('胡滚hz')%处理后频谱

function y=myditfft(x)

%本程序对输入序列实现DIT-FFT基2算法，点数取大于等于长度的2的幂次

%------------------------------------

%

myditfft.c

%------------------------------------

m=nextpow2(x)

%求的x长度对应的2的最低幂次m

N=2^m

if length(x)<N

x=[x,zeros(1,N-length(x))]

%若的长度不是2的幂，补0到2的整数幂

end

nxd=bin2dec(fliplr(dec2bin([1:N]-1,m)))+1

%求1：2^m数列的倒序

y=x(nxd)

%将倒序排列作为的初始值

for mm=1:m

%将DFT做m次基2分解，从左到右，对每次分解作DFT运算乎耐

Nmr=2^mm

u=1

%旋转因子u初始化

WN=exp(-i*2*pi/Nmr)

%本次分解的基本DFT因子WN＝exp(-i*2*pi/Nmr)

for j=1:Nmr/2

%本次跨越间隔内的各次碟形运算

for k=j:Nmr:N

%本次碟形运算的岁租春跨型租越间隔为Nmr=2^mm

kp=k+Nmr/2

%确定碟形运算的对应单元下标

t=y(kp)*u

%碟形运算的乘积项

y(kp)=y(k)-t

%碟形运算的加法项

y(k)=y(k)+t

end

u=u*WN

%修改旋转因子，多乘一个基本DFT因子WN

end

end

y(1)和y(2)是分离出来的实部知漏和虚部信号QPSK信号的包络就是y(1)和y(2)的模求出这个QPSK信号包络的均值和方差均值应该是信号的能量，也就是1方差可以通过计算加噪后的能量，然后减去信搭桐烂号的能量得轮指到a=1/sqrt(2)b=-1/sqrt(2)ser=zeros(1,11

---