数字信号处理实验

于 2023-06-06 16:33:53 发布
阅读量716 收藏 3
点赞数
文章标签： 信号处理
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_40490150/article/details/131070743
版权

1.模拟信号的数字处理

clc;clear up;close all;

f0=13;

t=0:0.002:1;

x=cos(2*pi*f0*t);

subplot(2,2,1);

plot(t,x);

f1=10;

T1=1/f1;

nT1=0:T1:1;

x1=cos(2*pi*f0*nT1);

n=0:length(nT1)-1;

subplot(2,2,2);

stem(n,x1);

clc;clear up;close all;

f0=13;

t=0:0.002:1;

x=cos(2*pi*f0*t);

subplot(3,1,1);

plot(t,x);

f1=10;

T1=1/f1;

nT1=(0:T1:1)';

x1=cos(2*pi*f0*nT1);

n=0:length(nT1)-1;

subplot(3,1,2);

stem(n,x1);

t1 = linspace(-0.5,1.5,500)';

ya=sinc((1/T1)*t1(:,ones(size(nT1)))-(1/T1)*nT1(:,ones(size(t1)))')*x1;

subplot(3,1,3);

plot(t1,ya);

clc;

t=(0:0.002:1)';

xa=cos(80*pi*t)+sin(100*pi*t);

subplot(3,1,1);

plot(t,xa);

fs=200;

T=1/fs;

nT=(0:T:1)';

n=0:length(nT)-1;

x=cos(80*pi*nT)+sin(100*pi*nT);

subplot(3,1,2);

stem(n,x);

ya=sinc((1/T)*t(:,ones(size(nT))')-(1/T)*nT(:,ones(size(t)))')*x;

subplot(3,1,3)

plot(t,ya);

clc;

t=(0:0.002:1)';

xa=cos(80*pi*t)+sin(100*pi*t);

subplot(3,1,1);

plot(t,xa);

fs=150;

T=1/fs;

nT=(0:T:1)';

n=0:length(nT)-1;

x=cos(80*pi*nT)+sin(100*pi*nT);

subplot(3,1,2);

stem(n,x);

tic

% for i=1:size(t)

%     for j=1:size(nT)

%         a(j)=t(i)-nT(j);

%     end

%     ya(i)=sinc((1/T)*a)*x;

% end

% toc

% ya=zeros(1,100);

% n=1;

% while n<=fs+1

%     ya=ya+x(n)*sinc((t-nT(n))/T);

%     n=n+1;

% end

ya=sinc((1/T)*t(:,ones(size(nT))')-(1/T)*nT(:,ones(size(t)))')*x;

subplot(3,1,3)

plot(t,ya);

2离散信号的频域分析

clc;

x=[1,2,3,4];

w=linspace(-pi,pi,500);

[H,w]=freqz(x,1,w);

subplot(2,1,1);

plot(w/pi,abs(H));

subplot(2,1,2);

plot(w/pi,angle(H));

clc;

x=[1,2,3,4];

%DTFT

H=freqz(x,1,8,'whole');

figure;

subplot(2,1,1);

plot(0:7,abs(H));

subplot(2,1,2);

plot(0:7,angle(H));

%FFT

figure;

H1=fft(x,8);

subplot(2,1,1);

plot(0:7,abs(H1));

subplot(2,1,2);

plot(0:7,angle(H1));

clc;

fs=60;

T=1/fs;

nT=(0:T:1)';

n=0:length(nT)-1;

x=cos(80*pi*nT)+sin(100*pi*nT)+exp(1i*150*pi*nT);

y=fft(x);

plot(n,abs(y));

%若为10点DFT

% y=fft(x,10);

% plot(0:9,abs(y));

x1 = [1, 2, 3, 4, 2, 3, 4, 1];

x2 = [4, 3, 2, 1, 1, 2, 3, 2];

X1=fft(x1,8);

X2=fft(x2,8);

Y=ifft(X1.*X2);

C=cconv(x1,x2,8);

isequal(round(Y),C)

  1. 离散时间系统的频域分析

close all;

clear up;

clc;

fs=2000;

T=1/fs;

nT=(0:T:1)';

n=0:length(nT)-1;

x1=[1 1 1 1];

x2=[1 2 3 4 4 3 2 1];

x3=[4 3 2 1 1 2 3 4];

x4=cos(pi/16*n);

x5=sin(pi/32*n);

x6=cos(1000*pi*nT)+cos(1200*pi*nT)+cos(2000*pi*nT);

x7=sin(100*pi*nT)+sin(200*pi*nT)+sin(300*pi*nT);

y1=fft(x1,4);

subplot(3,4,1);

stem(0:3,abs(y1));

subplot(3,4,2);

stem(0:3,angle(y1));

y2=fft(x2,8);

subplot(3,4,3);

stem(0:7,abs(y2));

subplot(3,4,4);

stem(0:7,angle(y2));

y3=fft(x3,8);

subplot(3,4,5);

stem(0:7,abs(y3));

subplot(3,4,6);

stem(0:7,angle(y3));

y4=fft(x4,8);

subplot(3,4,7);

stem(0:7,abs(y4));

subplot(3,4,8);

stem(0:7,angle(y4));

y5=fft(x5,8);

subplot(3,4,7);

stem(0:7,abs(y5));

subplot(3,4,8);

stem(0:7,angle(y5));

y6=fft(x5,8);

subplot(3,4,9);

stem(0:7,abs(y6));

subplot(3,4,10);

stem(0:7,angle(y6));

close all;

clear up;

clc;

fs=2000;

T=1/fs;

nT=(0:T:1)';

n=0:length(nT)-1;

x4=cos(pi/16*n);

x5=sin(pi/32*n);

x8=x4+x5;

y1=fft(x8,32);

subplot(3,4,1);

stem(0:31,abs(y1));

subplot(3,4,2);

stem(0:31,angle(y1));

close all;

clear up;

clc;

fs=2000;

T=1/fs;

nT=0:T:1;

n=0:length(nT)-1;

x1=[1 1 1 1];

x2=[1 2 3 4 4 3 2 1];

x3=[4 3 2 1 1 2 3 4];

x4=cos(pi/16*n);

x5=sin(pi/32*n);

x6=cos(1000*pi*nT)+cos(1200*pi*nT)+cos(2000*pi*nT);

x7=sin(100*pi*nT)+sin(200*pi*nT)+sin(300*pi*nT);

x8=x4+x5;

y6=fft(x6,555);

y7=fft(x7,555);

x = x6 + 1i * x7;  %构造信号

%% 计算16点x(n)的FFT

f = fft(x,555);

%% 利用对称性质求出x6(n)和x7(n)的傅里叶变换

N=555;

kk = conj(f); %取X的共轭

X6 = 0.5*(f+[kk(1),flip(kk(2:555))]);

X7 = -1i*0.5*(f-[kk(1),flip(kk(2:555))]);

figure;

plot((0:554), abs(X6),'r'); %取模，表示信号的频域幅值

figure;

plot((0:554), abs(y6),'r'); %取模，表示信号的频域幅值

4.FFT 的计算机实现

用三种不同的 DFT 程序计算 x(n) = R8 (n) 的傅里叶变换 X m ,并比较三种程序计算机运

行时间。

  1. 用 for 语句的 M 函数文件 dft1.m，用循环变量逐点计算 X(k)；

dft1.m

 function[Am,pha]=dft1(x)

N=length(x);

w=exp(-j*2*pi/N);

for k=1:N

sum=0;

for n=1:N

sum=sum+x(n)*w^((k-1)*(n-1));

end

Am(k)=abs(sum);

pha(k)=angle(sum);

end

sl1.m

N=256;

x=[ones(1,8),zeros(1,N-8)];

t=cputime;

[Am1,pha1]=dft1(x);

t1=cputime-t

n=[0:(length(x)-1)];

w=(2*pi/length(x))*n;

figure(1)

subplot(2,1,1), plot(w,Am1,'b'); grid;

title('Magnitude part');

xlabel('frequency in radians');

ylabel('|X(exp(jw))|');

subplot(2,1,2), plot(w,pha1,'r'); grid;

title('Phase Part');

xlabel('frequency in radians');

ylabel('argX[exp(jw)]/radians');

  1. 编写时域分解FFT 算法或频域分解FFT 算法的M函数文件dft2.m；

function [Am,pha]=dft2(x)

N=length(x); %输入数据的长度

M=log2(N); %蝶形运算的级数

indexs=zeros(1,N); %记录反序码值的数组

for n=0:N-1

    opposite=zeros(1,M); %记录反序码的数组

value=0; %记录反序码的值

    for m=1:M

        s=floor(n/2^(m-1)); %向下取整

        opposite(m)=mod(s,2); %计算余数

        value=value+opposite(m)*2^(M-m);

    end

    indexs(n+1)=value+1;

end

%%蝶形运算

%%W^J*(2^M-L)

%%M=log2(N)

%%J=0,1,2...(2^(L-1))

%%L=1,2 ,3...M

sz=zeros(1,N); %定义输出数组

for n=1:N

   sz(n)=x(indexs(n));

end

for L=1:M %按照 蝶形运行的级数依次计算

    d=2^(L-1); %两个点的距离

    for J=0:d-1

        p=J*(2^(M-L));

        W=exp(-i*2*pi*p/N); %转换因子公式，其中i代表虚部

        for K=J+1:2^L:N %找到每一个蝶形的左上点

            t1=sz(K);

            t2=sz(K+d);

            sz(K)=t1+t2*W; %蝶形运算公式

            sz(K+d)=t1-t2*W;

        end

    end

    Am=abs(sz);

    pha=angle(sz);

end

Xh.m

clc

close

clear

fs=2300;

N=4096;

n=0:N-1;

x=cos(1100*pi*n/fs)+cos(1300*pi*n/fs)+cos(2000*pi*n/fs);

t=cputime;

[Am,pha]=dft2(x);

t2=cputime-t

n=[0:(length(x)-1)];

w=(2*pi/length(x))*n;

figure(2)

subplot(211);

plot(w,Am,'b');grid;

title('Magnitude part');

xlabel('frequency in radians');

ylabel('|X(exp(jw))|');

subplot(212);

plot(w,pha,'r');grid;         %绘图

title('Phase Part');

xlabel('frequency in radians');

ylabel('argX[exp(jw)]/radians');

3、调用FFT库函数，直接计算X(k)；

Dft3.m

function[Am,pha]=dft3(x)

Xk=fft(x);

Am=abs(Xk);

pha=angle(Xk);

sl2.m

N=256;

x=[ones(1,8),zeros(1,N-8)];

t=cputime;

[Am3,pha3]=dft3(x);

t3=cputime-t

n=[0:(length(x)-1)];

w=(2*pi/length(x))*n;

figure(3)

subplot(2,1,1), plot(w,Am3,'b'); grid;

title('Magnitude part');

xlabel('frequency in radians');

ylabel('|X(exp(jw))|');

subplot(2,1,2), plot(w,pha3,'r'); grid;

title('Phase Part');

xlabel('frequency in radians');

ylabel('argX[exp(jw)]/radians');

设计IIR 数字滤波器

clc;clear;close all;

T=0.2*pi;Fs=1/T;

wpz=0.2;    

wsz=0.3;    

rp = 1;       

rs = 15;       

wp=2/T*tan(wpz*pi/2);ws=2/T*tan(wsz*pi/2);

[N, wc] = buttord(wp, ws, rp, rs, 's');

[B, A] = butter(N, wc, 'low','s');

[Bz, Az] = bilinear(B, A,Fs);  

[H, w] = freqz(Bz, Az,25,pi);

wp1 = 0.2;

ws1 = 0.3;

Rp1 = 2;

Rs1 = 30;

WP=2/T*tan(wp1*pi/2);WS=2/T*tan(ws1*pi/2);

[n, wc] = cheb1ord(WP,WS, Rp1, Rs1, 's');

[b, a] = cheby1(n, Rp1, wc,'low', 's');

[Bz1, Az1] = bilinear(b, a,Fs);  

[H1, w1] = freqz(Bz1, Az1,25,pi);

figure

subplot(2,2,1)

plot(w/pi,20*log10(abs(H)));

xlabel('\omega/pi');

ylabel('幅度/(dB)');

title('巴特沃斯低通IIR数字滤波器幅频响应特性曲线');

grid on;

subplot(2,2,2)

plot(w/pi,angle(H));

xlabel('\omega/pi');

ylabel('相位');

title('巴特沃斯低通IIR数字滤波器幅频响应特性曲线');

grid on;

subplot(2,2,3)

plot(w1/pi,20*log10(abs(H1)));

xlabel('\omega/pi');

ylabel('幅度/(dB)');

title('切比雪夫低通IIR数字滤波器幅频响应特性曲线');

grid on;

subplot(2,2,4)

plot(w1/pi,angle(H1));

xlabel('\omega/pi');

ylabel('相位');

title('切比雪夫低通IIR数字滤波器幅频响应特性曲线');

grid on;

clc;clear;close all;

T=0.2*pi;Fs=1/T;

wpz=0.2;    

wsz=0.3;    

  • rp = 1;       

rs = 15;       

wp=2/T*tan(wpz*pi/2);ws=2/T*tan(wsz*pi/2);

[N, wc] = buttord(wp, ws, rp, rs, 's');

[B, A] = butter(N, wc, 'low','s');

[Bz, Az] = bilinear(B, A,Fs);  

[H, w] = freqz(Bz, Az,25,pi);

wp1 = 0.2;

ws1 = 0.3;

Rp1 = 2;

Rs1 = 30;

WP=2/T*tan(wp1*pi/2);WS=2/T*tan(ws1*pi/2);

[n, wc] = cheb1ord(WP,WS, Rp1, Rs1, 's');

[b, a] = cheby1(n, Rp1, wc,'low', 's');

[Bz1, Az1] = bilinear(b, a,Fs);  

[H1, w1] = freqz(Bz1, Az1,25,pi);

x=[-4,-2,0,-4,-6,-4,-2,-4,-6,-6,-4,-4,-6,-6,-2,6,12,8,0,-16,-38,...

-60,-84,-90,-66,-32,-4,-2,-4,8,12,12,10,6,6,6,4,0,0,0,0,0,-2,...

-4,0,0,0,-2,-2,0,0,-2,-2,-2,-2,0];

 b= filter(Bz,Az, x);

 q= filter(Bz1,Az1,x);

n = 1:length(x);

subplot(3,1,1);

plot(n, x);

title('原始信号');

xlabel('采样点');

ylabel('幅度');

subplot(3,1,2);

plot(n, b);

title('巴特沃斯滤波后的信号');

xlabel('采样点');

ylabel('幅度');

subplot(3,1,3);

plot(n, q);

title('切比雪夫滤波后的信号');

xlabel('采样点');

ylabel('幅度');

clc;clear;close all;

Fs=44100;

wpz=0.2;    

wsz=0.3;    

rp = 1;       

rs = 15;       

wp=2*Fs*tan(wpz*pi/2);ws=2*Fs*tan(wsz*pi/2);

[N, wc] = buttord(wp, ws, rp, rs, 's');

[B, A] = butter(N, wc, 'low','s');

[Bz, Az] = bilinear(B, A,Fs);  

[H, w] = freqz(Bz, Az,25,pi);

wp1 = 0.2;

ws1 = 0.3;

Rp1 = 2;

Rs1 = 30;

WP=2*Fs*tan(wp1*pi/2);WS=2*Fs*tan(ws1*pi/2);

[n, wc] = cheb1ord(WP,WS, Rp1, Rs1, 's');

[b, a] = cheby1(n, Rp1, wc,'low', 's');

[Bz1, Az1] = bilinear(b, a,Fs);  

[H1, w1] = freqz(Bz1, Az1,25,pi);

load('xas.mat');

 b= filter(Bz,Az, xas);

  q= filter(Bz1,Az1,xas);

n = 1:length(xas);

subplot(3,1,1);

plot(n, real(xas));

%sound( real(xas),Fs);

title('原始语言信号');

xlabel('采样点');

ylabel('幅度');

subplot(3,1,2);

plot(n, real(b));

sound( real(b),Fs);

title('巴特沃斯滤波后的语言信号');

xlabel('采样点');

ylabel('幅度');

subplot(3,1,3);

plot(n, real(q));

title('切比雪夫滤波后的语言信号');

xlabel('采样点');

ylabel('幅度');

%sound( real(q),Fs);

用窗函数法设计 FIR 数字滤波器

% 产生窗函数并画图

N = 64; % 窗长度

n = 0:N-1;

w_rect = ones(1,N);

w_hann = 0.5*(1-cos(2*pi*n/(N-1)));

w_hamming = 0.54-0.46*cos(2*pi*n/(N-1));

w_blackman = 0.42-0.5*cos(2*pi*n/(N-1))+0.08*cos(4*pi*n/(N-1));

figure;

plot(n,w_rect,'r',n,w_hann,'g',n,w_hamming,'b',n,w_blackman,'m');

title('Window Functions');

xlabel('Sample Index');

ylabel('Amplitude');

legend('Rectangular','Hanning','Hamming','Blackman');

% 在命令窗口输入以下命令即可生成并绘制矩形窗、汉宁窗、汉明窗和布莱克曼窗函数：

function [w] = window(type, N)

% type: 窗类型，可选 'rect', 'hanning', 'hamming', 'blackman'

% N: 窗长度

    switch type

        case 'rect'

            w = ones(N,1);

        case 'hanning'

            w = 0.5 - 0.5*cos(2*pi*(0:N-1)/(N-1));

        case 'hamming'

            w = 0.54 - 0.46*cos(2*pi*(0:N-1)/(N-1));

        case 'blackman'

            w = 0.42 - 0.5*cos(2*pi*(0:N-1)/(N-1)) + 0.08*cos(4*pi*(0:N-1)/(N-1));

        otherwise

            error('Unsupported window type!');

    end

    plot(w);

end

clc

close all

clear up

load('xas.mat');

% 滤波器参数

a = 0.4;

N1 = 15; % 窗长度为15

N2 = 33; % 窗长度为33

% 理想幅频特性

w = linspace(0, pi, 1000);

Hd = exp(-1j*w*a).*(w <= 0.4*pi);

% 矩形窗设计滤波器

h1_rect = fir1(N1-1, a/pi);

h2_rect = fir1(N2-1, a/pi);

H1_rect = freqz(h1_rect, 1, w);

H2_rect = freqz(h2_rect, 1, w);

[H1,w1]= freqz(h1_rect, 1, w);

[H2,w2]= freqz(h2_rect, 1, w);

% 汉宁窗设计滤波器

w_hann1 = hann(N1)';

w_hann2 = hann(N2)';

h1_hann = fir1(N1-1, a/pi, w_hann1);

h2_hann = fir1(N2-1, a/pi, w_hann2);

H1_hann = freqz(h1_hann, 1, w);

H2_hann = freqz(h2_hann, 1, w);

[H3,w3] = freqz(h1_hann, 1, w);

[H4,w4]= freqz(h2_hann, 1, w);

% 汉明窗设计滤波器

w_hamming1 = hamming(N1)';

w_hamming2 = hamming(N2)';

h1_hamming = fir1(N1-1, a/pi, w_hamming1);

h2_hamming = fir1(N2-1, a/pi, w_hamming2);

 H1_hamming = freqz(h1_hamming, 1, w);

H2_hamming = freqz(h2_hamming, 1, w);

[H5,w5]= freqz(h1_hamming, 1, w);

[H6,w6] = freqz(h2_hamming, 1, w);

%%

%语音信号恢复

Fs=44100;

load('xas.mat');

x = xas;

y = filter(h2_hamming, 1, x);

sound(real(y), Fs);

n = 1:length(xas);

figure(3);

subplot(3,1,1);

plot(n, real(xas));

subplot(3,1,2);

plot(n, real(y));

subplot(3,1,3);

plot(n, abs(fftshift(fft(y))));

%%

% 布莱克曼窗设计滤波器

w_blackman1 = blackman(N1)';

w_blackman2 = blackman(N2)';

h1_blackman = fir1(N1-1, a/pi, w_blackman1);

h2_blackman = fir1(N2-1, a/pi, w_blackman2);

H1_blackman = freqz(h1_blackman, 1, w);

H2_blackman = freqz(h2_blackman, 1, w);

[H7,w7] = freqz(h1_blackman, 1, w);

[H8,w8] = freqz(h2_blackman, 1, w);

% Kaiser窗设计滤波器

beta = 8.5;

w_kaiser1 = kaiser(N1, beta)';

w_kaiser2 = kaiser(N2, beta)';

h1_kaiser = fir1(N1-1, a/pi, w_kaiser1);

h2_kaiser = fir1(N2-1, a/pi, w_kaiser2);

H1_kaiser = freqz(h1_kaiser, 1, w);

H2_kaiser = freqz(h2_kaiser, 1, w);

[H9,w9] = freqz(h2_kaiser, 1, w);

[H10,w10] = freqz(h2_kaiser, 1, w);

% 绘图

figure(1);

subplot(2,2,1);

plot(w1/pi, 20*log10(abs(H1)), 'r', w3/pi, 20*log10(abs(H3)), 'g', w5/pi, 20*log10(abs(H5)), 'b', w7/pi, 20*log10(abs(H7)), 'm', w9/pi, 20*log10(abs(H9)), 'k');

title('损耗函数曲线(N=15)');

xlabel('\omega/pi');

ylabel('幅度/（dB）');

legend('Rectangular', 'Hanning', 'Hamming', 'Blackman', 'Kaiser');

subplot(2,2,2);

plot(w2/pi,20*log10(abs(H1)), 'r', w4/pi, 20*log10(abs(H4)), 'g', w6/pi, 20*log10(abs(H6)), 'b', w8/pi, 20*log10(abs(H8)), 'm', w10/pi, 20*log10(abs(H10)), 'k');

title('损耗函数曲线(N=33)');

xlabel('\omega/pi');

ylabel('幅度/（dB）');

legend('Rectangular', 'Hanning', 'Hamming', 'Blackman', 'Kaiser');

subplot(2,2,3);

plot(w1/pi, angle(H1), 'r', w3/pi, angle(H3), 'g', w5/pi, angle(H5), 'b', w7/pi, angle(H7), 'm', w9/pi, angle(H9), 'k');

title('相频特性曲线(N=15)');

xlabel('w/pi');

ylabel('Phase (rad)');

legend('Rectangular', 'Hanning', 'Hamming', 'Blackman', 'Kaiser');

subplot(2,2,4);

plot(w2/pi, angle(H2), 'r', w4/pi, angle(H4), 'g', w6/pi, angle(H6), 'b', w8/pi, angle(H8), 'm', w10/pi, angle(H10), 'k');

title('相频特性曲线(N=33)');

xlabel('w/pi');

ylabel('Phase (rad)');

legend('Rectangular', 'Hanning', 'Hamming', 'Blackman', 'Kaiser');

figure(2);

subplot(1,2,1);

plot(w/pi, abs(H1_rect), 'r', w/pi, abs(H1_hann), 'g', w/pi, abs(H1_hamming), 'b', w/pi, abs(H1_blackman), 'm', w/pi, abs(H1_kaiser), 'k');

title('损耗函数曲线(N=15)');

xlabel('\omega/pi');

ylabel('幅度/（dB）');

legend('Rectangular', 'Hanning', 'Hamming', 'Blackman', 'Kaiser');

subplot(1,2,2);

plot(w/pi, abs(H2_rect), 'r', w/pi, abs(H2_hann), 'g', w/pi, abs(H2_hamming), 'b', w/pi, abs(H2_blackman), 'm', w/pi, abs(H2_kaiser), 'k');

title('损耗函数曲线(N=33)');

xlabel('\omega/pi');

ylabel('幅度/（dB）');

legend('Rectangular', 'Hanning', 'Hamming', 'Blackman', 'Kaiser');

%%

%语音信号恢复

Fs=44100;

load('xas.mat');

x = xas;

y = filter(h2_hamming, 1, x);

sound(real(y), Fs);

n = 1:length(xas);

subplot(3,1,1);

plot(n, real(xas));

subplot(3,1,2);

plot(n, real(y));

subplot(3,1,3);

plot(n, abs(fftshift(fft(y))));

%%

x=[-4 -2 0 -4 -6 -4 -2 -4 -6 -6 -4 -4 -6 -6 -2 6 12 8 0 -16 -38 -60 -84 -90 66 -32 -4 -2 -4 8 12 12 10 6 6 6 4 0 0 0 0 0 -2 -4 0 0 0 -2 -2 0 0 -2 -2 -2 -2 0];

 b= filter(h2_hamming, 1,x);

n = 1:length(x);

figure(4)

subplot(2,1,1);

plot(n,x);

title('原始信号');

xlabel('采样点');

ylabel('幅度');

subplot(2,1,2);

plot(n, b);

title('汉宁窗设计滤波器滤波后的信号');

%语音信号恢复

Fs=44100;

 load('xas.mat');

 x = xas;

% y = filter(h2_hamming, 1, x);

% sound(real(y), Fs);

% n = 1:length(xas);

% figure(3);

% subplot(3,1,1);

% plot(n, real(xas));

% subplot(3,1,2);

% plot(n, real(y));

% subplot(3,1,3);

% plot(n, abs(fftshift(fft(y))));

c1=filter(h2_rect,1,xas);

d1=filter(h2_hann,1,xas);

e1= filter(h2_hamming, 1,xas);

f1=filter(h2_blackman,1,xas);

g1=filter(h2_kaiser,1,xas);

n = 1:length(x);

figure

subplot(3,4,1);

plot(n,real(xas));

subplot(3,4,2);

plot(n, abs(fftshift(fft(xas))));

title('原始信号');

xlabel('采样点');

ylabel('幅度');

subplot(3,4,3);

plot(n, real(c1));

subplot(3,4,4);

plot(n, abs(fftshift(fft(c1))));

title('Rectangular设计滤波器滤波后的信号');

subplot(3,4,5);

plot(n, real(d1));

subplot(3,4,6);

plot(n, abs(fftshift(fft(d1))));

title('Hanning设计滤波器滤波后的信号');

xlabel('采样点');

ylabel('幅度');

subplot(3,4,7);

plot(n, real(e1));

subplot(3,4,8);

plot(n, abs(fftshift(fft(e1))));

title('Hamming设计滤波器滤波后的信号');

subplot(3,4,9);

plot(n, real(f1));

subplot(3,4,10);

plot(n, abs(fftshift(fft(f1))));

title('Blackman设计滤波器滤波后的信号');

xlabel('采样点');

ylabel('幅度');

subplot(3,4,11);

plot(n, real(g2));

subplot(3,4,12);

plot(n, abs(fftshift(fft(g2))));

title('Kaiser设计滤波器滤波后的信号');

sound(real(g1), Fs);

fp=2*pi*1.5*10^3;fs=2*pi*3*10^3;rs=50;Fs=2*pi*1.5*10^4;

wp=2*pi*fp/Fs;

ws=2*pi*fs/Fs;

Bt=ws-wp;

alph=0.5842*(rs-21)^0.4+0.07886*(rs-21);

M=ceil((rs-8)/2.285/Bt);

wc=(wp+ws)/2/pi;

hn=fir1(M,wc,kaiser(M+1,alph));

%% 绘制幅频响应曲线

[H, w] = freqz(hn, 1, 512, Fs);

figure;

subplot(2,1,1);

plot(w, 20*log10(abs(H)));

title('幅频响应曲线');

xlabel('频率(rad/sec)');

ylabel('幅频(dB)');

%% 绘制相频响应曲线

subplot(2,1,2);

plot(w, angle(H));

title('相频响应曲线');

xlabel('频率(rad/sec)');

ylabel('相位(rad)');

wpl=0.55*pi;

wpu=0.7*pi;

wsl=0.45*pi;

wsu=0.8*pi;rs=40;

wc=[(wpl+wsl)/2,(wpu+wsu)/2];

Bt=wpl-wsl;

N=ceil(6.2*pi/Bt);

hn=fir1(N-1,wc/pi,hanning(N));

subplot(2,2,1);yn='h(n)';

stem(hn);

xlabel('n');

ylabel('h(n)');

subplot(2,2,2);A=1;

[H,w]=freqz(hn);

plot(w/pi,20*log10(abs(H)));

xlabel('\omega/pi');

ylabel('幅度/(dB)');

title('损耗函数曲线');

信*
关注
  • 0
    点赞
  • 踩
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
数字信号处理_实验一_信号、系统及系统响应
weixin_48023487的博客
04-11 3940
数字信号处理试验一：介绍了MATLAB生成基本离散信号，求系统响应，DTFT，模拟采样定理
数字信号处理实验报告.docx
12-01
用MATLAB工具实现数字信号处理实验，实验一 离散时间信号和系统响应；实验二 用FFT对信号作频谱分析；实验三 用双线性变换法设计IIR数字滤波器；实验四 用窗函数法设计FIR数字滤波器；适用于数字信号处理课程学习者...
数字信号处理实验matlab仿真源代码
04-14
ex020100 信号合成 ex020200 信号合成 ex020300 复数序列的信号合成 ex020400 奇偶合成 ex020500 卷积计算 ex020600 卷积的图解 ex020700 卷积计算 ex02070b 卷积计算 ex020800 互相关计算 ex020700 卷积计算 ex020900 解差分方程 ex021000 解差分方程 例3.1~例3.2 求离散付利叶变换 ex030300 例3.1中x(n)=(0.9)n 的频谱曲线绘制 ex030400 用矩阵-向量乘法求有限长序列的DTFT ex030500 x(n)=(exp(jπ/3))n 的频谱及其周期性 ex030600 x(n)=2n 的离散付利叶变换及其共轭对称性 ex030700 DFDT 线性性(3.5)的验证 ex030800 DFDT 时域移位性(3.6)的验证 ex030900 DFDT 频域移位性(3.7)的验证 ex031000 DFDT 共轭性(3.8)的验证 ex031100 DFDT 折叠性(3.9)的验证 ex031200 DFDT 对称性(3.10)的验证 ex031300 脉冲函数为h(n)=(0.9)nu(n) 的系统的频谱曲线 ex031400 系统稳态输出的计算 ex031500 以差分方程表示的系统的频谱函数和稳态输出的计算 ex031600 以差分方程表示的滤波器的频谱函数的计算 ex031700 求付利叶变换及绘制曲线 ex031800 不同采样频率对频谱曲线的影响 ex031900 用例3.18a中的x(n)重构x(t) ex032000 用例3.18b中的x(n)重构x(t) ex032100 用ZOH和FOH把例3.18中的x(n)重构为x(t) ex032200 用spline函数把例3.18中的x1(n)和x2(n)重构为xa(t) 例4.1~例4.3 求z变换(须用MATLAB中的symbolic工具箱) 例4.4~例4.6 用其他方法求z变换 例4.7 求z反变换 ex040800 检验residuez函数 ex040900 求z反变换 ex041000 求不带复数的z反变换 ex04100a 求不带复数的z反变换 ex041100 由差分方程求零-极点及频率响应 例4.12 由离散传递函数求脉冲过度函数 ex041300 由差分方程求系统函数及脉冲响应 ex050200 宽度L周期N的周期性方波的离散付利叶曲线 ex050500 不同的离散付利叶采样密度对应的时域曲线 ex050600 离散付利叶变换计算 ex050700 例5.6取不同周期所得离散付利叶曲线 ex050800 高密度和高分辨频谱的差别 ex050900 循环折叠特性的检验 ex051000 循环奇偶分解特性(5.34)的检验 ex051100 循环移位特性的检验 ex051200 循环移位特性的检验 ex051300 循环卷积的计算 ex051400 循环卷积的计算 ex051500 周期N对循环卷积的影响 ex051600 循环卷积和线性卷积的比较 ex051700 周期N对循环卷积的影响 ex051800 用重叠保留法计算循环卷积 ex051900 用重叠保留法计算循环卷积 ex052000 四点FFT计算 ex052100 1<N<2048点FFT执行时间的比较 ex052200 1<L<150点快速卷积与FFT的执行时间比较 ex060100 级联形式转换 ex060200 并联形式转换 ex060300 混合形式转换 例6.4 线性相位系统的级联形式 例6.5 线性相位系统的无复数级联形式 ex060600 给定h(n),求其频率采样形式 ex060700 求频率采样形式并与线性相位形式相比较 例6.8 由差分方程求格型形式 例6.9 由全极点形式求格型形式 例6.10 由零-极点形式求梯形-格型形式 例7.1~例7.2 滤波器相对指标与绝对指标的转换 ex070300 振幅响应和幅度响应 ex070400 1-型线性相位FIR滤波器 ex070500 2-型线性相位FIR滤波器 ex070600 3-型线性相位FIR滤波器 ex070700 4-型线性相位FIR滤波器 ex070800 低通滤波器设计 - 哈明窗 ex070900 低通滤波器设计 - 凯泽窗 ex071000 带通滤波器设计 - 布莱克曼窗 ex071100 带阻滤波器设计 - 凯泽窗 ex071200 差分器设计 - 哈明窗
数字信号处理实验1
weixin_45662974的博客
07-19 4380
文章目录1.题目一2.题目二3.题目三4.题目四5.题目五 1.题目一 分别用MATLAB和Python画出信号x(n)=，其中分别作图，取a=0.5和a=1.5，绘制图形如下： n=-2:1:20; %x(n)信号的范围 n0=0; a1=0.5; %|a|<1 a2=1.5; %|a|>1 u=stepfun(n,n0); .
数字信号处理实验集合
weixin_44003265的博客
12-02 1万+
实验一 信号、系统及系统响应 一、实验目的 1、 熟悉连续信号经理想采样前后的频谱变化关系，加深对时域采样定理的理解。 2、 熟悉时域离散系统的时域特性。 3、 利用卷积方法观察分析系统的时域特性。 4、 掌握序列傅里叶变换的计算机实现方法，利用序列的傅里叶变换对连续信号、离散信号及系统响应进行频域分析。 二、实验代码 %dsp11 clear all; close all; figure(1) a1 = FF(444.128,222,222,1000); n = 0:49; [b1,K,D] = dt
数字信号处理实验.pdf
07-14
数字信号处理实验.pdf数字信号处理实验.pdf数字信号处理实验.pdf数字信号处理实验.pdf数字信号处理实验.pdf数字信号处理实验.pdf
北京理工大学数字信号处理实验报告.docx
03-22
北京理工大学数字信号处理实验报告
数字信号处理实验_数字信号处理_
10-02
信号的基本运算和变换，传递函数求系统冲激响应和阶跃响应
数字信号处理实验（刘舒帆版）答案
11-27
数字信号处理实验（刘舒帆版）答案。1、本实验教材中所有例题程序均例入此文件夹，程序名称与例题名称均一一对应。如： 例2-1，程序名为MS2_1。 2、个别例题可能不止一段程序，则教材上第1段程序与例题名称均对应，...
数字信号处理（高西全）实验课matlab程序
04-20
四川大学通信工程的专业课《数字信号处理》实验代码，matlab代码。配套书籍为《数字信号处理》高西全版，代码文件按“p123.m”格式命名，标明该段代码是书本第几页的实现代码。可供参考。
数字信号处理实验汇总
05-21
数字信号处理实验汇总，包括 离散信号产生和实现 、离散系统的变换域分析、 FFT算法的应用、无限冲激响应数字滤波器设计、有限冲激响应数字滤波器设计。
数字信号处理实验报告汇总（实验一到实验五）
01-01
数字信号处理实验报告，自己上机自己截图。仅供参考。
数字信号处理实验(一)
热门推荐
PiggyGaGa的博客
08-01 2万+
实验目的本次实验目的为：在matlab环境下产生几种基本的数字信号，并对这些基本的信号进行运算和变换，同时利用程序结果对采样定理进行验证，深刻理解采样定理。通过自己录制音频信号并对不同的音频信号进行不同处理，加深理解音频信号中声道的原理，以及混声、回声的形成原理。实验内容 用matlab产生单位脉冲信号，单位阶跃信号，矩形信号，正弦信号，余弦信号，指数信号，产生并观察f(x)=sinc(x)函数的波
数字信号处理实验2
weixin_45662974的博客
07-19 3738
目录1.题目一2.题目二3.题目三4.题目四 1.题目一 import matplotlib.pyplot as plt import numpy as np n=np.linspace(-4,4,9) #生成-5<n<5内的9个点 y1=np.zeros((1,9)) n=n.tolist() #将array转换为list y1=y1.tolist() y1=y1[0] index1=n.index(-3) index2=n.index(4.
数字信号处理综合实验——Matlab实现DTMF信号的产生与提取
cxcx2001的博客
07-05 8355
（1）通过查阅资料，自学电话中DTMF信号的产生与检测方法。 （2）利用Matlab实现DTMF信号的产生与提取。 自行选择数字（0-9），通过数字方法产生该数字的双频信号；拟定采样频率为8000Hz。根据该信号的频率特性，确定滤波器类型及技术指标，将两个单音分别提取出来。画出数字滤波器输入、输出信号的波形及频谱。 （3）利用Matlab GUI设计用户界面，实现拨号界面，包括10个数字键、显示框（显示按下的数字），画出相应的波形。 （4）自拟实验方案及具体实验步骤（要求写出设计方案）。......
数字信号处理实验6
weixin_45662974的博客
07-19 1万+
目录1.题目一2.题目二3.题目三4.题目四5.题目五6.题目六7.题目七 1.题目一 已知有限长序列x(n) = [4,0,3,0,2,0,1]，求x(n)右移2位成为新的向量y(n)，并画出循环移位的中间过程。(利用matlab作答) 先建立x(n)的序列，将其进行周期延拓后再进行移位操作，最后取主值序列即可得到循环移位后的序列。结果如下： %question 1 clear; xn=[4,0,3,0,2,0,1];%建立xn序列 Nx=length(xn); nx=0:Nx-1; nx1=-.
数字信号处理实验一
weixin_41210618的博客
10-13 3103
补充2019/10/26 (1)序列的加、减、乘、除和乘方运算。输入A=[1 2 3 4]，B=[3 4 5 6]，起点n=0，求C=A+B，D=A-B，E=A.*B，F=A./B，G=A.^B，并用stem语句画出A，B，C，D，E，F，G。(2)用MATLAB实现下列序列： 并分别用EXCLE文件、.txt、.mat格式保存两个序列。（3）编写MATLAB程序，导入上述保存文件数据，在同一个图形窗口中显示两组数据，分别添加标注，添加图的名称，轴名称和网格线。（4）编写一个函数实现专家评分，评分规则：给
现代数字信号处理实验
最新发布
12-13
现代数字信号处理实验通常包括以下内容： 1. 正弦信号的相关积累检测实验：通过MATLAB等工具，生成正弦信号并进行相关积累检测，了解正弦信号的基本特性。 2. 含噪谐波信号的周期图功率谱估计实验：通过MATLAB等工具，生成含噪谐波信号并进行周期图功率谱估计，了解谐波信号的频谱特性。 3. 相关正态分布离散随机过程的产生实验：通过MATLAB等工具，生成相关正态分布离散随机过程并进行分析，了解相关随机过程的基本特性。 4. 确定性最小二乘滤波器的实现实验：通过MATLAB等工具，实现确定性最小二乘滤波器并进行分析，了解滤波器的基本原理和应用。 5. 基于AR模型的谐波信号预测与参数估计实验：通过MATLAB等工具，实现AR模型并进行谐波信号预测和参数估计，了解AR模型的基本原理和应用。 6. 线性信道均衡中的LMS梯度实验：通过MATLAB等工具，实现LMS梯度算法并进行线性信道均衡，了解LMS算法的基本原理和应用。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值