怎样用MATLAB绘制信号的幅度谱
== TA α n=0:50; %定义序列的长度是 50 A=1; %设置信号有关的参数 a=0.4; T=1; %采样率 w0=2.0734; x=A*exp(-a*n*T).*sin(w0*n*T); %pi 是 MATLAB 定义的π，信号乘可采用“ .*” close all %清除已经绘制的 x(n)图形 subplot(3,1,1);stem(x); %绘制 x(n)的图形 title(‘理想采样信号序列 ’); k=-25:25; W=(pi/12.5)*k; X=x*(exp(-j*pi/12.5)).^(n’*k); magX=abs(X); %绘制 x(n)的幅度谱 subplot(3,1,2);stem(magX);title(‘理想采样信号序列的幅度谱 ’); angX=angle(X); %绘制 x(n)的相位谱 subplot(3,1,3);stem(angX) ; title (‘理想采样信号序列的相位谱 ’) 2、单位脉冲序列 在 MatLab 中，这一函数可以用 zeros 函数实现： n=1:50; %定义序列的长度是 50 x=zeros(1,50); %注意： MATLAB 中数组下标从 1 开始 x(1)=1; close all; subplot(3,1,1);stem(x);title(‘单位冲击信号序列 ’); k=-25:25; X=x*(exp(-j*pi/12.5)).^(n’*k); magX=abs(X); %绘制 x(n)的幅度谱 subplot(3,1,2);stem(magX);title(‘单位冲击信号的幅度谱 ’); angX=angle(X); %绘制 x(n)的相位谱 subplot(3,1,3);stem(angX) ; title (‘单位冲击信号的相位谱 ’) 3、矩形序列 n=1:50 x=sign(sign(10-n)+1); close all; subplot(3,1,1);stem(x);title(‘单位冲击信号序列 ’); k=-25:25; X=x*(exp(-j*pi/25)).^(n’*k); magX=abs(X); %绘制 x(n)的幅度谱 subplot(3,1,2);stem(magX);title(‘单位冲击信号的幅度谱 ’); angX=angle(X); %绘制 x(n)的相位谱 subplot(3,1,3);stem(angX) ; title (‘单位冲击信号的相位谱 ’) 4、特定冲击串： )3()2(5.2)1(5.2)()( ?+?+?+= nnnnnx δδδδ n=1:50; %定义序列的长度是 50 x=zeros(1,50); %注意： MATLAB 中数组下标从 1 开始 x(1)=1;x(2)=2.5;x(3)=2.5;x(4)=1; close all; subplot(3,1,1);stem(x);title(‘单位冲击信号序列 ’); k=-25:25; X=x*(exp(-j*pi/12.5)).^(n’*k); magX=abs(X); %绘制 x(n)的幅度谱 subplot(3,1,2);stem(magX);title(‘单位冲击信号的幅度谱 ’); angX=angle(X); %绘制 x(n)的相位谱 subplot(3,1,3);stem(angX) ; title (‘单位冲击信号的相位谱 ’) 5、卷积计算： ∑ +∞ ?∞= ?=?= m mnhmxnhnxny )()()()()( 在 MATLAB 中。提供了卷积函数 conv，即 y=conv(x,h)，调用十分方便。例如： 系统： )3()2(5.2)1(5.2)()( ?+?+?+= nnnnnh b δδδδ 信号： 500),sin()( 0 <≤?= ? nnTAetx nT a α n=1:50; %定义序列的长度是 50 hb=zeros(1,50); %注意： MATLAB 中数组下标从 1 开始 hb(1)=1;hb(2)=2.5;hb(3)=2.5;hb(4)=1; close all; subplot(3,1,1);stem(hb);title(‘系统 hb[n]’); m=1:50; %定义序列的长度是 50 A=444.128; %设置信号有关的参数 a=50*sqrt(2.0)*pi; T=0.001; %采样率 w0=50*sqrt(2.0)*pi; x=A*exp(-a*m*T).*sin(w0*m*T); %pi 是 MATLAB 定义的π，信号乘可采用“ .*” subplot(3,1,2);stem(x);title(‘输入信号 x[n]’); y=conv(x,hb); subplot(3,1,3);stem(y);title(‘输出信号 y[n]’); 6、卷积定律验证 k=-25:25; X=x*(exp(-j*pi/12.5)).^(n’*k); magX=abs(X); %绘制 x(n)的幅度谱 subplot(3,2,1);stem(magX);title(‘输入信号的幅度谱 ’); angX=angle(X); %绘制 x(n)的相位谱 subplot(3,2,2);stem(angX) ; title (‘输入信号的相位谱 ’) Hb=hb*(exp(-j*pi/12.5)).^(n’*k); magHb=abs(Hb); %绘制 hb(n)的幅度谱 subplot(3,2,3);stem(magHb);title(‘系统响应的幅度谱 ’); angHb=angle(Hb); %绘制 hb(n)的相位谱 subplot(3,2,4);stem(angHb) ; title (‘系统响应的相位谱 ’) n=1:99; k=1:99; Y=y*(exp(-j*pi/12.5)).^(n’*k); magY=abs(Y); %绘制 y(n)的幅度谱 subplot(3,2,5);stem(magY);title(‘输出信号的幅度谱 ’); angY=angle(Y); %绘制 y(n)的相位谱 subplot(3,2,6);stem(angY) ; title (‘输出信号的相位谱 ’) %以下将验证的结果显示 XHb=X.*Hb; Subplot(2,1,1);stem(abs(XHb));title(‘x(n)的幅度谱与 hb(n)幅度谱相乘 ’); Subplot(2,1,2);stem(abs(Y);title(‘y(n)的幅度谱 ’); axis([0,60,0,8000])绛旓細瑕缁樺埗E=1锛宼=1锛孴=2鐨勫懆鏈熺煩褰㈣剦鍐淇″彿鐨涓夎鍑芥暟褰㈠紡鐨勫箙搴﹂璋卞浘锛屼綘鍙互浣跨敤MATLAB鐨凢FT锛團ast Fourier Transform锛夊嚱鏁般備互涓嬫槸涓涓彲鑳界殑绋嬪簭锛歮atlab澶嶅埗浠ｇ爜 鍙傛暟瀹氫箟 E = 1; % 骞呭害 t = 1; % 鏃堕棿鍋忕Щ T = 2; % 鍛ㄦ湡 Fs = 1000; % 閲囨牱棰戠巼 t_sample = 0:1/Fs:T; % 鏃堕棿...
绛旓細鍦∕ATLAB涓粯鍒朵俊鍙娉㈠舰鍥鹃氬父娑夊強浠ヤ笅姝ラ锛氱敓鎴愭垨瀵煎叆鏁版嵁锛氶鍏堬紝浣犻渶瑕佹湁淇″彿鐨鏁版嵁銆傝繖浜涙暟鎹彲浠ユ槸瀹炴椂閲囬泦鐨勶紝涔熷彲浠ユ槸浠庢枃浠朵腑瀵煎叆鐨勩傚鏋滀綘鏈変俊鍙锋暟鎹枃浠讹紝鍙互浣跨敤MATLAB鐨勫鍏ュ悜瀵兼垨浣跨敤鍑芥暟濡俽eadtable銆乺eadcsv銆乺ead鍌呴噷鍙跺彉鎹㈢瓑鏉ュ姞杞芥暟鎹傞澶勭悊鏁版嵁锛氭牴鎹渶瑕侊紝浣犲彲鑳介渶瑕佸鏁版嵁杩涜棰勫鐞...
绛旓細matlab 瀹氫箟淇″彿鐨鍙傛暟 fs = 1000; % 閲囨牱鐜囦负1000Hz f_cutoff = 100; % 鎴棰戠巼涓100Hz 鍒涘缓婊ゆ尝鍣 N = 1000; % FIR婊ゆ尝鍣ㄧ殑闃舵暟 b = fir1(N, f_cutoff/(fs/2)); % 鐢熸垚FIR婊ゆ尝鍣ㄧ郴鏁 璁＄畻婊ゆ尝鍣鐨勫箙棰戝搷搴 freqz(b, 1, 1024, fs);璁惧畾鍥惧舰鐨勬爣棰樺拰鍧愭爣杞存爣绛 tit...
绛旓細Y=fft(y,512);F =10*f*[0:256]/512;fp=2*sqrt(Y.*conj(Y));锛呭箙搴﹁氨 xp=angle(Y); 锛呯浉浣嶈氨 gl=abs(Y).^2; %鍔熺巼璋 magif=ifft2(abs(f2));%骞呭害閲嶆瀯 pha=angle(f2);%鍙栫浉浣 phaif=ifft2(exp(j*pha));%鐩镐綅閲嶆瀯 鏁版嵁鏄痻(i)锛屽叡N涓偣锛岄噰鏍烽鐜囨槸fsample dt = 1...
绛旓細缁樺埗淇″彿x(n)鐨勫箙搴﹁氨鍜岀浉浣嶈氨锛歬=-25锛25銆乄=(pi/12.5)*k銆俋=x*(exp(-j*pi/12.5)).^(n鈥*k)銆俶agX=abs(X); %缁樺埗x(n)鐨勫箙搴﹁氨锛宻ubplot(3,1,2);stem(magX);title(鈥樼悊鎯抽噰鏍蜂俊鍙峰簭鍒楃殑骞呭害璋扁)銆俛ngX=angle(X); %缁樺埗x(n)鐨勭浉浣嶈氨銆俿ubplot(3,1,3);stem(angX)...
绛旓細缁樺埗淇″彿x(n)鐨勫箙搴﹁氨鍜岀浉浣嶈氨锛歬=-25锛25銆乄=(pi/12.5)*k銆俋=x*(exp(-j*pi/12.5)).^(n鈥*k)銆俶agX=abs(X); %缁樺埗x(n)鐨勫箙4102搴﹁氨锛宻ubplot(3,1,2);stem(magX);title(鈥樼悊鎯抽噰鏍蜂俊鍙峰簭鍒楃殑骞呭害璋扁)銆
绛旓細y=fft(x,N); %瀵淇″彿杩涜蹇烣ourier鍙樻崲 mag=abs(y)*2/N; %姹傚彇Fourier鍙樻崲鐨勬尟骞;*2/N杞彉涓虹湡瀹炲箙鍊 f=n*Fs/N;subplot(3,1,2)plot(f(1:N/2),mag(1:N/2)); %缁樺嚭Nyquist棰戠巼涔嬪墠闅忛鐜囧彉鍖栫殑鎸箙 xlabel('棰戠巼/Hz');ylabel('鎸箙');title('骞呭艰氨');grid on;p=...
绛旓細t=-10*pi:0.1:10*pi;y=sin(t)./t;plot(t,y);xp=angle(Y); 锛呯浉浣嶈氨 gl=abs(Y).^2; %鍔熺巼璋 magif=ifft2(abs(f2));%骞呭害閲嶆瀯 pha=angle(f2);%鍙栫浉浣 phaif=ifft2(exp(j*pha));%鐩镐綅閲嶆瀯
绛旓細50 A=1; %璁剧疆淇″彿鏈夊叧鐨勫弬鏁 a=0.4; T=1; %閲囨牱鐜 w0=2.0734; x=A*exp(-a*n*T).*sin(w0*n*T); %pi 鏄 MATLAB 瀹氫箟鐨勏锛屼俊鍙蜂箻鍙噰鐢ㄢ .*鈥 close all %娓呴櫎宸茬粡缁樺埗鐨 x(n)鍥惧舰 subplot(3,1,1);stem(x); %缁樺埗 x(n)鐨勫浘褰 title(鈥樼悊鎯抽噰鏍蜂俊鍙峰簭鍒 鈥);...
绛旓細鏈夎姹傦紝灏辨槸妯潗鏍囧拰绾靛潗鏍囧悜閲忕殑闀垮害瑕佷竴鑷达紝杩欐牱鎵嶈兘瀵瑰簲鐨勭敾鍥俱傛瘮濡傛椂闂寸偣鏈100涓紝閭ｄ箞瀵瑰簲鐨勫箙鍊肩偣涔熻鏈100涓傚惁鍒欙紝灏变細鍑虹幇浣犺鐨勯敊璇傚笇鏈涜兘澶熷府鍒版ゼ涓粇~~
