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

MATLAB技术论坛

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

[提问] 有偿,求大神帮忙,隐函数绘图

[复制链接]
发表于 2017-9-28 22:12:20 | 显示全部楼层 |阅读模式
40贝壳
clear
% clc
close all
syms p1 p2 r v o x c beta O X k w;% v表示 泊松比,  c表示b/a

c=1;
beta=0.195;
v=0.5;
H=[(p1^2-v)*cosh(p1*x) -(p2^2+v)*cos(p2*x); p1*(p1^2-o)*sinh(p1*x) p2*(p2^2+o)*sin(p2*x)]
eq1=det(H);
O=2-v-(1-beta)*k;
X=(pi/2)*c;
R=w^2/pi^4;
P1=(sqrt(sqrt((1-k/2)^2+(r-1))+(1-k/2)));
P2=(sqrt(sqrt((1-k/2)^2+(r-1))-(1-k/2)));
eq2 = subs(eq1,[p1 p2],[P1 P2])
eq3 = subs(eq2,[o x r],[O X R])
ezplot(eq3,[0 6 0 35]);
view(90,-90)
然后P2这个方程里面可能会小于零,函数虚数部分什么的,我画出来的是图1,然后我在P2=real(sqrt(sqrt((1-k/2)^2+(r-1))-(1-k/2)));加一个real提取实部画出来的是图2,正确的是图3红色线标出的那部分,图形还是有差别的,我搞不明白是不是实数部分提取的不对?我MATLAB不怎么好,大神帮帮忙呗,解决了还有偿

图1

图1

图2

图2

图3

图3

发表于 2017-9-29 11:47:54 | 显示全部楼层
你觉得你应该把问题描述清楚,最好有原始的书面公式,如果单纯只看你的程序,没有大的问题,不过你的程序是不是与书上的意思完全表达一致,那就另当别论了
回复

使用道具 举报

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

本版积分规则

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

GMT+8, 2019-8-21 23:53 , Processed in 0.098795 second(s), 12 queries , Gzip On, MemCache On.

Powered by Discuz! X3.4

© 2001-2017 Comsenz Inc.

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