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));