Save array in matrix from for loop

Question

Karl Zammit 2021-12-11

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1608915-save-array-in-matrix-from-for-loop

评论： Karl Zammit 2021-12-12

I am trying to save two arrays (called pstar A and pstarB) in matrix form such that I can than plot each value as a point against the value of x, however, I am being presented with the error "Unable to perform assignment because the left and right sides have a different number of elements."

Any clues on how to do so and how to obtain the plot thereafter?

for x=xo:((xf-xo)/divnum):xf
    % radial coordinate
    rA = x*bA;
    rB = x*bB; 
    % pressure due to rotation and radial variability
    pabsA = (rho*(swirl^2)*(omegaA^2)*((rA^2)-(r0A^2))/2)-((rho*(avolflowA.^2))/(8*(pi^2)*(sA^2)*((rA^2)-(r0A^2))));
    pabsB = (rho*(swirl^2)*(omegaB^2)*((rB^2)-(r0B^2))/2)-((rho*(avolflowB.^2))/(8*(pi^2)*(sB^2)*((rB^2)-(r0B^2))));
    % pstar
    pstarA = pabsA / (rho*(omegaA^2)*(bA^2));
    pstarB = pabsB / (rho*(omegaB^2)*(bB^2));
    psA = [pstarA, psA+1];
    psB = [pstarB, psB+1];
    counter = counter+1;
end

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

请先登录，再进行评论。

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

Answer 1

KSSV 2021-12-11

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1608915-save-array-in-matrix-from-for-loop#answer_852950

编辑：KSSV 2021-12-11

在 MATLAB Online 中打开

X=xo:((xf-xo)/divnum):xf
N = length(X) ; 
pstarA = cell(N,1) ; 
pstarB = cell(N,1) ; 
for i = 1:N
    x = X(i) ; 
    % radial coordinate    
    rA = x*bA;
    rB = x*bB; 
    % pressure due to rotation and radial variability
    pabsA = (rho*(swirl^2)*(omegaA^2)*((rA^2)-(r0A^2))/2)-((rho*(avolflowA.^2))/(8*(pi^2)*(sA^2)*((rA^2)-(r0A^2))));
    pabsB = (rho*(swirl^2)*(omegaB^2)*((rB^2)-(r0B^2))/2)-((rho*(avolflowB.^2))/(8*(pi^2)*(sB^2)*((rB^2)-(r0B^2))));
    % pstar
    pstarA{i} = pabsA / (rho*(omegaA^2)*(bA^2));
    pstarB{i} = pabsB / (rho*(omegaB^2)*(bB^2));
end

Also try:

X=xo:((xf-xo)/divnum):xf
N = length(X) ; 
pstarA = zeros(N,1) ; 
pstarB = zeros(N,1) ; 
for i = 1:N
    x = X(i) ; 
    % radial coordinate    
    rA = x*bA;
    rB = x*bB; 
    % pressure due to rotation and radial variability
    pabsA = (rho*(swirl^2)*(omegaA^2)*((rA^2)-(r0A^2))/2)-((rho*(avolflowA.^2))/(8*(pi^2)*(sA^2)*((rA^2)-(r0A^2))));
    pabsB = (rho*(swirl^2)*(omegaB^2)*((rB^2)-(r0B^2))/2)-((rho*(avolflowB.^2))/(8*(pi^2)*(sB^2)*((rB^2)-(r0B^2))));
    % pstar
    pstarA(i) = pabsA / (rho*(omegaA^2)*(bA^2));
    pstarB(i) = pabsB / (rho*(omegaB^2)*(bB^2));
end