use of If else condition

Question

Akshay Pratap Singh 2019-2-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/446568-use-of-if-else-condition

关闭： MATLAB Answer Bot 2021-8-20

在 MATLAB Online 中打开

Hi all, I am trying to plot a function in which If else condition is required.

I wrote a following code, May any body correct it...

code:

clear all
clc
format longEng
syms  z1 z2 z
phi=(pi/180)*39;
delta=(pi/180)*26;
gma=18.4;
h=4;
h1=1.91;
h2=0.08;
L=h+h1+h2;
beta=1;
alfa=1;
z1=z-h;
z2=z-h-h1;
Ra1=-1; 
Ra2=-(alfa*(z2/h2))^0.5;
Rp1=3*(beta*(1-(z1/(h1))))^0.5;
Rp2=3*(alfa*(z2/h2))^0.5;
delma1=0.5*(1-Ra1)*delta;
delma2=-0.5*(1-Ra2)*delta;
delmp1=0.5*(Rp1-1)*delta;
delmp2=0.5*(Rp2-1)*delta;
ka1=(2*cos(phi)^2)/(cos(phi)^2*(1+Ra1)+cos(delma1)*(1-Ra1)*(1+sqrt((sin(phi+delma1)*sin(phi))/cos(delma1)))^2);
ka2=(2*cos(phi)^2)/(cos(phi)^2*(1+Ra2)+cos(delma2)*(1-Ra2)*(1+sqrt((sin(phi+delma2)*sin(phi))/cos(delma2)))^2);
kp11=(2*cos(phi)^2)/(cos(phi)^2*(1+Rp1)+cos(delmp1)*(1-Rp1)*(1+sqrt((sin(phi+delmp1)*sin(phi))/cos(delmp1)))^2);
kp12=1+0.5*(Rp1-1)*((cos(phi)^2/(cos(delmp1)*(-sqrt((sin(phi+delmp1)*sin(phi))/cos(delmp1))+1)^2))-1);
if Rp1<1 
    k=kp11;
else
    k=kp12;
end
kp21=(2*cos(phi)^2)/(cos(phi)^2*(1+Rp2)+cos(delmp2)*(1-Rp2)*(1+sqrt((sin(phi+delmp2)*sin(phi))/cos(delmp2)))^2);
kp22=1+0.5*(Rp2-1)*((cos(phi)^2/(cos(delmp2)*(-sqrt((sin(phi+delmp2)*sin(phi))/cos(delmp2))+1)^2))-1);
if Rp2<1 
    k=kp21;
else
    k=kp22;
end
 
pa1=ka1; 
pp1=k;
pp2=k;
pa2=-ka2;
fplot(pa1,[0 h+h1])
hold on
fplot(pp1,[h h+h1])
hold on
fplot(pp2,[h+h1 L])
hold on
fplot(pa2,[h+h1 L])
hold off
grid on

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Geoff Hayes 2019-2-23

Akshay - are Rp1 and Rp2 symbolic expressions? Do you need to convert it to a double (numeric value) before trying to compare either with 1?

此问题已关闭。

Answer 1

Naman Bhaia 2019-2-26

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/446568-use-of-if-else-condition#answer_362729

Hey Akshay,

I have looked at the question and I believe it is similar to a previously answered question, I hope this link helps.

3 个评论
显示 1更早的评论隐藏 1更早的评论

Naman Bhaia 2019-2-26

So I havent looked at what the code is trying to do but to resolve the issue, define the value for 'z' somewhere and then do one of the following for comparing the values: