Matlab 复数 绘图我很想知道!根据我下面个程序,我能不能用H1画出它的频谱图?我不是很清楚我的矩阵A求的对吗?这个是公式.最后的H1(w)就是我要根据它来做频谱图close all;clear,clcsyms w;M=[5.31,0,0,0

来源:学生作业帮助网 编辑:作业帮 时间:2024/11/14 14:47:13

Matlab 复数 绘图我很想知道!根据我下面个程序,我能不能用H1画出它的频谱图?我不是很清楚我的矩阵A求的对吗?这个是公式.最后的H1(w)就是我要根据它来做频谱图close all;clear,clcsyms w;M=[5.31,0,0,0
Matlab 复数 绘图
我很想知道!根据我下面个程序,我能不能用H1画出它的频谱图?
我不是很清楚我的矩阵A求的对吗?




这个是公式.最后的H1(w)就是我要根据它来做频谱图


close all;clear,clc
syms w;
M=[5.31,0,0,0
    0,24.138,0,0
    0,0,10.45,0
    0,0,0,15.602]
K=[309.95,-309.93,0,0
    -309.95,-159.94,-150.01,0
    0,-150.01,261.55,-111.54
    0,0,111.54,811.91
    ]
C=[4000,-4000,0,0
    -4000,8000,-4000,0
    0,-4000,8000,-4000
    0,0,-4000,6698.3
    ]
B=[0,0
    0,0
    0,0
    700.37,2698.3]
A=-w^2*M+j*w*C+K
H=inv(A)*B*[1;j*w];
H1=H(1);


这个是我最终要得到的图形

Matlab 复数 绘图我很想知道!根据我下面个程序,我能不能用H1画出它的频谱图?我不是很清楚我的矩阵A求的对吗?这个是公式.最后的H1(w)就是我要根据它来做频谱图close all;clear,clcsyms w;M=[5.31,0,0,0
可以画图,但画出的图的确和你提供的不一样,估计是算法有问题吧:
clear all;clc;
syms w
M=[5.31,0,0,0;0,24.138,0,0;0,0,10.45,0;0,0,0,15.602];
K=[309.95,-309.93,0,0; -309.95,-159.94,-150.01,0;0,...
-150.01,261.55,-111.54;0,0,111.54,811.91];
C=[4000,-4000,0,0;-4000,8000,-4000,0;0,-4000,8000,...
-4000;0,0,-4000,6698.3];
B=[0,0;0,0;0,0;700.37,2698.3];
A=-w^2*M+j*w*C+K;
H=inv(A)*B*[1;j*w];
H1=H(1);
H1=subs(H1,w,{0:0.2:20});
w1=0:0.2:20;
plot(w1,real(H1));