nihe.m

来自「该方法用来计算变压器波过程」· M 代码 · 共 20 行

format long
k=0;
t=0;
tf=1.2;
tb=50;
for p=150:0.1:200
    for p1=0.01:0.001:0.02
        for t1=0.1:0.01:0.4
            i=i+1
            if ((abs(0.3*(p^(1/(1-p)))-p^(p/(1-p)))-exp(-p1*t1)+exp(-p1*t1*p))<0.1)&(abs(0.9*(p^(1/(1-p))-p^(p/(1-p)))-exp(-p1*(t1+0.6*tf))+exp(-p1*(t1+0.6*tf)*p))<0.01)&(abs(exp(-p1*tb)-exp(-p*p1*tb)-0.5*(p^(1/(1-p))-p^(p/(1-p))))<0.01)
                k=k+1
                a(k)=p
                b(k)=p1
                c(k)=t1
                disp(t1);disp(p);disp(p1);
            end
        end
    end
end