24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1132  |  回复: 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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 266分,求材料相关专业调剂 +12 哇呼哼呼哼 2026-03-30 14/700 2026-04-03 14:11 by 欣喜777
[硕博家园] 求老师收留 +9 lllq123 2026-04-03 9/450 2026-04-03 13:48 by 呼吸都是减肥
[考研] 求调剂,一志愿郑州大学材料与化工专硕,英二数二342分,求老师收留 +11 v12abo 2026-04-02 12/600 2026-04-03 11:33 by chenqifeng666
[考研] 324求调剂 +12 想上学求调 2026-04-02 12/600 2026-04-03 11:15 by linyelide
[考研] 材料考研调剂 +10 Gs大王 2026-04-02 10/500 2026-04-03 09:47 by 遗忘消失的灆
[考研] 279求调剂 +5 傅文秋 2026-04-02 5/250 2026-04-02 18:10 by 笔落锦州
[考研] 282求调剂 +18 ycy1201 2026-04-01 20/1000 2026-04-02 16:39 by liumengping
[考研] 找调剂 +6 楚乔乔 2026-04-01 6/300 2026-04-02 12:06 by xiaoranmu
[考研] 070300化学求调剂 +14 小黄鸭宝 2026-03-30 14/700 2026-04-02 10:07 by oooqiao
[考研] 367求调剂 +8 芋泥啵啵… 2026-03-28 8/400 2026-04-02 09:20 by olim
[考研] 286求调剂 +5 Sa67890. 2026-04-01 7/350 2026-04-01 19:50 by 6781022
[考研] 求0861交通运输专硕or材料专硕调剂 +4 勒布朗@ 2026-03-31 4/200 2026-04-01 09:54 by 一只好果子?
[考研] 0856调剂 +7 曲听筠 2026-03-30 7/350 2026-04-01 08:51 by JourneyLucky
[考研] 材料工程专硕求调剂 +10 hyl3153942 2026-03-29 10/500 2026-03-31 16:31 by hypershenger
[考研] 物理学调剂 +4 小羊36 2026-03-30 4/200 2026-03-31 16:16 by lishahe
[考研] 一志愿大连理工大学材料求调剂 +6 Gymno 2026-03-30 6/300 2026-03-31 07:26 by 无际的草原
[考研] 348求调剂 +6 小懒虫不懒了 2026-03-28 6/300 2026-03-30 10:29 by Evan_Liu
[考研] 343求调剂 +6 爱羁绊 2026-03-29 6/300 2026-03-29 12:00 by 无际的草原
[考研] 085405 考的11408求各位老师带走 +3 Qiu学ing 2026-03-28 3/150 2026-03-28 09:19 by 乐呵呵的追梦人
[考研] 272求调剂 +7 脚滑的守法公民 2026-03-27 7/350 2026-03-27 17:23 by laoshidan
信息提示
请填处理意见