24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1138  |  回复: 14
当前主题已经存档。
本帖产生 1 个 仿真EPI ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

feibao

银虫 (小有名气)

[交流] 【求助】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
回复此楼

» 猜你喜欢

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

anyuezhiji

银虫 (正式写手)

星空行者

feibao(金币+5): 2010-04-18 12:13
Out of memory. Type HELP MEMORY for your options

内存不足
引用回帖:
Matlab 中数组必须占用连续分配的内存段,当无法为新建的数组分配连续的内存段的时候,”Out of Memory” 就会出现。由于反复分配和释放数组会使可用的连续内存段减少,因此当 Matlab 刚刚启动时其连续内存最多,此时往往可以新建非常大的数组,这一点可以用命令 feature(’memstats’)(在 7.0 版本以上)看出。如果现实的最大连续内存段很小,但实际可用内存(非连续的)仍旧很多,则表明内存中碎片太多了。此时可以考虑用 pack 命令,pack 命令的作用就是将所有内存中的数组写入硬盘,然后重新建立这些数组,以减少内存碎片。此外,在命令行或者程序中都可以使用 clear 命令,随时减少不必要的内存。

你用diff对符号函数微分时可能占很大的内存 而你的内存已经无法为新建的数组分配连续的内存段
建议符号函数计算出结果简化下
尤其在用diff算微分前最好先把f简化下 如f=simple(f)
另外释放掉一些已经可以不用的变量 如 clear x1 x2

引用回帖:
out of Memory问题处理方法一般如下:

1. 增加虚拟内存
2. 采用PACK
3. 条件允许的话采用64位系统
4. 采用3GB switch,增加虚拟内存
5. 优化程序,减少变量
6. 如果必有必要,不要启动java虚拟机,采用matlab -nojvm启动
7. 关闭Matlab Server

[ Last edited by anyuezhiji on 2010-4-15 at 20:08 ]
暗月下没有留下风的痕迹,但它已经寂然飘逝。。By<暗月之寂>:tiger38:
4楼2010-04-15 19:14:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 15 个回答

lingwei004

金虫 (小有名气)

feibao(金币+1): 2010-04-18 12:12
数据太大了吧。你0.1.txt里面有多少个数据?估计是这个地方问题
2楼2010-04-15 16:59:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feibao

银虫 (小有名气)

引用回帖:
Originally posted by lingwei004 at 2010-04-15 16:59:47:
数据太大了吧。你0.1.txt里面有多少个数据?估计是这个地方问题

大约100多个数据。我若取其中一个数据赋值v=7.9842590e-001   u= 7.7820931e-001 仍然不起作用  谢谢 请帮忙指点指点
3楼2010-04-15 17:20:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anyuezhiji

银虫 (正式写手)

星空行者

feibao(金币+2): 2010-04-18 12:13
测试了下

上面的代码在用diff算微分前
f的表达式的长度为927553字节

MATLAB的表达式非汉字 一个字节就对应一个字符
f的表达式长度就有92万多字符

暗月下没有留下风的痕迹,但它已经寂然飘逝。。By<暗月之寂>:tiger38:
5楼2010-04-15 19:36:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 26调剂 086003 +3 失活的细胞 2026-04-04 3/150 2026-04-04 02:24 by 毕业摆摊卖煎饼
[考研] 266求调剂 +6 学员97LZgn 2026-04-03 6/300 2026-04-04 00:00 by userper
[考研] 一志愿北京科技大学材料工程085601,求调剂 +14 cdyw 2026-04-02 15/750 2026-04-03 21:49 by qlm5820
[考研] 338求调剂 +4 zzz,,r 2026-04-03 4/200 2026-04-03 16:39 by lijunpoly
[考研] 化学调剂 +9 艾志恒 2026-04-03 10/500 2026-04-03 15:51 by 晴空210210
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +9 哇呼哼呼哼 2026-04-02 9/450 2026-04-03 12:05 by 1753564080
[考研] 一志愿安徽大学0817化学工程与技术,求调剂 +14 我不是只因 2026-04-02 15/750 2026-04-03 09:49 by 蓝云思雨
[考研] 复试调剂 +3 bvzz 2026-04-01 3/150 2026-04-03 09:47 by 蓝云思雨
[考研] 295求调剂 +7 愿旅途永远坦然 2026-04-02 7/350 2026-04-03 08:22 by fangshan711
[考研] 材料340分调剂 +7 夏夜晚风_long 2026-04-02 9/450 2026-04-02 21:20 by dongzh2009
[考研] 264分,求任意工科调剂 +4 zzlqwq 2026-03-29 5/250 2026-04-02 17:17 by 何曾几何
[考研] 304求调剂 +12 素年祭语 2026-03-31 15/750 2026-04-01 22:41 by peike
[考研] 379求调剂 +3 ?苦瓜不苦 2026-04-01 3/150 2026-04-01 20:09 by 暮云清寒
[考研] 265求调剂 +11 yelck 2026-04-01 12/600 2026-04-01 19:12 by 549790059
[考研] 生物与医药考研调剂 +5 铁憨憨123425 2026-03-31 5/250 2026-04-01 18:01 by syh9288
[考研] 调剂 +5 好好读书。 2026-03-28 7/350 2026-04-01 15:32 by 王亮_大连医科大
[考研] 08工科,295,接受跨专业调剂 +6 lmnlzy 2026-03-31 6/300 2026-04-01 11:02 by 逆水乘风
[考研] 318求调剂 +8 七忆77 2026-04-01 8/400 2026-04-01 10:37 by Jaylen.
[有机交流] 甲基亚磺磺酸钠和甲基磺酸酯反应机理 10+3 kaobao456 2026-03-29 4/200 2026-03-30 23:16 by nBu锂
[考研] 356求调剂 +3 gysy?s?a 2026-03-28 3/150 2026-03-29 00:33 by 544594351
信息提示
请填处理意见