rdaqi.asv

来自「这是一个基于大气环境评价的PSO程序」· ASV 代码 · 共 25 行

function SR=rdaqi(a)
a=input('put into data please:');
w=[0,0,0];m=[0,0,0];
    R(1)=1/(1+73.34*exp(-0.3391*a(1)/0.02));
    R(2)=1/(1+73.34*exp(-0.3391*a(2)/0.015));
R(3)=1/(1+73.34*exp(-0.3391*a(3)/0.05));
for i=1:3;        
if R(i)>=0&R(i)<=0.0215
    m(i)=0.1+R(i)*0.2/0.0215;w(i)=(m(i)/2)^2;
elseif R(i)>=0.0215&R(i)<=0.0383
    m(i)=0.3+(R(i)-0.0215)*0.2/(0.0383-0.0215);
    w(i)=(m(i)/2)^2;
elseif R(i)>=0.0383&R(i)<=0.1245
    m(i)=0.5+(R(i)-0.0383)*0.2/(0.1245-0.0383);
    w(i)=1-((m(i)-1)/2)^2;
elseif R(i)>=0.1245&R(i)<=0.4506
    m(i)=0.7+(R(i)-0.1245)*0.2/(0.4506-0.1245);
    w(i)=1-((m(i)-1)/2)^2;
elseif R(i)>=0.4506&R(i)<=0.939
    m(i)
    w(i)=3;
end
SR=(w(1)*R(1)+w(2)*R(2)+w(3)*R(3))/8;
disp(SR);
end