how to display those value of M, which is only less than 1?

Question

SAHIL SAHOO 2022-12-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1869992-how-to-display-those-value-of-m-which-is-only-less-than-1

评论： Zahrah Walid 2022-12-3

ti = 0;
tf = 70E-8;
tspan=[ti tf];
KC = 1E-4;
for j = 1:100
y0= [(10e-2)*rand(2,1); ((-3.14).*rand(1,1) + (3.14).*rand(1,1));
      (10e-2)*rand(2,1); ((-3.14).*rand(1,1) + (3.14).*rand(1,1));
      (10e-2)*rand(2,1); ((-3.14).*rand(1,1) + (3.14).*rand(1,1));
      (10e-2)*rand(2,1); ((-3.14).*rand(1,1) + (3.14).*rand(1,1));
      (10e-2)*rand(2,1); ((-3.14).*rand(1,1) + (3.14).*rand(1,1));
      (10e-2)*rand(2,1); ((-3.14).*rand(1,1) + (3.14).*rand(1,1));
      (10e-2)*rand(2,1); ((-3.14).*rand(1,1) + (3.14).*rand(1,1));
      (10e-2)*rand(2,1); ((-3.14).*rand(1,1) + (3.14).*rand(1,1));
      (10e-2)*rand(2,1); ((-3.14).*rand(1,1) + (3.14).*rand(1,1));
      (10e-2)*rand(2,1); ((-3.14).*rand(1,1) + (3.14).*rand(1,1));
      (10e-2)*rand(2,1); ((-3.14).*rand(1,1) + (3.14).*rand(1,1));
      (10e-2)*rand(2,1); ((-3.14).*rand(1,1) + (3.14).*rand(1,1));
      (10e-2)*rand(2,1); ((-3.14).*rand(1,1) + (3.14).*rand(1,1));
      (10e-2)*rand(2,1); ((-3.14).*rand(1,1) + (3.14).*rand(1,1));
      (10e-2)*rand(2,1); ((-3.14).*rand(1,1) + (3.14).*rand(1,1));
      (10e-2)*rand(2,1); ((-3.14).*rand(1,1) + (3.14).*rand(1,1));
      (10e-2)*rand(2,1); ((-3.14).*rand(1,1) + (3.14).*rand(1,1));
      (10e-2)*rand(2,1); ((-3.14).*rand(1,1) + (3.14).*rand(1,1));
      (10e-2)*rand(2,1); ((-3.14).*rand(1,1) + (3.14).*rand(1,1));
      (10e-2)*rand(2,1); ((-3.14).*rand(1,1) + (3.14).*rand(1,1));
      ((-3.14).*rand(20,1) + (3.14).*rand(20,1))];
yita_mn = [
    0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1;
    1 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
    0 1 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
    0 0 1 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
    0 0 0 1 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
    0 0 0 0 1 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0;
    0 0 0 0 0 1 0 1 0 0 0 0 0 0 0 0 0 0 0 0;
    0 0 0 0 0 0 1 0 1 0 0 0 0 0 0 0 0 0 0 0;
    0 0 0 0 0 0 0 1 0 1 0 0 0 0 0 0 0 0 0 0;
    0 0 0 0 0 0 0 0 1 0 1 0 0 0 0 0 0 0 0 0;
    0 0 0 0 0 0 0 0 0 1 0 1 0 0 0 0 0 0 0 0;
    0 0 0 0 0 0 0 0 0 0 1 0 1 0 0 0 0 0 0 0;
    0 0 0 0 0 0 0 0 0 0 0 1 0 1 0 0 0 0 0 0;
    0 0 0 0 0 0 0 0 0 0 0 0 1 0 1 0 0 0 0 0;
    0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 1 0 0 0 0;
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 1 0 0 0;
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 1 0 0;
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 1 0;
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 1;
    1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0
    ]*(KC);
