MATLAB: 3D surface plots

3d plots

Best Answer

I have made few modifications. In your program, you are not saving the output of P for each diameter and gear ratio.
clear all
close all
clc
%WIND TURBINE SYSTEM PARAMETERS
% D = 16.5:0.1:20.5;    %Rotor diameter in m
%GR= 14.6:0.2:20.6;   %Gear Ratio of the gearbox connecting rotor to generator
beta=0;   % Angle of attack for maximum efficiency of the turbine
%AIR PROPERTIES
rho = 1.27;    % Density of air kg/m^3
Vw = [0.001:0.01:35];
Wgen = 1800;   %Synchronous Speed of the Generator in rpm
i=0;
for D=16.5:0.1:20.5
    i=i+1;
    j=0;
    for GR=14.6:0.2:20.6
        j=j+1;
    Vt= ((Wgen.*pi/30)./GR).*(D/2);
    Lambda =  Vt./Vw;
    %Heier Approach to Calculate the Coefficient of Performance (Cp)
    %Matrix of Coefficient
    C=[0.5 116 0.4 0 0 6 21 0.08 0.035];
    Lambda_i = ((1./(Lambda+(C(8).*beta)))-(C(9)./((beta.*beta.*beta)+1))).^(-1);
    Cp=0.5*((C(2)./Lambda_i)-(C(3).*beta)-C(6)).*exp(-C(7)./Lambda_i);
    %Calculation of Power
    P(i,j,:)=0.5*rho*((pi/4)*(D*D/2).*(Vw.^3)).*(Cp);
    end
end
[xx,yy]=meshgrid(16.5:0.1:20.5,14.6:0.2:20.6);
surf(xx,yy,nanmean(P,3)','FaceAlpha',0.5,'EdgeColor','none')
xlabel('Rotor diameter')
ylabel('Gear Ratio')
colorbar
Related Question