| 查看: 1395 | 回复: 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 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
找导师
已经有4人回复
-
我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
5楼2010-04-15 19:36:01
2楼2010-04-15 16:59:47
3楼2010-04-15 17:20:51
4楼2010-04-15 19:14:57
10