0

load-ascii コマンドを使用して (x,y) を読み取るようにこのプログラムを変更するにはどうすればよいですか?

n=0;
sum_x = 0;
sum_y = 0;
sum_x2 = 0;
sum_xy = 0;
disp('This program performs a least-squares fit of an');
disp('input data set to a straight line. Enter the name');
disp('of the file containing the input (x,y) pairs:  ');
filename = input('   ','s');
[fid,msg] = fopen(filename,'rt');
if fid<0
    disp(msg);
else
    [in,count]=fscanf(fid, '%g %g',2);
    while ~feof(fid)
        x=in(1);
        y=in(2);
        n=n+1;
        sum_x=sum_x+x;
        sum_y=sum_y+y;
        sum_x2=sum_x2+x.^2;
        sum_xy=sum_xy+x*y;
        [in,count] = fscanf(fid, '%f',[1 2]);
    end
    fclose(fid);
    x_bar = sum_x / n;
    y_bar = sum_y / n;
    slope = (sum_xy - sum_x*y_bar) / (sum_x2 - sum_x*x_bar);
    y_int = y_bar - slope * x_bar;
    fprintf('Regression coefficients for the least-squares line:\n');
    fprintf('Slope (m)      =%12.3f\n',slope);
    fprintf('Intercept (b)  =%12.3f\n',y_int);
    fprintf('No of points   =%12d\n',n);
end
4

1 に答える 1

0

私があなたを正しければ、おそらくあなたはこれを望んでいます：

x=load('-ascii',filename);
sumx =sum(x(:,1));
sumy =sum(x(:,2));
sumx2=sum(x(:,1).^2);
sumy2=sum(x(:,2).^2);
sumxy=sum(x(:,1).*x(:,2));
n    =size(x,1);

これらも試してみてください

sums =sum(x);
sum2s=sum(x.^2);
sumxy=sum(prod(x,2));
means=mean(x);

という機能があることもご存知ですか？

[slope,int]=regress(x(:,2), [x(:,1) ones(n,1)])

これはあなたの仕事を不要にするかもしれません:)

編集

以下に入力した内容を考えると、データは別々の行にあるのではなく、スペースで区切られた数字のみで構成されているようです。したがって、これが正しい場合は、loadステートメントの後に次の行を挿入する必要があります

x=reshape(x, 2, length(x)/2);

ASCII ファイルを読み取るときは、正確な形式が重要であることに注意してください。数字の数が奇数の場合、これは失敗し、ファイルに改行文字がある場合は誤った値が返されることに注意してください。

あなたが書いた残りのコードはすべて次のようになると思います。

x=load('-ascii',filename);
[slope,int]=regress(x(:,2), [x(:,1) ones(size(x,1),1)])
于 2010-05-01T01:15:31.737 に答える