Re: [問題]for迴圈和while的程式問題

看板MATLAB作者HAK (no pain no gain)時間12年前 (2013/05/10 14:24)推噓2(2推 0噓 5→)

留言7則, 2人參與討論串3/3 (看更多)

※ 引述《profyang (prof)》之銘言： : ※ 引述《HAK (no pain no gain)》之銘言： : : Matlab新手想請問關於for迴圈與while組合的程式問題 : : 我的問題在於目前有一個4x4的A矩陣 : : X1的初始值為[1;0;0;0],新的X1=A*X1 : : 我想要求n次之後X1中的元素值小於0.00001，並且記錄每一個X1的值 : : 因此寫了下面的程式: : : X1(1:4,1)=[1;0;0;0]; : : for j=1:100 : : while X1(1:4,1)>0.00001 : : n(j,1:4)=X1'; : : X1=A*X1; : : end : : end : : 但是卻進不了迴圈 : : 不知可否請先進點撥一下問題出在哪裡呢 : X1=[1;0;0;0];X2=X1;n=[]; : while sum(abs(X2))>0 : n=[n,X1]; : X1=A*X1; : X2 = X1; : X2(X2<0.00001) = 0; : end : 應該會動吧我用A=0.1試過好像會動不好意思又有問題了假如說我的A矩陣有5個並且已經存成3維矩陣的形式想要對這5個A矩陣做同樣的事情並且把根據不同矩陣所求得的n1存成3維矩陣n2的形式因此寫了這樣的程式: for k=1:5; M=A(:,:,k); X1=[1;0;0;0];X1_temp=X1;n1=[]; while sum(abs(X1_temp))>0 n1=[n1,X1]; X1=M*X1; X1_temp = X1; X1_temp(X1_temp<0.001) = 0; end n2(:,:,k)=n1; end 但是執行的時候出現??? Subscripted assignment dimension mismatch.的訊息可否請問該如何修改呢? -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.116.134.111

推

profyang

05/10 16:57, , 1^F

05/10 16:57, 1^F

→

profyang

05/10 16:57, , 2^F

05/10 16:57, 2^F

→

profyang

05/10 16:58, , 3^F

05/10 16:58, 3^F