程序:
%霍尔茨法
w=0:0.0001:2; %设置固有频率范围以及精度0.0001
x=zeros(100,20001);
x(1,:)=[ones(1,20001)]; %令振幅X1=1
sum=[zeros(1,20001)]
for i=1:99 %求解各阶振幅
sum=sum+x(i,:)
x(i+1,:)=x(i,:)-sum.*w.^2+x(1,:);
end
sums=[zeros(1,20001)];
for i=1:100
程序: %霍尔茨法 w=0:0.0001:2; %设置固有频率范围以及精度0.0001 x=zeros(100,20001); x(1,:)=[ones(1,20001)]; % 令 振幅X1=1 sum=[zeros(1,20001)] for i=1:99 %求解各阶振幅 sum=sum+x(i,:) x(i+1,:)=x(i,:)-sum.*w.^2+x(1,:); end
程序:
%霍尔茨法
w=0:0.0001:2; %设置固有频率范围以及精度0.0001
x=zeros(100,20001);
x(1,:)=[ones(1,20001)]; %令振幅X1=1
sum=[zeros(1,20001)]
for i=1:99 %求解各阶振幅
sum=sum+x(i,:)
x(i+1,:)=x(i,:)-sum.*w.^2+x(1,:);
end
sums=[zeros(1,20001)];
for i=1:100