普遍的なコードを作成する

Question

私はtheta関数を計算するための普遍的な方程式の作成に取り組んできました。私が持っている現在の方程式は普遍的ではありません。これは最大2テタまで有効です。私がthetaの計算を2以上にしなければならない場合は、括弧内の数字を手動で変更する必要があります。theta()。

ユニバーサル機能を設定する方法はありますか？ thata(j) = equationのようなものです。私は以下の2つのコードを貼り付けました。一つは動作していますが普遍的ではなく、もう一つは演奏しようとしたが、まだエラーで終わったものです。

コード1：

for iter = 1:num_iters 


    theta(1) = theta(1) - (alpha/m)*((X*theta-y)')*X(:,1), 
    theta(2) = theta(2) - (alpha/m)*((X*theta-y)')*X(:,2); 
    theta_hist(iter,1) = theta(1); 
    theta_hist(iter,2) = theta(2); 

end 

コード2：

num_cols = size(X,2) 
for inter = 1:num_iters 
    for j = 1:num_cols 
     theta(j) = theta(j) - (alpha/m)*((X*theta-y)')*X(:,j), 
    end, 
end, 

このコードの課題は、そのコードがループ二回通過する場合（i = 1、j = 2） 。関数は、thetaの新しい値を使用しています。これは計算を台無しにしています。

ご意見やご提案をいただければ幸いです。この問題の

詳しい説明は以下のリンクに記載されている： https://docs.google.com/document/d/1XwAVV1OBN9BhQ7n60F2oEUoT263t_2hcAAuCaroR-2g/edit?usp=sharing

Answer 1

あなたはあなたのコードをベクトル化すべきです。あなたは、マトリックスdimenstionsおよび/またはサンプルの入力を提供していないが、私はそれを正しく推測している場合、以下はあなたが必要なものを与える必要があります。

% Some random inputs 
alpha = 1; 
m = 1; 
num_iters = 15; 
N = 5; k = 12; 
X = randn(k,N); 
y = randn(k,1); 

theta = zeros(N,1); 
theta_hist = zeros(N,num_iters); 
for i = 1:num_iters 
    theta = theta - (alpha/m)*X'*(X*theta-y); 
    theta_hist(:,i) = theta; 
end 

Answer 2

複数の反復の後、私は自分の質問に答えることができました。すべての助けを https://github.com/Boniface316/Octave/tree/master/GradientDescent_ThetaCalculator

ありがとう：

コードは、以下のリンクにgithubのに掲載されています。