声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1400|回复: 5

请有时间的高手帮忙运行一下代码,thanx!

[复制链接]
发表于 2006-3-15 18:00 | 显示全部楼层 |阅读模式

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

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

x
%subrutine<BR>function C_b1 = element_Cb1(z0,z1,t,bb)<BR>syms Xi <BR>psi = bb;
<P>N1 = 1-10*Xi^3+15*Xi^4-6*Xi^5;<BR>N2 = 10*Xi^3-15*Xi^4+6*Xi^5;<BR>N3 = Xi-6*Xi^3+8*Xi^4-3*Xi^5;<BR>N4 = -4*Xi^3+7*Xi^4-3*Xi^5;<BR>N5 = 1/2*(Xi^2-3*Xi^3+3*Xi^4-Xi^5);<BR>N6 = 1/2*(1*Xi^3-2*Xi^4+Xi^5);</P>
<P>N  = [N1,N3,N5,N2,N4,N6];</P>
<P>NN1 = diff(N,Xi,1);<BR>cc = (psi'*(NN1'*NN1)*psi)*(N'*N);<BR>C_b1 = int(cc,Xi,0,1);<BR>C_b1 = C_b1;<BR>C_b1 = double(C_b1);</P>
<P>return</P>
<P>% main<BR>clear all; clc;<BR>format long;</P>
<P>n = 10<BR>m = 100;<BR>zz = [0,1:n]';<BR>delta_t = 0.1;</P>
<P>C_b1_u = zeros(3*n+3);<BR>psi_q_try = ones(3*(n+1),1);</P>
<P>for j = 1:m<BR>    <BR>    t = delta_t*j;<BR>       <BR>        for k = 1:n<BR>            bb = psi_q_try(3*k-2:3*(k+1));<BR>            C_b1 = element_Cb1(zz(k),zz(k+1),t,bb);<BR>            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;<BR>        end <BR>end<BR><BR>在同一设置下运行多次,所能执行的循环数不一致!</P>
回复
分享到:

使用道具 举报

发表于 2006-3-15 20:23 | 显示全部楼层

回复:(kevinhcj)请有时间的高手帮忙运行一下代码,...

运行到j和k都等于9的时候停止
发表于 2006-3-15 20:27 | 显示全部楼层
不好意思楼主 我还想问一下 你说的那个矩阵积分的方法 算4阶矩阵的时候 10分钟都积不出来 是这样的么
 楼主| 发表于 2006-3-15 21:39 | 显示全部楼层
int(FUN,X,a,b)<BR>其中FUN可以是矩阵,使用就跟单个函数一样。至于所需的时间,看矩阵的内容而定了,要不你重新发个帖子,把代码帖出来,让我看看。
发表于 2006-3-16 09:57 | 显示全部楼层

回复:(kevinhcj)请有时间的高手帮忙运行一下代码,...

<DIV class=quote><B>以下是引用<I>kevinhcj</I>在2006-3-15 18:00:47的发言:</B><BR>
<P><BR>在同一设置下运行多次,所能执行的循环数不一致!</P></DIV>
<P><BR>试了两次,结果都和2楼一样,不存在循环数不一致问题</P>
发表于 2006-3-16 10:00 | 显示全部楼层

回复:(kevinhcj)请有时间的高手帮忙运行一下代码,...

另外有一个疑问,你的子程序中根本和z0,z1,t这几个变量无关,为什么还要传递呢?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-26 10:17 , Processed in 0.059295 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表