最小二乗法の計算方法について，配列インデックスに関するエラーの解決方法．

正義金田

2023 5 23

0 个回答

1 次查看（30 天）

0 个投票

fx=10^-2:10^-1    % 範囲指定
bx=p(fx)
l1=log10(bx);     
l2=log10(fx);        
MY=mean(l1);          % l1結果の平均
MX=mean(l2);          % l2結果の平均
SXY=0;                % SXY = 0 
for a=1:length(fx)　　　　
    SXY=SXY+(l1(a)-MY)*(l2(a)-MX); 　%最小二乗法による回帰直線の傾き
end
SXY=SXY/length(fx); 
z=SXY/var(l2)    

以上のプログラムでエラーが出ます．（配列インデックスは正の整数または logical 値でなければなりません。）

このエラーを解決するにはどのようにすればよいでしょうか．

3 个评论
显示 1更早的评论隐藏 1更早的评论

正義金田 2023-5-24

在 MATLAB Online 中打开

Hishida様

コメントありがとうございます．

　仰る通り，範囲指定をしている冒頭のfxで（配列インデックスは正の整数または logical 値でなければなりません）というエラーメッセージが出てきます．

　下記コードのように，範囲を１よりも大きな値にすることでz=〇〇といったデータを得られるのですが，１よりも小さい範囲でのzの値を求めたく思っています．（z=〇〇の形すら出力されません）

ご検討のほど，よろしくお願いいたします．

また不備がありましたら，ご連絡いただければと思います．

fx=10^1:10^2    % 範囲指定
bx=p(fx)
l1=log10(bx);     
l2=log10(fx);         
MY=mean(l1);          % l1結果の平均
MX=mean(l2);          % l2結果の平均
SXY=0;                % SXY = 0 
for a=10:20:length(fx)
    SXY=SXY+(l1(a)-MY)*(l2(a)-MX);   %最小二乗法による回帰直線の傾き
end
SXY=SXY/length(fx); 
z=SXY/var(l2)       