24小时热门版块排行榜    

查看: 1205  |  回复: 11
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

gengbiaolu

铜虫 (正式写手)

[交流] 【求助】求高人帮我优化程序,以提高速度。

我编的这个程序太需要时间了,我算了一个晚上也没有结果,求高人帮我优化程序,以提高速度,万分感谢!!!
clear all
global N c v
N=1000;c=1.8;v=1;
Y0=[1,zeros(1,N)]';
[t,YY]=ode45(@lgbquantum,[0,100],Y0);
b=abs(YY).*abs(YY);
s1=0
for n=0:N
   s1=s1+n*b(:,n+1);
end
plot(t,s1/N)
hold on
s2=0
  for n=0:N
    s2=s2+(N-n)*b(:,n+1);
  end  
plot(t,s2/N)

function Yd=lgbquantum(t,YY)
global N c v
n=0:N;
H1=c*(n.*(n-1)+(N-n).*(N-n-1))/(2*N);
H2=-v/2*sqrt(n.*(N-n+1));
Yd=-i*(diag(H1)+diag(H2(2:N+1),-1)+diag(H2(2:N+1),1))*YY;

[ Last edited by gengbiaolu on 2010-10-10 at 22:55 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijinfeng042

木虫 (小有名气)

Matlab

引用回帖:
Originally posted by gengbiaolu at 2010-10-13 22:19:30:
我要算的东西大概意思是:
H是矩阵,其对角元为:H(n,n)=(c/(2N))*[n(n-1)+(N-n)(N-n-1)]
非对角元:H(n,n-1)=H(n-1,n)= - (v/2)*[n(N-n+1)]^0.5
其它矩阵元为零。
通过时间的一价微分方程   i d(a_n(t))/dt=H ...

看你你的方程反而疑惑了
>> dsolve('Dy=H*y/i')

ans =

C2*(1/exp(H*t*i))

难道不是
还是
D(y1,y2...)=H y  ?????  N个方程?
工作了,偶尔会上来~可以关注新浪微博 @云是风的梦_Matlab
9楼2010-10-14 13:54:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 12 个回答

lijinfeng042

木虫 (小有名气)

Matlab


hiqun(金币+1):感谢专业应助 2010-10-10 12:15:39
这是一个不完整程序? ode87 没有这个函数的 是否你自己定义的? lgbquantum的w没定义
工作了,偶尔会上来~可以关注新浪微博 @云是风的梦_Matlab
2楼2010-10-10 12:06:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gengbiaolu

铜虫 (正式写手)

robert2020:建议虫友使用“引用回复该帖”,方便对方收到你的信息。 2010-10-11 12:56:34
就用ode45吧,其实ode87也有了,精度高些。w是多余的,已去掉了。

[ Last edited by gengbiaolu on 2010-10-10 at 22:57 ]
3楼2010-10-10 22:53:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijinfeng042

木虫 (小有名气)

Matlab


nono2009(金币+1):鼓励交流。 2010-10-12 16:18:22
gengbiaolu(金币+10):谢谢你的帮助!Matlab真的没有办法优化了吗? 2010-10-12 17:27:31
引用回帖:
Originally posted by gengbiaolu at 2010-10-10 22:53:37:
就用ode45吧,其实ode87也有了,精度高些。w是多余的,已去掉了。

[ Last edited by gengbiaolu on 2010-10-10 at 22:57 ]

唉 还是提高不了 100的时候测试 基本上是调用函数花时间的
ode方程不存在刚性问题
循环优化不了
如果没办法就建议把它做成C/C++ 的dll调用试试 我没有混合编程的能力
工作了,偶尔会上来~可以关注新浪微博 @云是风的梦_Matlab
4楼2010-10-12 13:41:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见