| 查看: 1343 | 回复: 14 | |||
| 当前主题已经存档。 | |||
| 本帖产生 1 个 仿真EPI ，点击这里进行查看 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】MATLAB程序出错求助 已有5人参与
|
|||
|
请高手指教 多谢 下面这个程序在MATLAB7.0版本中运行时总是出现如下问题： ??? Error using ==> minus Out of memory. Type HELP MEMORY for your options. Error in ==> sym.sym>char2sym at 395 p = cumsum((S == '(') - (S == ')')); Error in ==> sym.sym at 92 S = char2sym(x); Error in ==> sym.maple at 92 result = sym(result); Error in ==> sym.diff at 50 R = reshape(maple('map','diff',S(  ,x),size(S));但在MATLAB7.10版本中没有显示出错，却一直不出结果 不知怎么回事 请帮忙了 clear all >> clear all format short syms u v m=1; s=1+0.1; a=0.98; b=0.99; n1=1.48; n2=-a*n1; n3=b*n1; r=(n1^2-n2^2)/(n1^2-n3^2) w1=v^2*r-u^2; w=v^2-u^2; ne=sqrt((w^2/v^2)*(n1^2-n3^2)+n3^2); F1=-besselj(m+1,u)+m./u.*besselj(m,u); F=F1./(u.*besselj(m,u)); M1=besseli(m+1,w1)+m./w1.*besseli(m,w1); M=M1./(w1.*besseli(m,w1)); N1=-besselk(m+1,w1)+m./w1.*besselk(m,w1); N=N1./(w1.*besselk(m,w1)); Q=m*(1./u.^2+1./w1.^2); R=m*(1./w1.^2-1./w.^2)/s; S1=besseli(m+1,w1*s)+m./w1./s.*besseli(m,w1*s); S=S1./(w1.*besseli(m,w1*s)); T1=-besselk(m+1,w*s)+m./w./s.*besselk(m,w*s); T=T1./(w.*besselk(m,w*s)); X1=-besselk(m+1,w1*s)+m./w1./s.*besselk(m,w1*s); X=X1./(w1.*besselk(m,w1*s)); I1m=besseli(m,w1); K1m=besselk(m,w1); I2m=besseli(m,w1*s); K2m=besselk(m,w1*s); x1=Q.*ne.*I1m x2=Q.*ne.*K1m; x3=I1m.*(F-M); x4=K1m.*(F-N); x5=R*ne.*I2m; x6=R.*ne.*K2m; x7=I2m.*(-S-T); x8=K2m.*(-X-T); x9=I1m.*(n1^2.*F-n2^2.*M); x10=K1m.*(n1^2.*F-n2^2.*N); x11=Q.*ne.*I1m; x12=Q.*ne.*K1m; x13=I2m.*(-S.*n2^2-n3^2.*T); x14=K2m.*(-n2^2.*X-n3^2.*T); x15=R.*ne.*I2m; x16=R.*ne.*K2m; f=det([x1,x2,x3,x4;x5,x6,x7,x8;x9,x10,x11,x12;x13,x14,x15,x16]); fv=diff(f,v); fu=diff(f,u); f1=-fv/fu; f2=1+(u/v)^2*(1-2*v/u*f1); Q=load('e:\0.1.txt'); for i=1:length(Q(:,1)); v=Q(i,1); u=Q(i,3) f22=eval(f2) f11=eval(f1) hold on plot(v,f11,f22) end |
回复此楼» 猜你喜欢
今天系统多次维护，明天很可能放榜！
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有74人回复
-
岩土工程学报什么时候才能终审完呐
已经有4人回复
-
时间戳又变了8-15
已经有25人回复
-
2027广东省杰青
已经有6人回复
-
快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
3楼2010-04-15 17:20:51
2楼2010-04-15 16:59:47
4楼2010-04-15 19:14:57
5楼2010-04-15 19:36:01
25
