MATLAB: Please help, the value of a function is changing with the domain.

functionsplotting

I am required to create a graph of 4 different equations, y1, y2, y3, y4 shown below, where the only difference is this R value. When I start with a domain of theta going from 0-90 degrees (the first graph) you can clearly see a distinction between all 4 equations. When I change the domain from 0-180 degrees, all 4 equations result in identical outputs that overlap. It seems as I get closer to a domain of 180 degrees, the equations become more and more similar. Does anyone know what is causing this problem? Thanks in advance.

if true
R1 = 3;
R2 = 3.5;
R3 = 4;
R4 = 10;
theta = 0:1:180;
y1 = ((pi/2)*sin(theta*pi/180)*(1+ ((cos(theta*pi/180))/(((R1^2) - (sin(theta*pi/180).^2)).^0.5))));
y2 = ((pi/2)*sin(theta*pi/180)*(1+ ((cos(theta*pi/180))/(((R2^2) - (sin(theta*pi/180).^2)).^0.5))));
y3 = ((pi/2)*sin(theta*pi/180)*(1+ ((cos(theta*pi/180))/(((R3^2) - (sin(theta*pi/180).^2)).^0.5))));
y4 = ((pi/2)*sin(theta*pi/180)*(1+ ((cos(theta*pi/180))/(((R4^2) - (sin(theta*pi/180).^2)).^0.5))));
figure; hold all;
plot(theta, y1, theta, y2, theta, y3,theta, y4)

Code:

theta = 0:1:90 :

theta = 0:1:180 :

theta = 0:1:170 :

Best Answer

Hello Joshua, For each variable, try this with extra dots in the locations shown. This is because it appears that you want term-by-term multiplication or division of your vector elements, not some kind of matrix multiplication or division.

y1 = ((pi/2)*sin(theta*pi/180).*(1+ ((cos(theta*pi/180))./(((R1^2) - (sin(theta*pi/180 ...
                              ^                         ^

Interesting effect.

Related Solutions

MATLAB: Solving non linear equation using fsolve

You are not passing theta_n and other parameters to root2d(). One fix is to make root2d a Nested Function. Also, you must choose a different initial point. phi0 = [10,10] worked for me.

function solve_it()
    beta_a = 10;
    alpha_n = 10 ;
    eta= 18;
    i = 18;
    theta_i = 10;
    theta_n = 10 ;
    
    fun = @root2d;
    phi0 = [10,10];
    phi = fsolve(fun,phi0);
    phi_n = phi(1),
    phi_i = phi(2),
    function F = root2d(phi)
    F = zeros(2,1);
    F(1) = (pi*sin((pi*(phi(1) + theta_n))/180)*cos((pi*phi(2))/180)*(cos((pi*theta_n)/180) + tan((pi*i)/180)*tan((pi*theta_i)/180)))/(180*sin((pi*phi(1))/180)*(cos((pi*(phi(1) + theta_n))/180)*cos((phi(2)*pi)/180) + sin((phi(2)*pi)/180)*tan((theta_i*pi)/180))^2) - (pi*cos((pi*phi(1))/180)*(cos((pi*theta_n)/180) + tan((pi*i)/180)*tan((pi*theta_i)/180)))/(180*sin((phi(1)*pi)/180)^2*(cos((pi*(phi(1) + theta_n))/180)*cos((pi*phi(2))/180) + sin((pi*phi(2))/180)*tan((pi*theta_i)/180)));
    F(2) = ((cos((pi*theta_n)/180) + tan((pi*i)/180)*tan((pi*theta_i)/180))*((pi*cos((pi*(phi(1) + theta_n))/180)*sin((pi*phi(2))/180))/180 - (pi*cos((pi*phi(2))/180)*tan((pi*theta_i)/180))/180))/(sin((pi*phi(1))/180)*(cos((pi*(phi(1) + theta_n))/180)*cos((phi(2)*pi)/180) + sin((phi(2)*pi)/180)*tan((theta_i*pi)/180))^2);
    end
end

MATLAB: Got some problems when trying to plot with a shifting angle in 2D Plotting

The posted code moves the start points from (0,0) to (5,5) and the end points of each line are located on the circle around (0,0) with radius 20. I guess that you want a circle around (5,5) with radius 15:

hold on;grid;
ang=30;
theta=0;
for i=0:(360/ang)
  xl = [5, 5 + 15*cos(theta*pi/180)];
  yl = [5, 5 + 15*sin(theta*pi/180)];
  plot(xl,yl,'-r')
  theta = (i+1)*ang;
end

Related Question