《量化投资:以MATLAB为工具》

MATLAB技术论坛

 找回密码
 注册账号
查看: 414|回复: 0
收起左侧

[提问] 隐函数画图,想把数据调出来

[复制链接]
发表于 2017-10-19 20:51:58 | 显示全部楼层 |阅读模式
10贝壳
clear
clc
syms p1 p2 pr pa m r1 r2 p h hc hd w x v ld k c ra b o O k H;

m=1;
h=0;
hc=0;
hd=0;
v=0.5;
c=1;
b=0;
ld=0;
O=(1+ld)*w^2/pi^4;
x=pi*c/2;
k=(v+hd)/(1+hc);
o=(2-(v+hd)/(1+hc))*m^2-(p*(1-b)-h)/(1+hc);
r1=sqrt(m^2-(p-h)/(2+2*hc)+sqrt(((m^2)-(p-h)/(2+2*hc))^2-(m^4-O/(1+hc))));
r2=sqrt(-(m^2-(p-h)/(2+2*hc))+sqrt(((m^2)-(p-h)/(2+2*hc))^2-(m^4-O/(1+hc))));
H=[(r1^2-k*m^2)*cosh(r1*x) -(r2^2+k*m^2)*cos(r2*x);r1*(r1^2-o)*sinh(r1*x) r2*(r2^2+o)*sin(r2*x)];
eqE=(det(H));
subplot(211)
ezplot(eqE,[0 6 0 35]);
h=findobj(gca,'type','patch');%把ezplot隐函数绘图中的数据以矩阵形式调用出来再用plot画图
x=get(h,'xdata');
y=get(h,'ydata');
plot(x,y,'r')
subplot(212)
axis([0 6 0 35])

view(90,-90)
hold on
用ezplot可以把图出来,但是我想要里面的数据,用findobj命令最后再用plot画图,但是画不出来,有哪位大神知道原因吗

您需要登录后才可以回帖 登录 | 注册账号

本版积分规则

QQ|网站地图|MATLAB技术论坛|Simulink仿真论坛 ( 蜀ICP备19014457号 

GMT+8, 2019-5-23 17:18 , Processed in 0.069764 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.4

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表