Nonlinear system (has anyone seen this equation form?)
1 次查看(过去 30 天)
显示 更早的评论
Has anyone seen this form of a nonlinear equation with respect to X, but linear with respect to Y & Z? I provided a contour plot within the region for all 3 variables between -2 & 2. The plot is actually Z*conjugate(Z) so that the magnitude is above ZERO. If I am correct I may have seen this before in orbital mechanics, but my undergrad years are a bit behind me. I generated this eq. during a parametric analysis of a 2nd order oscillator. I would like to come up to speed about the physical phenomena this equation might describe.
I would provide a PDF of my results, but I do not see that I can upload a file. So below is my MATLAB code. If you can provide info how to attach or upload a file please respond.
alpha=-2:0.004:2;
rho=10^-6;
for k=1:length(alpha)-1
rho(k+1)=rho(k)*10^(12/(length(alpha)-1));
end
ones(1:k+1)=1;
X=fliplr(-1*ones'*rho);
X=[X ones'*rho];
Y=alpha'*ones;
Y=[Y Y];
z=((1-Y).*((1-Y).^2+1).*((1+Y).^2-1).^2.*X-(1+Y).*((1+Y).^2+1).*((1-Y).^2-1).^2)./(4*(1-Y).*(1+Y).*((1-Y).*(1+(1+Y).^2)-(1+Y).*((1-Y).^2+1).*X));
Z=z.*conj(z);
mesh(X,Y,Z);
[m n]=size(Z);
axis([X(1,1) X(m,n) Y(1,1) Y(m,n) 0 2 0 2])
2 个评论
bym
2012-11-25
ones()
is built in function which you are overwriting in your code, not a good idea
回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Systems of Nonlinear Equations 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!