| 查看: 779 | 回复: 4 | |||
jiasen1980金虫 (正式写手)
|
[交流]
【求助】那位高手帮看下MATLAB程序出错了
|
|
%FWT on the NLS equation i*du/dt+alpha*d^2u/dx^2+nu*|u|^2*u=0 %x=linspace(-3,3,1000);u=exp(2*x.^2);[c,L]=wavedec(u,3,'sym6');A=wrcoef('a',c,L,'sym6',3); %plot(x,u,'x',x,A,'r') %-----Claculate the wavelet------- xmax=16; level=3; alpha=1/2; nu=1; deltat=5*10^(-5); maxtime=2; iter=3; wav='sym6'; [phi,psi,xval]=wavefun(wav,iter); hval=xval(2)-xval(1); MX=max(xval); s=1; % mm=0; % jump=10^(s+1); flag=1; for t=0:deltat:maxtime if flag==1 flag=0; thresh=2^(level-3)*10^-5; phin=2^(level/2)*phi; psin=2^(level/2)*psi; lenphi=length(phin); shift=floor(lenphi/MX); N=xmax*2^level-(MX-1); lenPHI=(N-1)*shift+lenphi; x=linspace(-xmax/2,xmax/2,lenPHI); h=x(2)-x(1); %--------------Clculate the initial approximation of u-------------- if t==0 ua=3*sech(x); end t=0; for L=1:N aold(1,L)=h*ua((L-1)*shift+1  L-1)*shift+lenphi)*phin';end for L=1:N b(1,L)=h*ua((L-1)*shift+1 L-1)*shift+lenphi)*psin';end u0=zeros(1,lenPHI); for L=1:N u0(1,(L-1)*shift+1 L-1)*shift+lenphi)=u0(1,(L-1)*shift+1L-1)*shift+lenphi)+aold(1,L)*phin;end figure plot(x,ua,'k',x,u0,'r'); title(strcat('xmax=',num2str(xmax),',alpha=',num2str(alpha),',\nu=',num2str(nu),',maxtime=',num2str(maxtime),',deltat=',num2str(deltat),',level=',num2str(level),',iter=',num2str(iter),',time=',num2str(t))); %----Calculate KM Matrix for sym6 load sym6KMburg3; for L=1:N for K=1:N if abs(L-K) end end end %------------- load sym6LambdaNLS3; end % U=zeros(floor(maxtime/(jump*deltat))+1,lenPHI); %----Calculate C Matrix----------------- C=zeros(N,N); for L=1:N for K=1:N if abs(L-K) if abs(m-K) if abs(n-K) end end end end end end end % % for L=1:N % % for K=1:N % % if abs(L-K) % % end % % end % % end % %a=a+i*a*Mat; anon=aold+i*deltat*aold*nu*C; alin=anon+i*deltat*anon*alpha*KM; anew=alin; aold=anew; if max(abs(aold))>50 u=zeros(1,lenPHI); for L=1:N u(1,(L-1)*shift+1 L-1)*shift+lenphi)=u(1,(L-1)*shift+1L-1)*shift+lenphi)+abs(aold(1,L))*phi;end figure plot(x,u0,'r',x,u,'k'); title(strcat('xmax=',num2str(xmax),',alpha=',num2str(alpha),',\nu=',num2str(nu),',maxtime=',num2str(maxtime),',deltat=',num2str(deltat),',level=',num2str(level),',iter=',num2str(iter),',time=',num2str(t))); break; end u=zeros(1,lenPHI); for L=1:N u(1,(L-1)*shift+1 L-1)*shift+lenphi)=u(1,(L-1)*shift+1L-1)*shift+lenphi)+abs(aold(1,L))*phin;end for L=1:N b(1,L)=h*u((L-1)*shift+1 L-1)*shift+lenphi)*psin';end if max(abs(b))>thresh max(abs(b)); level=level+1; flag=1; ut=interp(u,2); ua=ut(1:length(ut)-1); % a=interp(a,2); end % if t == (mm*jump*deltat) % u=zeros(1,lenPHI); % for L=1:N % u(1,(L-1)*shift+1 L-1)*shift+lenphi)=u(1,(L-1)*shift+1L-1)*shift+lenphi)+abs(aold(1,L))*phin;% end % mm=mm+1; % U(mm,  =u;% end if t==deltat*10^(s+1)|t==deltat*10^(s+2)|t==deltat*10^(s+3)|t==deltat*10^(s+4)|t==0.5|t==0.75|t==1.5|t==1.75|t==2u=zeros(1,lenPHI); for L=1:N u(1,(L-1)*shift+1 L-1)*shift+lenphi)=u(1,(L-1)*shift+1L-1)*shift+lenphi)+abs(aold(1,L))*phin;end figure plot(x,u0,'r',x,u,'k'); title(strcat('xmax=',num2str(xmax),',alpha=',num2str(alpha),',\nu=',num2str(nu),',maxtime=',num2str(maxtime),',deltat=',num2str(deltat),',level=',num2str(level),',iter=',num2str(iter),',time=',num2str(t))); end end % figure % S=size(U); % % tp=[0:1:S(1)-1]; % tp=linspace(0,maxtime,S(1)); % [X,T]=meshgrid(x,tp); % mesh(X,T,U); % title(strcat('xmax=',num2str(xmax),',alpha=',num2str(alpha),',\nu=',num2str(nu),',maxtime=',num2str(maxtime),',deltat=',num2str(deltat),',level=',num2str(level),',iter=',num2str(iter),',time=',num2str(t))); ??? Error using ==> load Unable to read file sym6KMburg3: No such file or directory. |
回复此楼» 猜你喜欢
一志愿中南化学337求调剂
已经有5人回复
-
281求调剂（0805）
已经有23人回复
-
265求调剂
已经有9人回复
-
279分求调剂 一志愿211
已经有10人回复
-
招收调剂硕士
已经有12人回复
-
296求调剂
已经有8人回复
-
招收博士1-2人
已经有3人回复
-
085700资源与环境308求调剂
已经有11人回复
-
274求调剂
已经有8人回复
-
290求调剂
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助matlab一程序运行的问题，帮忙看看哪里不对
已经有4人回复
- matlab程序运行出错，求指点修改
已经有4人回复
- 请教关于matlab矩阵转换问题
已经有6人回复
- 用matlab求解方程出问题，请帮忙看看
已经有3人回复
- 关于Matlab作图的问题
已经有14人回复
- matlab调试问题
已经有3人回复
- 求助一简单matlab程序的出错问题，谢谢
已经有3人回复
- matlab程序运行总是出现下面的错误，请各位高手指点一下，谢谢了！急 急急
已经有5人回复
- 采用matlab中nlinfit拟合，出现问题了，不知道是哪错了，该怎么解决，求高手指点一下
已经有9人回复
- matlab中运行提示出错，请求高手指点
已经有8人回复
- 【求助】matlab程序运行出错，是怎么回事啊【已完结】
已经有6人回复
- 【求助】大家看看，matlab安装过程的截图，找出错
已经有14人回复
lixuke2005
铁杆木虫 (著名写手)
- 应助: 12 (小学生)
- 金币: 307.1
- 散金: 5990
- 红花: 11
- 帖子: 2214
- 在线: 1583.9小时
- 虫号: 548552
- 注册: 2008-04-19
- 专业: 艺术学其他学科
2楼2010-09-08 20:32:07
3楼2010-09-09 09:41:00
4楼2010-09-09 10:33:48
yqx1985
木虫 (著名写手)
云中仙
- 应助: 34 (小学生)
- 贵宾: 0.032
- 金币: 1608.8
- 散金: 1182
- 红花: 33
- 帖子: 1544
- 在线: 620.3小时
- 虫号: 550008
- 注册: 2008-04-22
- 性别: GG
- 专业: 药物设计与药物信息
5楼2010-09-09 17:46:01
