| 查看: 1389 | 回复: 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 |
回复此楼» 猜你喜欢
中青基了要发朋友圈吗？
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
国社科申报系统有变化
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有6人回复
feibao(金币+3): 2010-04-18 12:16
feibao(金币+2): 2010-04-21 13:58
feibao(金币+2): 2010-04-21 13:58
|
f以及求导后的fv，fu在内存中都是一个变量的形式存储了一串很长的表达式 要求是存在连续的内存中 也许牵扯到在跨段和跨页的问题 在内存装载变量和赋值计算都要消耗内存 分开后可以分成几部分存 不要求是连续内存 计算也是分开的 这样算完一部分就能清空一部分 然后再加载后面的计算 实际对最大连续内存的需求不到原来1/24 单个变量所占内存太大处理起来可能出错 上面我做的处理没有把解析式加到一起再赋值 而是把每个解析式赋值 再把赋值的结果相加 [ Last edited by anyuezhiji on 2010-4-17 at 14:45 ] |
13楼2010-04-17 14:37:05
2楼2010-04-15 16:59:47
3楼2010-04-15 17:20:51
4楼2010-04-15 19:14:57
30