MATLAB: My equation gives only one value for 1×200 vector variables

plotting

%Initial Values and Constants n_i = 1; n_t = 1.33; theta_i = linspace(0,pi/2,200); theta_t = asin(n_i/n_t*sin(theta_i));

%Calculate results R_p = ((n_t*cos(theta_i)-cos(theta_t))/(cos(theta_t)+n_t*cos(theta_i)))^2;

My equation R_p give one single value of 0.0029. I want to plot this function for theta_i, but I can't. What should I do?

Best Answer

You forgot the . symbol when doing elementwise operations. Try this instead:

R_p = ((n_t*cos(theta_i)-cos(theta_t))./(cos(theta_t)+n_t*cos(theta_i))).^2;

Related Solutions

MATLAB: I have a file 1×200 struct how can I distribute it to the worker (32 worker)

Just create the struct normally before the parfor() loop, and inside the parfor() loop, index the struct by the parfor variable. MATLAB will automatically distribute only the needed portion to each worker.

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

Best Answer

Related Solutions

MATLAB: I have a file 1×200 struct how can I distribute it to the worker (32 worker)

MATLAB: Solving non linear equation using fsolve

Related Question