Results as a loop

Question

Mustafa Vural 2020-6-20

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/551815-results-as-a-loop

编辑： KALYAN ACHARJYA 2020-6-21

I want to make my results as a loop.

When I change the parameter b_A or T_A from 1:5 to 1:3, I want that my results also change. But with this code, I need to change manually the result code, every time I change the loop of the parameter. Can someone help me plz?

This code is to estimate a 3 parameter weilbul distribution.

clear all
n   = 1000;
t0  = 0.5;
b_A = 1:5;
T_A = 1:5;
LowerBound= [0 0 0];
rng('shuffle');
data    = zeros(n,length(b_A),length(T_A));
params  = zeros(length(b_A),length(T_A));
data2P = zeros(n,length(b_A),length(T_A));
params2p = zeros(length(b_A),4,length(T_A));
for k= T_A
    for i= b_A
        data(:,i,k) = wblrnd(i,k, [n,1]) + t0;
        data2P(:,i,k) = wblrnd(b_A(i),T_A(k), [n,1]);
        start       = [i k t0];
        custompdf   = @(x,a,b,c) (x>c).*(b/a).*(((x-c)/a).^(b-1)).*exp(-((x-c)/a).^b);
        opt         = statset('MaxIter',1e5,'MaxFunEvals',1e5,'FunValCheck','off');
        params(i,1:3,k) = mle(data(:,i,k),'pdf',custompdf,'start',start,'Options',opt,'LowerBound',LowerBound,'UpperBound',[Inf Inf min(data(:,i,k))])
        params(i,4,k)   = i;
        params(i,5,k)   = k;
        params(i,6,k)   = t0;
        
        params2p(i,1:2,k) = wblfit(data2P(:,i,k));
        params2p(i,3,k)   = i;
        params2p(i,4,k)   = k;
    end
end
result(1:5,:)= params(:,:,1);
result(6:10,:)= params(:,:,2);
result(11:15,:)= params(:,:,3);
result(16:20,:)= params(:,:,4);
result(21:25,:)= params(:,:,5);

2 个评论
显示无隐藏无

Image Analyst 2020-6-20

If you don't want to change them by changing your code, how do you want to change them? Where would you change these values if not in your code? Do you want to make a list of the starting and stopping values in a text file or somewhere, that your code read in? Do you want your code to ask the user? I have no idea - you forgot to say.

Mustafa Vural 2020-6-21

在 MATLAB Online 中打开

Oh I am sorry,

at the end you see:

result(1:5,:)= params(:,:,1);
result(6:10,:)= params(:,:,2);
result(11:15,:)= params(:,:,3);
result(16:20,:)= params(:,:,4);
result(21:25,:)= params(:,:,5);

this is specially for the parameters b=1:5 and T=1:5. You see the lines:

result(1:5,:)

result(6:10,:) and so on.

But when I change for example T to 1:3, I need to change the code like:

result(1:3,:)

(result(4:6,:) and so on.

Now I want to create a code, that it changes automatically when I change my parameter b and T.

Have you any idea?

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

KALYAN ACHARJYA 2020-6-21

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/551815-results-as-a-loop#answer_454891

编辑：KALYAN ACHARJYA 2020-6-21

在 MATLAB Online 中打开

%Last section

data1=0;
results=[];
for i=1:length(T)
result(b+data1,:)= params(:,:,T(i));
data1=data1+5;
end 

2 个评论
显示无隐藏无

madhan ravi 2020-6-21

Preallocation?

KALYAN ACHARJYA 2020-6-21

编辑：KALYAN ACHARJYA 2020-6-21

Thanks Madhan, sorry again, I missed the same. But, I have no idea about exact size of results, its depends on b, and author said it can be change any sizes.

请先登录，再进行评论。

Results as a loop

2 个评论
显示无隐藏无

回答（1 个）

2 个评论
显示无隐藏无

另请参阅

类别

标签

Community Treasure Hunt

Results as a loop

2 个评论 显示 无隐藏 无

回答（1 个）

2 个评论 显示 无隐藏 无

另请参阅

类别

标签

Community Treasure Hunt

2 个评论
显示无隐藏无

2 个评论
显示无隐藏无