matlab全面模拟汽轮机组的轴心轨迹图，方便研究使用.rar

MATLAB 轴心轨迹图、幅值谱、相位谱、功率谱 2011-04-22 22:30:33| 分类： 专业相关 | 标签：知识点通信 |字号大中小 订阅 . 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); 程序一：两个波形图和轴心图 clear; close all; sf=800;%采样频率; N=512;%采样点数 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个图，当前激活第一个图 plot(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的波形系数和波峰系数 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));%平均绝对值 xp=max(x);%幅值也就是最大值 z=x.^2; d=mean(z); xrms=sqrtm(d);%有效值 Ft=xrms/xav; Fc=xp/xrms; disp(['x的波形系数是']); disp(num2str(Ft)); disp(['x的波峰系数是']); disp(num2str(Fc)); %y的波形系数和波峰系数 y=8*sin(60*pi*t+pi/2)+3*sin(120*pi*t+pi/4); x=y; xav=mean(abs(x));%平均绝对值 xp=max(x);%幅值也就是最大值 z=x.^2; d=mean(z); xrms=sqrtm(d);%有效值 Ft=xrms/xav; Fc=xp/xrms; disp(['y的波形系数是']); disp(num2str(Ft)); disp(['y的波峰系数是']); disp(num2str(Fc)); 程序三：求解幅频谱、相频谱、功率谱 %幅频谱、相频谱、功率谱 clc; fs=800;%设定采样频率 N=512;%采样点数 n=0:N-1; t=n/fs; y=8*sin(60*pi*t+pi/2)+3*sin(120*pi*t+pi/4); %进行FFT变换并做幅频谱 yy=fft(y,N);%进行fft变换 py1=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).^2; subplot(2,1,2); plot(f,power); grid on; title('B的功率谱'); xlabel('频率(Hz)'); ylabel('功率'); %相频图 pyy=angle(yy); Pyy=180/pi*pyy;%弧度变角度 figure(2); plot(f,Pyy(1:257)); gridon; title('B的相频谱'); xlabel('频率(Hz)'); ylabel('相位（度）'); 程序四：%相关系数和相干性 clc; closeall; dt=1/800; sf=800; t=0:1/800: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); %相关系数 N=128; a=xcorr(x,y,128); a=a(N+1:2*N+1); plot(t(1:N),a(1:N)); gridon; %求A、B的相干函数 figure(2); N=512; [cxy,w]=cohere(x,y,N/4,sf,hamming(N/4),0,'mean'); plot(w,cxy); grid; title('相干分析'); xlabel('频率(Hz)');