符号积分的一个奇怪的问题！

kevinhcj · 发表于 2006-3-14 15:03

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

我有如下循环 for i = 1:n C = C+element_C(N); end 其中 N为符号Xi的6×6的矩阵，element_C(N)主要就是对N在0～1之间积分，其余运算不涉及矩阵size的变化。但是改循环只能循环30～40次，而且各次执行其次数还不一致，出现的问题说计算失败那一次reshape时矩阵单元数目不一致！ 注; maple在做完符号积分后对结果有一次reshape（因为做之前将矩阵变成列向量，其后又还原为矩阵）,因此，这里的问题相当于6×6－－> 36×1－？－> 6×6 但是我在int.m＆maple.m里面读出被积函数的size一直没有变化 我在自己的计算中一致不涉及矩阵size的变化，难道这种计算又次数的限制（我用的是D版） 程序debuge到此，也有2天无法进行了，高手们给点建议。 thanx in advance!!

suffer · 发表于 2006-3-14 15:18

最好把详细的代码贴出来

kevinhcj · 发表于 2006-3-14 18:36

ok! %subrution function C_b1 = element_Cb2(z0,z1,t,bb) syms Xi psi = bb;
N1 = 1-10*Xi^3+15*Xi^4-6*Xi^5; N2 = 10*Xi^3-15*Xi^4+6*Xi^5; N3 = Xi-6*Xi^3+8*Xi^4-3*Xi^5; N4 = -4*Xi^3+7*Xi^4-3*Xi^5; N5 = 1/2*(Xi^2-3*Xi^3+3*Xi^4-Xi^5); N6 = 1/2*(1*Xi^3-2*Xi^4+Xi^5);
N = [N1,N3,N5,N2,N4,N6];
NN1 = diff(N,Xi,1); cc = (psi'*(NN1'*NN1)*psi)*(N'*N); C_b1 = int(cc,Xi,0,1); C_b1 = C_b1; C_b1 = double(C_b1);
return
% main clear all; clc; format long;
n = 10 m = 100; zz = [0,1:n]'; delta_t = 0.1;
C_b1_u = zeros(3*n+3); psi_q_try = ones(3*(n+1),1);
for j = 1:m t = delta_t*j; for k = 1:n bb = psi_q_try(3*k-2:3*(k+1)); C_b1 = element_Cb1(zz(k),zz(k+1),t,bb); C_b1_u(3*k-2:3*(k+1),3*k-2:3*(k+1)) = C_b1_u(3*k-2:3*(k+1),3*k-2:3*(k+1)) + C_b1; end end

happy · 发表于 2006-3-15 09:26

<DIV class=quote>以下是引用kevinhcj在2006-3-14 18:36:29的发言： 
ok! %subrution function C_b1 = element_Cb2(z0,z1,t,bb) syms Xi psi = bb;
N1 = 1-10*Xi^3+15*Xi^4-6*Xi^5; N2 = 10*Xi^3-15*Xi^4+6*Xi^5; N3 = Xi-6*Xi^3+8*Xi^4-3*Xi^5; N4 = -4*Xi^3+7*Xi^4-3*Xi^5; N5 = 1/2*(Xi^2-3*Xi^3+3*Xi^4-Xi^5); N6 = 1/2*(1*Xi^3-2*Xi^4+Xi^5);
N = [N1,N3,N5,N2,N4,N6];
NN1 = diff(N,Xi,1); cc = (psi'*(NN1'*NN1)*psi)*(N'*N); C_b1 = int(cc,Xi,0,1); C_b1 = C_b1; C_b1 = double(C_b1);
return
% main clear all; clc; format long;
n = 10 m = 100; zz = [0,1:n]'; delta_t = 0.1;
C_b1_u = zeros(3*n+3); psi_q_try = ones(3*(n+1),1);
for j = 1:m t = delta_t*j; for k = 1:n bb = psi_q_try(3*k-2:3*(k+1)); C_b1 = element_Cb1(zz(k),zz(k+1),t,bb); C_b1_u(3*k-2:3*(k+1),3*k-2:3*(k+1)) = C_b1_u(3*k-2:3*(k+1),3*k-2:3*(k+1)) + C_b1; end end </DIV>
红色标记的是不是应该是 C_b1 = element_Cb2(zz(k),zz(k+1),t,bb); 如果是的话，应该是你子程序中用int找不到积分

kevinhcj · 发表于 2006-3-15 10:46

thank u! 那里是笔误，对不起！我的子程序就是C_b1 = element_Cb1(zz(k),zz(k+1),t,bb); 这个问题的奇怪之处在于，输入参数的size一直没有变化，退一万步讲，如果参数size变化了，那就应该一个外循环也不能执行啊。而这个程序可以执行很多次循环而后出错。更有甚者，不改变任何设置，每一次运行的总循环次数n×（j－1）＋k不一致。所以我觉得是不是D版或者本身有循环数限制？ thanx again！

happy · 发表于 2006-3-16 10:18

<DIV class=quote>以下是引用kevinhcj在2006-3-15 10:46:00的发言： thank u! 那里是笔误，对不起！我的子程序就是C_b1 = element_Cb1(zz(k),zz(k+1),t,bb); 这个问题的奇怪之处在于，输入参数的size一直没有变化，退一万步讲，如果参数size变化了，那就应该一个外循环也不能执行啊。而这个程序可以执行很多次循环而后出错。更有甚者，不改变任何设置，每一次运行的总循环次数n×（j－1）＋k不一致。所以我觉得是不是D版或者本身有循环数限制？ thanx again！</DIV>
 你说的现象是很奇怪，不过我运行没有你说的总循环次数不一致的问题

kevinhcj · 发表于 2006-3-16 13:01

呵呵，这个问题另外一个高手也没有解决，不会是MATLAB的bug把？不过看来也不太像。

happy · 发表于 2006-3-17 07:46

<DIV class=quote>以下是引用kevinhcj在2006-3-16 13:01:52的发言： 呵呵，这个问题另外一个高手也没有解决，不会是MATLAB的bug把？不过看来也不太像。</DIV>
可以去mathwork问问看，那里专业，呵呵

kevinhcj · 发表于 2006-3-18 08:59

it's a good idea!

账号		自动登录	找回密码
密码			我要加入

符号积分的一个奇怪的问题！

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

回复：（kevinhcj）符号积分的一个奇怪的问题！

回复：（suffer）回复：（kevinhcj）符号积分的一个奇...

回复：（kevinhcj）回复：（suffer）回复：（kevinhcj...

回复：（kevinhcj）thank u!那里是笔误，对不起！我的...

回复：（kevinhcj）呵呵，这个问题另外一个高手也没有...