rdaqi.m

这是一个基于大气环境评价的PSO程序

function SR=rdaqi(a)
a=input('put into data please:');
w=[0,0,0];m=[0,0,0];
    R(1)=1/(1+73.3388*exp(-0.3391*a(1)/0.02));
    R(2)=1/(1+73.3388*exp(-0.3391*a(2)/0.015));
R(3)=1/(1+73.3388*exp(-0.3391*a(3)/0.05));
for i=1:3;        
if R(i)>=0&R(i)<=0.0188
    m(i)=0.1+R(i)*0.2/0.0215;w(i)=(m(i)/2)^2;
elseif R(i)>=0.0188&R(i)<=0.0335
    m(i)=0.3+(R(i)-0.0188)*0.2/(0.0335-0.0188);
    w(i)=(m(i)/2)^2;
elseif R(i)>=0.0335&R(i)<=0.1098
    m(i)=0.5+(R(i)-0.0335)*0.2/(0.1098-0.0335);
    w(i)=1-((m(i)-1)/2)^2;
elseif R(i)>=0.1098&R(i)<=0.3723
    m(i)=0.7+(R(i)-0.1098)*0.2/(0.3723-0.1098);
    w(i)=1-((m(i)-1)/2)^2;
elseif R(i)>=0.3723&R(i)<=0.9308
    m(i)=0.9+(R(i)-0.3723)*0.1/(0.9308-0.3723);
    w(i)=1-((m(i)-1)/2)^2;
end
SR=(w(1)*R(1)+w(2)*R(2)+w(3)*R(3))/sum(w);
disp(SR);
end