MATLAB-轴心轨迹图、幅值谱、相位谱、功率谱(转)

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)；程序1:两个波形和轴图清除；全部关闭；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)；子情节(2，1，1)；生成了% m*n个图表，当前激活了第一个图表图(t，x)；Xlabel(时间t单位s)；% x轴注释Ylabel(振幅x单位厘米)；% y轴注释标题(a的波形图)；%图形标题网格打开；%显示网格子情节(2，1，2)；图(t，y，g)；Xlabel(时间t单位s)；% x轴注释Ylabel(振幅y单位厘米)；% y轴注释标题(b的波形图)；%图形标题网格打开；%显示网格%轴轨迹图图(2)图(x，y，b)网格打开；程序2:计算X和y的波形系数和峰值系数%x y的波形系数和波峰系数clc全部关闭；sf=800。t=0:1/sf:1。x=6 * sin(60 * pi * t)8 * sin(120 * pi * t pi/6)；xav=平均值(绝对值(x)；平均绝对值%xp=最大值(x)；%振幅也是最大值z=x.2；d=平均值(z)；xrms=sqrtm(d)；%有效值Ft=xrms/xav。Fc=xp/xrms。disp(x的波形系数为)；disp(num 2 str(Ft)；Disp(x的波峰因数为)。disp(num 2 str(Fc)；波形系数和波峰系数为%yy=8 * sin(60 * pi * t pi/2)3 * sin(120 * pi * t pi/4)；x=y。xav=平均值(绝对值(x)；平均绝对值%xp=最大值(x)；%振幅也是最大值z=x.2；d=平均值(z)；xrms=sqrtm(d)；%有效值Ft=xrms/xav。Fc=xp/xrms。Disp(y的波形系数为)；disp(num 2 str(Ft)；Disp(y的波峰因数为)。disp(num 2 str(Fc)；程序3:求解振幅谱、相位谱和功率谱%振幅谱、相位谱、功率谱clcfs=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)；%快速傅里叶变换和幅度谱yy=fft(y，N)；% fft变换py1=YY/512；py2=py1 * 2；%单边频谱乘以2mag=ABS(py2)；%寻道幅度f=(0:256)* 800/512；子情节(2，1，1)；图(f，mag(1:257)；%制作频谱图Xlabel(频率(赫兹)；伊拉贝尔(振幅)；标题(b的振幅谱)；网格；%功率谱功率=mag(1:257)。2；子情节(2，1，2)；图(f，功率)；网格打开；标题(b的功率谱)；Xlabel(频率(赫兹)；伊拉贝尔(权力)；%相位频率图pyy=角度(YY)；Pyy=180/pi*pyy。%弧度改变角度图(2)；情节(f，Pyy(1:257)；格里登；标题(b的相位谱)；Xlabel(频率(赫兹)；伊拉贝尔(阶段(度)；程序4:%相关系数和一致性clcclosealldt=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=128a=xcorr(x，y，128)；a=a(N 1:2 * N 1)；图(t(1:N)，a(1:N)；格里登；%找到A和B的相干函数图(2)；N=512。cxy，w=内聚(x，y，N/4，sf，汉明(N/4)，0，均值)；绘图(w，cxy)；网格；标题(连贯性分析)；Xlabel(频率(赫兹)；下面是一个matlab帮助的例子t=0:0.001:0.6x=sin(2 * pi * 50 * t)sin(2 * pi * 120 * t)；y=x 2*randn(尺寸(t)；绘图(1000*t(1:50)，y(1:50)标题(被零均值随机噪声破坏的信号)xlabel(