Matlab 复数 绘图我很想知道!根据我下面个程序,我能不能用H1画出它的频谱图?我不是很清楚我的矩阵A求的对吗?这个是公式.最后的H1(w)就是我要根据它来做频谱图close all;clear,clcsyms w;M=[5.31,0,0,0
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/08 21:49:25
![Matlab 复数 绘图我很想知道!根据我下面个程序,我能不能用H1画出它的频谱图?我不是很清楚我的矩阵A求的对吗?这个是公式.最后的H1(w)就是我要根据它来做频谱图close all;clear,clcsyms w;M=[5.31,0,0,0](/uploads/image/z/3763471-31-1.jpg?t=Matlab+%E5%A4%8D%E6%95%B0+%E7%BB%98%E5%9B%BE%E6%88%91%E5%BE%88%E6%83%B3%E7%9F%A5%E9%81%93%21%E6%A0%B9%E6%8D%AE%E6%88%91%E4%B8%8B%E9%9D%A2%E4%B8%AA%E7%A8%8B%E5%BA%8F%2C%E6%88%91%E8%83%BD%E4%B8%8D%E8%83%BD%E7%94%A8H1%E7%94%BB%E5%87%BA%E5%AE%83%E7%9A%84%E9%A2%91%E8%B0%B1%E5%9B%BE%3F%E6%88%91%E4%B8%8D%E6%98%AF%E5%BE%88%E6%B8%85%E6%A5%9A%E6%88%91%E7%9A%84%E7%9F%A9%E9%98%B5A%E6%B1%82%E7%9A%84%E5%AF%B9%E5%90%97%3F%E8%BF%99%E4%B8%AA%E6%98%AF%E5%85%AC%E5%BC%8F.%E6%9C%80%E5%90%8E%E7%9A%84H1%28w%29%E5%B0%B1%E6%98%AF%E6%88%91%E8%A6%81%E6%A0%B9%E6%8D%AE%E5%AE%83%E6%9D%A5%E5%81%9A%E9%A2%91%E8%B0%B1%E5%9B%BEclose+all%3Bclear%2Cclcsyms+w%3BM%3D%5B5.31%2C0%2C0%2C0)
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));