1、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;xlabel(時(shí)間t 單位s; % x軸注解ylab

2、el(幅值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(2plot(x,y,bgrid 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=m

3、ean(abs(x;%平均絕對值xp=max(x;%幅值也就是最大值z=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=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;di

4、sp(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=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;%做

5、幅頻譜圖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(相位（度）;程序四：%相關(guān)系數(shù)和相干性clc;closeall;dt=1/800;sf=800;t=0:1/800:1;x

6、=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;%相關(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:50title(Signal Corrupted with Zero-Mean Random Noise xlabel(time (millisecondsY = fft(y,512;Pyy