MATLAB-軸心軌跡圖、幅值譜、相位譜、功率譜(轉(zhuǎn))(共4頁(yè))

1、精選優(yōu)質(zhì)文檔-傾情為你奉上MATLAB編程，題目如下x=6*sin(60*pi*t)+8*sin(120*pi*t+pi/6);y=8*sin(60*pi*t+pi/2)+3*sin(120*pi*t+pi/4);程序一：兩個(gè)波形圖和軸心圖clear;close all;sf=800;%采樣頻率;N=512;%采樣點(diǎn)數(shù)t=0:1/sf:1;x=6*sin(60*pi*t)+8*sin(120*pi*t+pi/6);y=8*sin(60*pi*t+pi/2)+3*sin(120*pi*t+pi/4);subplot(2,1,1); %生成m*n個(gè)圖，當(dāng)前激活第一個(gè)圖plot(t,x);xlabe

2、l('時(shí)間t 單位s'); % x軸注解ylabel('幅值x 單位cm'); % y軸注解title('A的波形圖'); % 圖形標(biāo)題grid on; % 顯示格線subplot(2,1,2);plot(t, y, 'g');xlabel('時(shí)間t 單位s'); % x軸注解ylabel('幅值y 單位cm '); % y軸注解title('B的波形圖'); % 圖形標(biāo)題grid on; % 顯示格線%軸心軌跡圖figure(2)plot(x,y,'b')grid

3、on;程序二：計(jì)算X、Y的波形系數(shù)和波峰系數(shù)%x y的波形系數(shù)和波峰系數(shù)clc;close all;sf=800;t=0:1/sf:1;x=6*sin(60*pi*t)+8*sin(120*pi*t+pi/6);xav=mean(abs(x);%平均絕對(duì)值xp=max(x);%幅值也就是最大值z(mì)=x.2;d=mean(z);xrms=sqrtm(d);%有效值Ft=xrms/xav;Fc=xp/xrms;disp('x的波形系數(shù)是');disp(num2str(Ft);disp('x的波峰系數(shù)是');disp(num2str(Fc);%y的波形系數(shù)和波峰系數(shù)y=

4、8*sin(60*pi*t+pi/2)+3*sin(120*pi*t+pi/4);x=y;xav=mean(abs(x);%平均絕對(duì)值xp=max(x);%幅值也就是最大值z(mì)=x.2;d=mean(z);xrms=sqrtm(d);%有效值Ft=xrms/xav;Fc=xp/xrms;disp('y的波形系數(shù)是');disp(num2str(Ft);disp('y的波峰系數(shù)是');disp(num2str(Fc);程序三：求解幅頻譜、相頻譜、功率譜%幅頻譜、相頻譜、功率譜clc;fs=800;%設(shè)定采樣頻率N=512;%采樣點(diǎn)數(shù)n=0:N-1;t=n/fs;y=

5、8*sin(60*pi*t+pi/2)+3*sin(120*pi*t+pi/4);%進(jìn)行FFT變換并做幅頻譜yy=fft(y,N);%進(jìn)行fft變換py1=yy/512;py2=py1*2;%單邊譜乘以2mag= abs(py2) ;%求幅值f=(0:256)*800/512;subplot(2,1,1);plot(f,mag(1:257);%做幅頻譜圖xlabel('頻率(Hz)');ylabel('幅值');title('B的幅頻譜');grid;%功率譜power=mag(1:257).2;subplot(2,1,2);plot(f,pow

6、er);grid on;title('B的功率譜');xlabel('頻率(Hz)');ylabel('功率'); %相頻圖pyy=angle(yy);Pyy=180/pi*pyy;%弧度變角度f(wàn)igure(2);plot(f,Pyy(1:257);gridon;title('B的相頻譜');xlabel('頻率(Hz)');ylabel('相位（度）');程序四：%相關(guān)系數(shù)和相干性clc;closeall;dt=1/800;sf=800;t=0:1/800:1;x=6*sin(60*pi*t)+8

7、*sin(120*pi*t+pi/6);y=8*sin(60*pi*t+pi/2)+3*sin(120*pi*t+pi/4);%相關(guān)系數(shù)N=128;a=xcorr(x,y,128);a=a(N+1:2*N+1);plot(t(1:N),a(1:N);gridon;%求A、B的相干函數(shù)figure(2);N=512;cxy,w=cohere(x,y,N/4,sf,hamming(N/4),0,'mean');plot(w,cxy);grid;title('相干分析');xlabel('頻率(Hz)');以下是matlab help的例子t = 0:0.001:0.6; x = sin(2*pi*50*t)+sin(2*pi*120*t); y = x + 2*randn(size(t); plot(1000*t(1:50),y(1:50) title('Signal Corrupted with Zero-Mean Random Noise') xlabel('time (milliseconds)'Y = fft(y,512);Pyy = Y.* conj(Y) /