cp matrix dimensions
B=2;landa=0:2:20;Cp=0.5*(116/landa+0.4*B)*exp(-21/landa);plot(lnda,Cp) closeCp=[ ];landa=(0:0.1:16);result=[];for B=0:3:30 for i=1:1:length(landa) landa_i=1/((1/(landa(i)+0.08*B))-(0.035/(B^3+1))); Cpi=0.22*((116/landa_i)-(0.4*B)-5)*exp(-12.5/landa_i); Cp=[Cp Cpi]; end if B==0 Cpmax=max(Cp); end result=[result Cp']; plot(landa,Cp,'b'); % mesh(landa,B);
hold on; Cp=[ ];endxlabel('Tip speed ratio');ylabel ('Performance corfficient Cp');%hold off;
[i1,j1]=max(result)xi=landa(j1);yi=i1;xyi=sortrows([xi' yi'],1)hold on x1=landa(min(j1):max(j1)) y1=interp1(xyi(:,1),xyi(:,2),x1); plot(x1,y1,'r') set(gca,'ylim',[-0.1 max(i1)])clc,clearc1 = 1;c2 = 2;c3 = 3;L = 5;syms landaf = (c1-c2*landa^4)*((sinh(landa*L))^2 - (sin(landa*L))^2) - (c1-c2*landa^4)*(cosh(landa*L)+cos(landa*L))^2 + c3*(cosh(landa*L)-cos(landa*L))*(sinh(landa*L)+sin(landa*L)) - (cosh(landa*L)-cos(landa*L))*(sinh(landa*L)+sin(landa*L));F = matlabFunction(f);x = linspace(-1,1)/3;y = F(x);[~,ix] = min(abs(y));plot(x,y,x(ix),y(ix),'*r')grid on
Best Answer