matlab 隐函数变成极坐标

qianxunzn

怎样将在直角坐标下画出来的隐函数图形转化成极坐标的，求帮助

guocong89

获得数据点，计算极坐标值再计算

只是这样好像只有半个圆

1. fig=ezplot('x^2+y^2-1')
   
2. x=get(fig,'xdata');
   
3. y=get(fig,'ydata');
   
4. r=sqrt(x.^2+y.^2);
   
5. theta=acos(y./r)；polar(theta,r)

复制代码

qianxunzn

fig=ezplot('x^2+y^2-1')
x=get(fig,'xdata');
y=get(fig,'ydata');
r=sqrt(x.^2+y.^2);
theta=acos(y./r)；polar(theta,r)
fig =
    1.520021972656250e+002
??? Error using ==> plus
Matrix dimensions must agree.

楼主你好 你的程序运行之后是这样
我吧我的程序给你看看  这个隐函数是在直角坐标下的图形，求转换为极坐标
clear all;close all
ezplot('100*((3/y^2-(1-2/y^2+3/y^4)*cos(2*x))^2+((1+2/y^2-3/y^4)*sin(2*x))^2)-25*(15+10/y^2*cos(2*x)-5*sqrt((3/y^2-(1-2/y^2+3/y^4)*cos(2*x))^2+((1+2/y^2-3/y^4)*sin(2*x))^2))-2.5',[0,2*pi]);

guocong89

1. clear all;close all
   
2. fig=ezplot('100*((3/y^2-(1-2/y^2+3/y^4)*cos(2*x))^2+((1+2/y^2-3/y^4)*sin(2*x))^2)-25*(15+10/y^2*cos(2*x)-5*sqrt((3/y^2-(1-2/y^2+3/y^4)*cos(2*x))^2+((1+2/y^2-3/y^4)*sin(2*x))^2))-2.5',[0,2*pi]);
   
3. x=get(fig,'xdata');
   
4. y=get(fig,'ydata');
   
5. r=sqrt(x.^2+y.^2);
   
6. theta=acos(y./r);
   
7. figure;
   
8. polar(theta,r)

复制代码

qianxunzn

楼主，第一个图对着呢，那个图的x坐标就是度数，总共360，每一度对应一个数字，这样画出来应该是一个闭合的曲线，第二个图好像不对

guocong89

楼主没把问题说清楚，把 polar(theta,r) 改成 polar(x,y) 即可

qianxunzn

clear all;close all
fig=ezplot('100*((3/y^2-(1-2/y^2+3/y^4)*cos(2*x))^2+((1+2/y^2-3/y^4)*sin(2*x))^2)-25*(15+10/y^2*cos(2*x)-5*sqrt((3/y^2-(1-2/y^2+3/y^4)*cos(2*x))^2+((1+2/y^2-3/y^4)*sin(2*x))^2))-2.5',[0,2*pi]);
x=get(fig,'xdata');
y=get(fig,'ydata');
figure;
polar(x,y)
出现??? Error using ==> polar
THETA and RHO must be the same size.
还是不行啊

qianxunzn

楼主是个大好人，给我解决了一个大问题，帮人帮到底，非常的感谢楼主

qianxunzn

楼主，有个新问题，就是下面这张图，一个x对应很多y值，我只想取大于1的那些y值，其他的都和昨天的一样，这个要怎么定义啊这个是新的程序
clear all;close all
ezplot('100*((3/(2*y^2)+1/2*(1-2/y^2+3/y^4)*cos(2*x))^2+1/4*((1+2/y^2-3/y^4)*sin(2*x))^2)-25*(7.5+5/y^2*cos(2*x)-5*sqrt((3/(2*y^2)+1/2*(1-2/y^2+3/y^4)*cos(2*x))^2+1/4*((1+2/y^2-3/y^4)*sin(2*x))^2))-2.5',[0,2*pi])
x=get(fig,'xdata');
y=get(fig,'ydata');
r=sqrt(x.^2+y.^2);
theta=acos(y./r);
figure;
polar(x,y)

qianxunzn

图片弄不上去了