MATLAB: Flow control with functions of matrices

if statementmatrix

theta is a matrix [0:0.01:0.5]*pi

I want to evaluate a function t using a flow control of values of theta:

if theta>pi*(1/2 - 0.01)
t=Inf
else
t = tan(theta)
end

How to achieve such a flow control when theta is a matrix? t is assumed to be a matrix with the same dimensions as theta.

Best Answer

t = inf(size(theta));
mask  theta <= pi*(1/2 - 0.01);
t(mask) = tan(theta(mask));

Related Solutions

MATLAB: How to define a function that evaluates in input

function y=T(theta_Degree)
theta=theta_Degree*pi/180;
y=[cos(theta)^2 sin(theta)^2 2*cos(theta)*sin(theta);sin(theta)^2 cos(theta)^2 -2*cos(theta)*sin(theta);-cos(theta)*sin(theta) cos(theta)*sin(theta) cos(theta)^2-sin(theta)^2];
end

Then type in command line

T(40)

and see the results.

MATLAB: I am getting empty sym for vpasolve

Since thetac is an array of 15 elements, you are trying to solve a system of 15 incompatible equations, therefore, there is no answer. You should solve the equation for each value of thetac.

syms theta1 theta2 gama M l
eqn = theta1 == 2*atan((2*tan(theta2/2)*tan(gama/2 - theta2/2) + (4*tan(theta2/2)^2 - 4*M^2*tan(theta2/2)^2 - 4*M^2*tan(gama/2 - theta2/2)^2 - M^2*gama^2*l^2 - 8*M^2*tan(theta2/2)*tan(gama/2 - theta2/2) - 4*M^2*gama*l*tan(theta2/2) - 4*M^2*gama*l*tan(gama/2 - theta2/2) - M^2*gama^2*l^2*tan(theta2/2)^2*tan(gama/2 - theta2/2)^2 + 2*M^2*gama^2*l^2*tan(theta2/2)*tan(gama/2 - theta2/2) + 4*M^2*gama*l*tan(theta2/2)*tan(gama/2 - theta2/2)^2 + 4*M^2*gama*l*tan(theta2/2)^2*tan(gama/2 - theta2/2) + 4)^(1/2) - 2)/(2*tan(theta2/2) + 2*M*tan(gama/2 - theta2/2) + 2*M*tan(theta2/2) + M*gama*l - M*gama*l*tan(theta2/2)*tan(gama/2 - theta2/2)) - (2*tan(theta2/2)*tan(gama/2 - theta2/2))/(2*tan(theta2/2) + 2*M*tan(gama/2 - theta2/2) + 2*M*tan(theta2/2) + M*gama*l - M*gama*l*tan(theta2/2)*tan(gama/2 - theta2/2)));
F=1.06;
gama=pi/F;
thetac=0:0.2:gama;  
l=0.2181818;
M=0.960615;  
for k=1:numel(thetac)
    theta2=thetac(k);
    seqn = subs(eqn);
    sol(k) = vpasolve(seqn)
end
alpha_1=sol./thetac;
plot(thetac,alpha_1)

Best Answer

Related Solutions

MATLAB: How to define a function that evaluates in input

MATLAB: I am getting empty sym for vpasolve

Related Question