I have an expression I need to plot. It depends on k_x, k_y, and k_z. I want to plot it with the expression (an eigenvalue) on the y-axis and the vector k=(k_x, k_y, k_z) on the x-axis, where the variables along the x-axis go: (pi, 0, pi) –> (0, 0, 0) –> (0, 0, pi) –> (pi, 0, pi). I realise k_y is always zero so the expression could be simplified, but it still leaves k_x and k_z changing independently of each other. I know this can be plotted because I'm essentially trying to replicate a figure in a paper I've seen. It's figure 2 in this: https://arxiv.org/pdf/1303.5784.pdf – I hope this makes sense, Thank you
MATLAB: How do you make a 2D plot where the x-axis represents more than one variable
plotting
Related Solutions
Hello,
I found your problem. It is the "axis equal" commands in the code for the figure 4 plots (3 instances). I donwloaded your code and ran it without the "axis equal" commands in the figure 4 plot code and got the following result:
If you look at the range of data on the X axis versus the range of data on the Y axis, you can see why. When you use axis equal, it was choosing the scale from the X axis and applying it to the Y axis, which means that the Y axis would be visually compressed.
Another thing that the axis equal command does is sets the axis mode settings to manual. Doing that freezes all axis limits at their current values, which is why your xlim and ylim commands had no effect on your plots. See the help page on "axis" for more info:
I figured out your problem by looking at the follwoing sections from the parameter descriptions, near the bottom of the "axis" help page:
Hope this helps =)
"...can't seem to figure out how to pass in the known variables of A to this section of code:"
The MATLAB documentation explains how:
Simply define your function with all of the required inputs:
function F = myfun(k0,kx,ky,kz,eE,eH,eO)M = [k0^2*eO-ky^2-kz^2, k0^2*eH+kx*ky, kx*kz; kx*ky-k0^2*eH, k0^2*eO-kx^2-kz^2, ky*kz; kx*kz, ky*kz, k0^2*eE-kx^2-ky^2];F = det(M);end>> k0 = 1.1;>> kx = 1.2;>> ky = 1.3;>> eE = 1.4;>> eH = 1.5;>> eO = 1.6;>> fun = @(kz)myfun(k0,kx,ky,kz,eE,eH,eO);>> kz = fzero(fun,pi)kz = 0.51807Check the function value is zero (within floating point tolerances):
>> fun(kz)ans = 1.4664e-16
Best Answer