MATLAB 求助帖

sgj821

clc
clear
L1=0.01;
w=1;
X=zeros(66,2);
R1=zeros(66,1);
L=[1.6;1.6;1.6;1.6;1.6;1.6;1.6;1.6;1.1;1.1;1.1;1.1;1.1;1.1;2;2;2;2;2;2;  1.6;1.6;1.6;1.6;1.6;1.6;1.6;1.6;1.6;1.6;1.6;1.6;  1.6;1.6;1.6;1.6;1.6;1.6;1.6;1.6;1.6;1.6;1.6;1.6;1.6;1.6;1.6;1.6;1.6;1.6;1.6;1.6;1.6;1.6;    1;1;1;1;1;1;1;1;1;1.2;1.2;1.2]
P=[40;40;40;40;80;80;80;80;40;40;40;80;80;80;40;40;40;80;80;80;8;8;8;8;15;15;15;15;40;40;40;40;8;8;8;8;8;8;8;8;8;8;8;8;8;8;8;8;8;5;10;20;30;40;10;20;30;10;20;30;10;20;30;10;20;30]
D=[0.1;0.13;0.2;0.4;0.1;0.13;0.2;0.4;0.1;0.2;0.4;0.1;0.2;0.4;0.1;0.2;0.4;0.1;0.2;0.4;0.1;0.13;0.2;0.4;0.1;0.13;0.2;0.4;0.1;0.1;0.1;0.1;0.008;0.008;0.008;0.008;0.008;0.008;0.008;0.008;0.008;0.009;0.01;0.011;0.012;0.013;0.014;0.015;0.016;0.008;0.008;0.008;0.008;0.008;0.008;0.008;0.008;0.008;0.008;0.008;0.01;0.01;0.01;0.01;0.01;0.01]
L2=[0.9;0.9;0.9;0.9;0.9;0.9;0.9;0.9;0.6;0.6;0.6;0.6;0.6;0.6;1.2;1.2;1.2;1.2;1.2;1.2;0.9;0.9;0.9;0.9;0.9;0.9;0.9;0.9;0.3;0.5;0.7;1.1;0.5;0.6;0.7;0.8;0.9;1;1.1;1.2;1.3;0.9;0.9;0.9;0.9;0.9;0.9;0.9;0.9;0.9;0.9;0.9;0.9;0.9;0.6;0.6;0.6;0.4;0.4;0.4;0.4;0.4;0.4;0.8;0.8;0.8]
R=[402;331;249;200;719;589;439;309;282;187;147;498;323;250;520;309;209;934;549;365;115;97;76;58;188;156;120;88;166;245;324;478;467;563;656;749;842;935;1029;1122;1215;777;721;673;631;591;558;529;504;571;1018;1806;2500;3168;680;1204;1664;455;803;1107;391;692;955;776;1380;1911]
for i=1:66
    %求R1
    m=fix(3*L(i,1)/(pi*D(i,1)));
    a=0;
    for j=1:m
        a=a+(1-tanh(j*pi*D(i,1)/L(i,1)))/(j*j)*sin(j*pi*(L2(i,1)-L1)/2/L(i,1))*sin(j*pi*(L2(i,1)+L1)/2/L(i,1))*sin(j*pi*L2(i,1)/L(i,1));
    end
    a=4/pi/pi*a;
    R1(i,1)=(L2(i,1)-L1)/(w*L(i,1))/((L1+L2(i,1))/L(i,1)*(1-L2(i,1)/L(i,1))-a);

    P(i,1)=log(P(i,1));

    R(i,1)=log(R(i,1))-log(R1(i,1));
end
% X=[P D];
R1
X=[ones(size(P)) P];
b=regress(R,X)

求大神看看这个程序有什么问题  为什么输出结果是虚数？？

lijun2007zhu

我粗略的看了一下，R1计算出来的数组存在负数，导致计算R(i,1)=log(R(i,1))-log(R1(i,1));时，log(R1(i,1))出现虚数，以至于之后的b=regress(R,X)产生虚数。

sgj821

回复 2# lijun2007zhu


   多谢大神帮忙