I am trying to plot the results of my Matlab programme using a for loop which operates from 1 to 10 in increments of 0.1. However I am getting the following error when I attempt to run my code:
Subscript indices must either be real positive integers or logicals.
Error in BEM (line 101) store_sigma(:,TSR)=sigma;
This is my entire code:
% Inputs R=0.4; % Radius of Rotor B=3; % Number of blades U=1.73; % Fluid velocity Rho=998; % Fluid Density N=9; % Number of Blade Elements Cp_estimate=0.5; % Estimate power coefficient Alpha_design=4; % Design alpha Cl_design=1; % Design lift coefficient angular_velocity=22.92; % Variables i=1; % Counter alpha_new=0; % Initial value for alpha new tolerance=0.00001; % Tolerance Value axial_induction=[0,0,0,0,0,0,0,0,0]; Check=1; % Initial check value axial_induction_old=0; % Initial value for old axial induction factor Cl=[1.3; 1.1; 1; 0.9; 0.86; 0.83; 0.8; 0.75; 0.5]; % Lift Coefficients Cd=[0.027; 0.024; 0.02; 0.019; 0.018; 0.016; 0.013; 0.012; 0.01]; % Drag Coefficients r_local=R/9*(1:9)'; r_over_R=r_local / R; for TSR=1:.1:10 % TSR from 1 to 10 disp(TSR) Check=1; TSR_local=r_over_R .* TSR; Phi=(2/3)*atan(1./TSR_local); C=((8*pi.*r_local) ./ (B.*Cl_design)).*(1-cos(Phi)); sigma=(B*C) ./ (pi.*r_local.*2); axial_induction= 1 ./ (((4.*(sin(Phi).^2)) ./ (sigma.*Cl_design.*cos(Phi)))+1); angular_induction= (1-(3*axial_induction)) ./ ((4.*axial_induction)-1); relative_wind = (U.*(1-axial_induction))./(angular_velocity*R).*(1+angular_induction); F=(2/pi) .* acos(exp(-(((B/2) .* (1-(r_over_R))) ./ ((r_over_R) .* sin(relative_wind))))); % Tip Loss Factor C_T=(sigma .* ((1-axial_induction).^2) .* ((Cl.*cos(relative_wind))+(Cd.*sin(relative_wind)))) ./ ((sin(relative_wind)).^2); while abs(Check)>=tolerance axial_induction_old = axial_induction; TSR_local = TSR .* (r_local./R); % Local Tip Speed Ratio Phi = (2/3) .* atan(1./TSR_local); % Angle of Relative Fluid for i=1:length(C_T) if C_T(i) > 0.96; axial_induction(i) = 1 / (((4*F(i)*cos(relative_wind(i))) / (sigma(i)*Cl(i)))-1); else axial_induction(i) = 1 / (1+(4*F(i)*(sin(relative_wind(i))^2)) / (sigma(i)*Cl(i)*cos(relative_wind(i)))); end; end; D=(8./(TSR.*9)).*(F.*(sin(Phi).^2).*(cos(Phi)-((TSR_local).*(sin(Phi)))).*(sin(Phi)+((TSR_local).*(cos(Phi)))).*(1-(Cd./Cl).*cot(Phi)).*(TSR_local.^2)); Cp=sum(D); Diff=axial_induction-axial_induction_old; Check=max(Diff(:)); Check end store_sigma(:,TSR)=sigma; store_Phi(:,TSR)=Phi; store_TSR_local(:,TSR)=TSR_local; store_axial_induction(:,TSR)=axial_induction; store_angular_induction(:,TSR)=angular_induction; store_relative_wind(:,TSR)=relative_wind; store_Check(:,TSR)=Check; store_Diff(:,TSR)=Diff; store_Cp(:,TSR)=Cp; store_TSR(:,TSR)=TSR; store_F(:,TSR)=F; end Does anybody have any idea of how I could index my stored values so that I can use an interger value for my incremental steps inside the for loop?
Thanks.