N = 20;
tp = 1E-12;
o = sort(10e2*rand(1,20),'ascend');
[T,Y]= ode45(@(t,y) rate_eq(t,y,yita_mn,N,o),tspan./tp,y0);
r = ((1/20).*( exp(i.*Y(:,3)) + exp(i.*Y(:,6)) + exp(i.*Y(:,9)) + exp(i.*Y(:,12)) + exp(i.*Y(:,15)) ...
    +exp(i.*Y(:,18)) +exp(i.*Y(:,21)) +exp(i.*Y(:,24)) + exp(i.*Y(:,27)) + exp(i.*Y(:,30)) + exp(i.*Y(:,33)) ...
    + exp(i.*Y(:,36)) + exp(i.*Y(:,39)) +exp(i.*Y(:,42)) + exp(i.*Y(:,45)) + exp(i.*Y(:,48)) + exp(i.*Y(:,51)) + exp(i.*Y(:,54))+ exp(i.*Y(:,57)) + exp(i.*Y(:,60))));
M(j) = max(abs(r))
end
% here the 100 values of M will be generate, I want to know how many out of
% 100 have the M value less than 1.000?
function dy = rate_eq(t,y,yita_mn,N,o)
dy = zeros(4*N,1);
dGdt = zeros(N,1);
dAdt = zeros(N,1);
dOdt = zeros(N,1);
P = 2;
a = 1;
T = 2000;
tp = 1E-12;
Gt = y(1:3:3*N-2);
At = y(2:3:3*N-1);
Ot = y(3:3:3*N-0);
k = 1E-4;
for i = 1:N
    dGdt(i) = (P - Gt(i) - (1 + 2.*Gt(i)).*((At(i)))^2)./T ;
    dAdt(i) = Gt(i)*(At(i));
    dOdt(i) = -a.*Gt(i) + o(1,i).*tp;  
    for j = 1:N
        dAdt(i) = dAdt(i) + yita_mn(i,j)*(At(j))*cos(Ot(j)-Ot(i));
        dOdt(i) = dOdt(i) + yita_mn(i,j)*((At(j)/At(i)))*sin(Ot(j)-Ot(i));
    end
end
dy(1:3:3*N-2) = dGdt;
dy(2:3:3*N-1) = dAdt;
dy(3:3:3*N-0) = dOdt;
n1 = (1:20)';
n2 = circshift(n1,-1);
n61 = n1 +60;
n62 = circshift(n61,-1);
n80 = circshift(n61,1);
j2 =  3*(1:20)-1;
j5 = circshift(j2,-1);
j8 = circshift(j2,-2);
j59 = circshift(j2,1);
dy(n61) = (o(1,n2).' - o(1,n1).').*tp + a.*(Gt(n2) - Gt(n1)) - (k).*(y(j2)./y(j5)).*sin(y(n61)) - (k).*(y( j5)./y(j2)).*sin(y(n61)) + (k).*(y(j8)./y(j5)).*sin(y(n62)) + (k).*(y(j59)./y(j2)).*sin(y(n80));  
end

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

请先登录，再进行评论。

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

Answer 1

Zahrah Walid 2022-12-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1869992-how-to-display-those-value-of-m-which-is-only-less-than-1#answer_1119002

在 MATLAB Online 中打开

I'm not sure if you want to replace the values higher than 1 by one specific value or you want to remove them completely; if it is the latter you can add this statment after the for loop:

M=M(int8(M)<1)

Take care that due to Mtlab precision the "0.9999999" can be displayed to you as "1.000" but when you do any logical opertation it deal by its actual value (because this case exists in your code). If you want to print actual values less than 1, you should use this statment instead:

M=M(M<1)

If your purpose was the former, then you can add if condition inside your for loop as follows:

if (max(abs(r))<1)
    M(j)=max(abs(r))
else
    M(j)=1; %the value you want
end

Care to not use this if condition if you want to remove the values higher than 1 as it will results in NAN elements in your array.

2 个评论
显示无隐藏无

SAHIL SAHOO 2022-12-3

the M value are like this 1.000 0.6453 0.2345 1.000 1.000 1.000 0.2353.............. I want to know only those term which are less than 1.000 for an example, from this set of array i wanted the 0.6453 0.2345 0.2353 only these.

in my program there are 100 value of M I wanted to know how many of them is less than 1.000.

Zahrah Walid 2022-12-3