fscan.m

来自「fscan 是在matlab里面的一个频率的扫描程序 可以对于一个特定的对象进行」· M 代码 · 共 34 行

function [Hr,F] = fscan( bus,line,mac_con,RLC )
%FSCAN Summary of this function goes here
%   Detailed explanation goes here
nbus = length(bus(:,1)); 
F=zeros(1,2+nbus);
for x=0.1:0.1:25
    Y=Yfmatrix(bus,line,mac_con,RLC,x);
    %if min(svd(Y))<0.005
     %   if F==[0]
      %      F=[x]
       % else F=[F;x];
        %end
    %end
    [v,D]=eig(Y);
    eigenvalues=zeros(1,nbus);
    for i=1:nbus
        eigenvalues(i)=abs(D(i,i));
    end
    mineigenvalue = min(svd(Y));
    if F==[0] F = [x,1/mineigenvalue,ones(1,rank(Y))./eigenvalues];
    else F = [F;x,1/mineigenvalue,ones(1,rank(Y))./eigenvalues];
    end
end
Hr = [0];
for i=2:249
    if (F(i,2)>F(i-1,2)) & (F(i,2)>F(i+1,2))
        if Hr == [0] Hr = [(i/10),F(i,2)];
        else Hr = [Hr;(i/10),F(i,2)];
        end
    end 
end 
plot(F(:,1),F(:,2));

return