interpolation to maximum points

Question

Alzapoa 2012-8-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/46656-interpolation-to-maximum-points

Hi all, I have below code to generate a group of graphs. I need to interpolate a line through their maximum points. could you please help? in addition I need to cut the graph below the -0.1 if it is possible. Thanks

Cp=[ ];
landa=(0:0.1:16);
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
    plot(landa,Cp,'b');
    % mesh(landa,B);
    hold on;
    Cp=[ ];
end
xlabel('Tip speed ratio');
ylabel ('Performance corfficient Cp');
hold off;

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

请先登录，再进行评论。

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

Answer 1

Azzi Abdelmalek 2012-8-24

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/46656-interpolation-to-maximum-points#answer_57063

在 MATLAB Online 中打开

 close
Cp=[ ];
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=[ ];
end
xlabel('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)])

2 个评论
显示无隐藏无

Alzapoa 2012-8-24

Hi Azzi, Thanks for the solution. if you look at the result plot the red line goes to landa=0.4382 and then back down which is not meant to do that. it should go through the max values independent from the order of plot. so what to do with that I hope I am clear.