(完整word版)MATLAB轴心轨迹图、幅值谱、相位谱、功率谱(转)

1、MATLAB 编程，题目如下 x=6*si n( 60*pi*t)+8*si n(120*pi*t+pi/6); y=8*si n( 60*pi*t+pi/2)+3*si n(120*pi*t+pi/4); 程序一：两个波形图和轴心图 clear; close all; sf=800;%采样频率； N=512;%采样点数 t=O:1/sf:1; x=6*si n( 60*pi*t)+8*si n(120*pi*t+pi/6); y=8*si n( 60*pi*t+pi/2)+3*si n(120*pi*t+pi/4); subplot(2,1,1); % 生成 m*n个图，当前激活第一个图 p

2、lot(t,x); xlabel( 时间 t 单位 s); % x轴注解 ylabel( 幅值 x单位 cm); % y 轴注解 title(A 的波形图);%图形标题 grid on; % 显示格线 subplot(2,1,2); plot(t, y, g); xlabel( 时间 t 单位 s); % x轴注解 ylabel( 幅值 y 单位 cm ); % y 轴注解 title(B 的波形图);%图形标题 grid on; % 显示格线 %由心轨迹图 figure(2) plot(x,y,b) grid on; 程序二：计算 X、Y 的波形系数和波峰系数 %x y的波形系数和波峰系数

3、clc; close all; sf=800; t=0:1/sf:1; x=6*si n( 60*pi*t)+8*si n(120*pi*t+pi/6); xav=mean(abs(x);% 平均绝对值 xp=max(x);%幅值也就是最大值 z=x.A2; d=mea n( z); xrms=sqrtm(d);% 有效值 Ft=xrms/xav; Fc=xp/xrms; disp(x 的波形系数是); disp( num2str(Ft); disp(x 的波峰系数是); disp( nu m2str(Fc); %丫的波形系数和波峰系数 y=8*si n(60*pi*t+pi/2)+3*si

4、n(120*pi*t+pi/4); x=y; xav=mean(abs(x);% 平均绝对值 xp=max(x);%幅值也就是最大值 z=x.A2; d=mea n( z); xrms=sqrtm(d);% 有效值 Ft=xrms/xav; Fc=xp/xrms; disp(y 的波形系数是); disp( num2str(Ft); disp(y 的波峰系数是); disp( nu m2str(Fc); 程序三：求解幅频谱、相频谱、功率谱 %畐频谱、相频谱、功率谱 clc; fs=800;%设定采样频率 N=512;%采样点数 n=0:N-1; t=n/fs; y=8*si n( 60*pi*

5、t+pi/2)+3*si n(120*pi*t+pi/4); %进行 FFT 变换并做幅频谱 yy=fft(y,N);% 进行 fft 变换 py仁 yy/512; py2=py1*2;%单边谱乘以 2 mag= 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).A2; subplot(2,1,2); plot(f,power); grid on; tit

6、le(B 的功率谱); xlabel(频率(Hz); ylabel(功率); %相频图 pyy=a ngle(yy); Pyy=180/pi*pyy;% 弧度变角度 figure(2); plot(f,Pyy(1:257); grido n; title(B 的相频谱); xlabel(频率(Hz); ylabel(相位(度); 程序四：嗨目关系数和相干性 clc; closeall; dt=1/800; sf=800; t=0:1/800:1; x=6*si n(60*pi*t)+8*si n(120*pi*t+pi/6); y=8*s in (60*pi*t+pi/2)+3*si n(12

7、0*pi*t+pi/4); %相关系数 N=128; a=xcorr(x,y,128); a=a(N+1:2*N+1); plot(t(1:N),a(1:N); grido n; %求 A、B 的相干函数 figure(2); N=512; cxy,w=cohere(x,y,N/4,sf,ham min g(N/4),0,mea n); plot(w,cxy); grid; title( 相干分析); xlabel(频率(Hz); 以下是 matlab help 的例子 t = 0:0.001:0.6; x = sin (2*pi*50*t)+si n(2*pi*120*t); y = x + 2*ra ndn( size(t); plot(1000*t(1:50),y(1:50) title(Sig nal Corrupted with Zero-Mea n Ra ndom Noise) xlabel(time (millisec on ds) Y = fft(y,512);