24小时热门版块排行榜    

查看: 393  |  回复: 3

1053507209

新虫 (初入文坛)

[求助] 画图求助

Sample Textfigure
>> nn=2:100;
>> tt=0:0.2:20;
>> [n,t]=meshgrid(nn,tt);
>> r=0.2;
>> x=0;
>> Om=1;
>> g=5;
>> Z1=gamma(2*n+1);
>> Z2=gamma(n+1);
>> Z3=(-exp(i*x)*tanh(r)).^n.*sqrt(Z1);
>> Z4=Z2.*2.^n.*sqrt(cosh(r));
>> f=Z3./Z4;
>> a=1/2*g*sqrt((Om/g).^2+8*(2*n+1));
>> b=1/2*g*sqrt((Om/g).^2+8*(2*n-1));
>> y1=exp(-i*Om.*t/2);
>> y2=sin(a.*t);
>> y3=sin(b.*t);
>> y4=cos(a.*t);
>> y5=cos(b.*t);
>> F1=-i*2*g*sqrt(n+1).*f.*y1.*y2./(sqrt(3)*a);
>> X1=n.*Om.*f.*y1.*(i*y3+2*b.*y5./Om)./(2*sqrt(3)*b.*(2*n-1));
>> X2=(n-1).*f./(sqrt(3)*(2*n-1));
>> F2=X1+X2;
>> F3=y1.*(f.*y4./sqrt(3)-i*Om.*f.*y2./(2*sqrt(3)*a));
>> F5=-2*i*g*sqrt(n).*f.*y1.*y2./(sqrt(3)*a);
>> F6=-i*g*sqrt(n).*f.*y1.*y3./(sqrt(3)*b);
>> X3=sqrt(n.*(n-1)).*Om.*f.*y1.*(i*y3+2*b.*y5./Om)./(2*sqrt(3)*b.*(2*n-1));
>> X4=f.*sqrt(n.*(n-1))./(sqrt(3)*(2*n-1));
>> F7=X3+X4;
>> H4= com(conj(F5).*F5+conj(F7).*F7-conj(F2).*F2-conj(F1).*F1)
>> plot(tt,H4)
我画的曲线图怎么也画不出?哪位高手能指点下啊?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ausir

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
请问H4的com()的定义是什么? 没有这个定义,什么都做不了。

在H4之前,全部F1,F2,F5,F7都正常。 H4的com里面的运算也没有问题。

因此只要给出com的定义,就可以正确完成绘图。
2楼2013-05-26 10:12:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

忆如沉

银虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
你的那个COM没有定义,还有,你是要画三维的图形吗,怎么感觉像,要是画三维的,应该是plot3.或者是suif,mesh,吧,希望可以帮到你啊
你若不离不弃,我便生死相依。
3楼2013-05-26 12:28:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

1053507209

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by ausir at 2013-05-26 10:12:11
请问H4的com()的定义是什么? 没有这个定义,什么都做不了。

在H4之前,全部F1,F2,F5,F7都正常。 H4的com里面的运算也没有问题。

因此只要给出com的定义,就可以正确完成绘图。

谢谢,是sum,写错了,是对n求和;画的图应该具有周期性,我得的并没不是,全是直线。
4楼2013-05-27 10:58:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 1053507209 的主题更新
信息提示
请填处理意见