MATLAB: What is wrong in this code?

griddatainterpolationoptimizationproblemwhile loop

Could you please help me to solve the problem? Loop cannot stop.

Attached you can find the file. The function is digdata and the other is code. Thank you.

Best Answer

In this loop:

while abs(res1)>0.01
  As   = vq.*(Us(j).^2)./(16.*(pi.^2).*msksi.*fs1(j));
  res1 = As1(i)-As; 
end

you calculate the same value in each iteration. As will never change its value. Therefore this must be an infinite loop.

Due to the missing comments, the purpose of this loop cannot be guessed. Therefore a reliable suggestion of improvement cannot be given.

Related Solutions

MATLAB: Image Help with loops

What do you mean misread? Is it homework? We're not supposed to just do your homework for you (so you don't get caught cheating), and if it's homework you're supposed to tag it as homework. Please read this http://www.mathworks.com/matlabcentral/answers/8626-how-do-i-get-help-on-homework-questions-on-matlab-answers

So I assume you know how to do a for loop, so give it a try.

for col = 1 : 256
    for row = 1 : 256
       % Code for you to finish....
    end
end

MATLAB: How can put the values inside tolerances

Firstly, you don't need those digdata, digdata2 functions. You're code is very slow because of them. You're loding the same mat files again and again for more than 40000 times. secondly, you should move those lines outside the inner loop, why do the same thing 201 times everytime?!.

So here's the cleaned up code, this is faster so it's much easier now to figure out what's your problem.

ms=1;
msksi=0.013;
Vs1=5.5;
As1=0.05:0.005:1.05;
fs1=0.5:0.005:1.5;
%load the data only once!!
m = load('myarray.mat');
myarray = m.myarray;
d = load('datas2.mat');
datas2 = d.datas2;
loop=0;
for i=1:length(As1)
    %do it here, why repeat the very same thing 200 times everytime!
    vq(i) = griddata(myarray(:,1),myarray(:,2),myarray(:,3),Vs1,As1(i));
    vq2(i) = griddata(datas2(:,1),datas2(:,2),datas2(:,3),Vs1,As1(i));
    for j=1:length(fs1)
        Us(j)=Vs1.*fs1(j);       
        As(i,j) = func_find_As(vq(i),Us(j),mksi,fs1(j));
        fs(i,j) = func_find_fs(vq2(i),Us(j),ms,As1(j));
        [res1(i,j), res2(i,j)] = calculateRes(As(i,j),As1,fs(i,j),fs1);
        [As1, fs1, Us] = checkTolerance(res1,res2);
    end
    loop=loop+1;
end

Best Answer

Related Solutions

MATLAB: Image Help with loops

MATLAB: How can put the values inside tolerances

Related Question