How to get all positive values in following code

Question

MIch 2022-9-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1809370-how-to-get-all-positive-values-in-following-code

移动：Rik 2022-9-22

Hello.

If possible, I need help with my code.

I wrote the code and I need to get all positive values for a variable named delta but it seems that my code isn't successful enough. Is it possible to save good T values (the ones that result with positive delta) and then when the loop is broken to continue from the place where the loop was broken???

%constants
var_g_T=1.2;
var_d_T=0.05;
C= 0.3;
P= [2.0 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.0 2.5 2.5];
b_r=14;
g_r=[1 2 2 3 4 5 6 6 7 7 8 8 9 10 11 12 12 13 14 14];
r_r=[6 1 7 2 3 4 5 14 5 13 7 9 10 11 12 13 14 8 1 9];
v_gr = size (g_r);
b_gr = v_gr(2);
v_rr = size (r_r);
b_rr = v_rr(2);
kv_g=[3232 5924 5924 3556 3738 2401 6109 6109 5223 5223 6093 6093 2484 3883 3707 5899 5899 2991 5199 5199];  
kv_r=[3233 996 1890 3556 2244 2401 1197 1874 1197 987 1890 1165 2484 2344 3707 987 1874 2991 996 1165];   
%for loop
for i = 1 : 5
    
    delta_neg=false;
    while delta_neg==0
        for m = 1 : 14
            T(i,m) = (var_g_T-var_d_T)*rand(1)+var_d_T;     
        end
        
        for n = 1:b_rr
            gpp(i,n)=P(1,g_r(1,n));
            gvp(i,n)=0.14*T(i,g_r(1,n))/((kv_g(i,n)/gpp(i,n))^0.02 -1);
            rpp(i,n)=P(1,r_r(1,n));
            rvp(i,n)=0.14*T(i,r_r(1,n))/((kv_r(i,n)/rpp(i,n))^0.02 -1);
           
            % is delta positive???
            delta(i,n)=rvp(i,n)-gvp(i,n)-C;
            disp(delta)
            if delta(i,n)>0
                delta_neg=true;
            else
                delta_neg=false;
                break
            end
        end
    end 
end

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Karim 2022-9-21

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1809370-how-to-get-all-positive-values-in-following-code#answer_1057940

在 MATLAB Online 中打开

Im modified the code a bit, see below for the adjustments and some coments in the code.

I stored the T values along with the i and n index into a table. Is this what you needed?

%constants
var_g_T=1.2;
var_d_T=0.05;
C= 0.3;
P= [2.0 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.5 2.0 2.5 2.5];
b_r=14;
g_r=[1 2 2 3 4 5 6 6 7 7 8 8 9 10 11 12 12 13 14 14];
r_r=[6 1 7 2 3 4 5 14 5 13 7 9 10 11 12 13 14 8 1 9];
v_gr = size (g_r);
b_gr = v_gr(2);
v_rr = size (r_r);
b_rr = v_rr(2);
kv_g=[3232 5924 5924 3556 3738 2401 6109 6109 5223 5223 6093 6093 2484 3883 3707 5899 5899 2991 5199 5199];  
kv_r=[3233 996 1890 3556 2244 2401 1197 1874 1197 987 1890 1165 2484 2344 3707 987 1874 2991 996 1165];   
numN = b_rr;
numI = 5;
% compute all values of "T"
T = (var_g_T-var_d_T).*rand(numI,numN) + var_d_T     
T = 5×20
    0.2096    0.8173    1.1229    1.0401    0.3960    1.0598    0.6261    0.7479    0.3235    0.0804    1.1333    1.1929    1.1431    1.1939    1.0399    0.6139    0.1151    1.1760    0.5305    0.3641
    0.9752    0.2826    0.8557    0.5500    0.9210    1.1690    1.0053    1.0433    1.1118    0.4117    0.6160    0.7454    1.0043    0.5711    0.2330    1.1377    0.8652    0.7427    1.1052    0.3710
    0.3143    0.8098    0.8303    0.2676    1.0149    0.5286    0.1634    1.1807    0.1821    0.1114    0.3258    0.2602    0.0637    0.2262    0.9717    0.6669    0.9608    0.5051    0.5073    0.5251
    0.2599    0.6098    0.7954    0.3595    0.2799    0.5029    0.0561    0.3552    0.8500    0.2601    0.1308    0.3171    0.5341    0.7645    1.1256    1.0907    0.6599    0.5355    0.7342    1.1104
    0.3396    0.1736    0.4916    0.9848    0.7131    0.6029    0.8343    0.6765    0.4516    0.5311    0.7465    0.4529    0.6683    0.2592    0.4811    0.8691    0.5817    0.4345    0.0748    0.5035
% intialize some matrices to store the temporary results
[gpp,gvp,rpp,rvp,delta] = deal(zeros(numI,numN));
n = 1:numN;
for i = 1:numI
    
    gpp(i,n) = P(1,g_r(n));
    gvp(i,n) = 0.14*T(i,g_r(n))./((kv_g(n)./gpp(i,n)).^0.02 -1);
    
    rpp(i,n) = P(1,r_r(n));
    rvp(i,n) = 0.14*T(i,r_r(n))./((kv_r(n)./rpp(n)).^0.02 -1);
    
    % evaluate all delta's
    delta(i,n) = rvp(i,n) - gvp(i,n) - C;
end
% extract the values from T that have a positive delta
T_positive = T( delta > 0 )
T_positive = 13×1
    0.2096
    0.2826
    0.4916
    1.0598
    0.1634
    0.3552
    0.7454
    0.3171
    1.1377
    1.0907
% get the i and n value for the positive T values
[i_value,n_value] = find(delta > 0 );
% store the data in a table
MyResult = table(T_positive,i_value,n_value)
MyResult = 13×3 table
    T_positive    i_value    n_value
    __________    _______    _______

     0.20958         1          1   
     0.28256         2          2   
     0.49162         5          3   
      1.0598         1          6   
     0.16338         3          7   
      0.3552         4          8   
     0.74545         2         12   
     0.31709         4         12   
      1.1377         2         16   
      1.0907         4         16   
     0.86914         5         16   
     0.65992         4         17   
     0.50514         3         18   

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

MIch 2022-9-22

移动：Rik 2022-9-22

Thanks for your answer but this is not what I was looking for. I need to get positive delta calculated from the gvp and rvp equations. T is the variable in those equation, and in this example there is 14 T variables (the highest value in g_r and r_r variables)

请先登录，再进行评论。

How to get all positive values in following code

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

回答（1 个）

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

另请参阅

类别

标签

Community Treasure Hunt

How to get all positive values in following code

0 个评论 显示 -2更早的评论隐藏 -2更早的评论

回答（1 个）

1 个评论 显示 -1更早的评论隐藏 -1更早的评论

另请参阅

类别

标签

Community Treasure Hunt

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

1 个评论
显示 -1更早的评论隐藏 -1更早的评论