how can i get the nonzero value of Tdim{i}(1),while i have to start value of k=2:12?

1 次查看（过去 30 天）

Deepesh Kumar Gupta 2021-10-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1468231-how-can-i-get-the-nonzero-value-of-tdim-i-1-while-i-have-to-start-value-of-k-2-12

here is the result:

here is the main code:

% calculation of Rs
 theta=0:30*(pi/180):330*(pi/180);
x=ro*cos(theta);
y=ro*sin(theta);
y1=0;
rs=sqrt((x-e).^2+(y-y1).^2);
D=zeros(1000,12);
for i= 1:1000
    for j=1:12
D(i,j)=rs(j)./ri(i);
    end
end
D;
Rs=num2cell(D,2);
% Value of angle calculation
Tdimnot=(Thotspot-Tinf)./(Tref-Tinf);
alpha=atan2((y-y1),(x-e))*180/pi;
for k=1:12
if alpha(k)<0
    alpha(k)=360+alpha(k);
else
    alpha(k)= alpha(k);
end
end
for k=1:12
g(k)=alpha(k)*pi/180;
end
%  calculation of  1000 matrix of (12,12)
N=1000;
C=cell(1,N);
for i=1:N
    A=zeros(12,12);
for k=1:12
    for j=1:12
        if j==1
             A(k,j)=1./(Rs{i}(k))+Bi(i).*log(Rs{i}(k));
        else
        A(k,j)=((j-1)*((Rs{i}(k)).^(j-1)+(Rs{i}(k)).^(1-j))./(Rs{i}(k))+Bi(i)*((Rs{i}(k)).^(j-1)-(Rs{i}(k)).^(1-j)))* cos((j-1).*g(k));
        end
    end
end
C{i}=A;
end
% Calculation of Y
N=1000;
    Y=cell(1,N);
     E=zeros(12,1);
 for i=1:N
E=(C{i})\(B{i});    
Y{i}=E;
 end
 % calculation of each 12 dimensionless temperature at pheripheral
 N=1000;
  Tdim=cell(1,N);
 
 for i=1:N
     F=zeros(12,1);
 for k=2:12
    F(k,:)=1+(Y{i}(1)).*log(Rs{i}(k))+(Y{i}(k))*((Rs{i}(k)).^(k-1)-(Rs{i}(k)).^(k-1)).*cos((k-1).*g(k));
 end
 Tdim{i}=F;
 