graphing 2 variables in one equation (can not be transposed)
显示 更早的评论
Hi i would like to plot c in terms of k in this equation.... X=Y*sqrt((k^2+(c*w)^2)/(((k-m*w^2)^2+(c*w)^2)))
I know all the values except for C and K so really im trying to do a simple 2D x-y plane graph, but cannot do so easily because I can not (or do not know how to) express k in terms of c or vice versa so that i can use the ezplot command or something similar, Ive been racking my brain over this for a few days now, trying different ways, different commands and they all give me some sort of error.
note that... m=2000 ; Y=0.2 ; w=157.08 ; X=0.1 ; ANY HELP WILL BE GREATLY APPRECIATED thank you
采纳的回答
Walter Roberson
2015-5-27
Square both sides to get rid of the sqrt(). Then multiply both sides by the denominator of the right hand side, to get two polynomials equal to each other. Distribute the Y^2 (that used to be Y but you squared) over the (k^2+(c*w)^2) to get Y^2*k^2 + Y^2*c^2*w^2 . Now you have c^2 times coefficients on each side, so subtract the Y^2*c^2*w^2 from both sides so that you end up with all the c^2 on the left and nothing in c on the right. Now you have of the form something1 * c^2 + something2 = something3 and that is clearly a quadratic so proceed to solve it: re-arrange to something1 * c^2 = something3 - something2, then divide through by the something1 to get c^2 = (something3 - something2) / something1 . Square root both sides, remembering to take both negative and positive square roots, and you have c in terms of +/- sqrt(an expression). Remember to plot both roots.
c1 = sqrt(-((X.^2 - Y.^2) .* (m.^2 .* w.^4 .* X.^2 - 2 .* k .* m .* w .^2 .* X.^2 + k.^2 .* X.^2 - k.^2 .* Y.^2))) ./ (w.*(X.^2 - Y.^2));
c = [c1(:),-c1(:)];
plot(k, c);
2 个评论
Eli Kroitor
2015-5-27
Walter Roberson
2015-5-27
If you have the Symbolic Toolbox then you can solve() the expression. For example,
syms X Y k c w m
solve( X==Y*sqrt((k^2+(c*w)^2)/(((k-m*w^2)^2+(c*w)^2))), c)
更多回答（1 个）
Murali Krishna
2015-5-27
try this.. create a function
function z = myfun(k,c)
m=2000 ; Y=0.2 ; w=157.08 ; X=0.1;
z = -X+Y*sqrt((k^2+(c*w)^2)/(((k-m*w^2)^2+(c*w)^2)));
end
ezplot(@(k,c)myfun(k,c),[kmin,kmax,cmin,cmax]) in mfile..
类别
在 帮助中心 和 File Exchange 中查找有关 Spline Postprocessing 的更多信息
产品
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
