jis.m

来自「一些数值计算用的MATLAB代码,适合初学MATLAB和数字计算的人使用,根据M」· M 代码 · 共 31 行

function [k,b,a2]=jis(a);
a0=(a/1200);
a2=a0-16;
b=[0:a2];
for i=1:(a2+1)       
    if b(i)<=5
       c(i)=b(i)*0.05;
       d(i)=b(i)*0.05;
    end
    if b(i)>5&b(i)<=20
       c(i)=5*0.05+(b(i)-5)*0.1;
       d(i)=b(i)*0.1;
    end
    if b(i)>20&b(i)<=50
       c(i)=5*0.05+15*0.1+(b(i)-20)*0.15;
       d(i)=b(i)*0.15;
    end
    if b(i)>50&b(i)<=200
       c(i)=5*0.05+15*0.1+30*0.15+(b(i)-50)*0.2;
       d(i)=b(i)*0.2;
    end
end
d=rot90(d,2);
k=c+d;
d=rot90(d,2);