newuniverse 发表于 2006-9-7 00:10

[求助]请达人看看我的问题在哪里

function df = tov(r,f)
%f是P&m相对于r的函数
%y是能量密度
df=zeros(2,1);
G=6.67e-11;
c=2.9979e8;
global y;（这里y可以赋任一值）
df(1)=-(G*f(2)*y/r^2)*((f(1)/(c^2*y)+1)*(1+4*pi*r^3*f(1)/(c^2*f(2)))*(1-2*G*f(2)/(c^2*r))^-1);
df(2)=4*pi*r^2*y;
对上面那个常微分方程组进行ode45求解（y可以赋任一值）
结果出现下面的错误
??? Error using ==> mrdivide
Matrix dimensions must agree.

Error in ==> tov at 8
df(1)=-(G*f(2)*y/r^2)*((f(1)/(c^2*y)+1)*(1+4*pi*r^3*f(1)/(c^2*f(2)))*(1-2*G*f(2)/(c^2*r))^-1);请问我的问题在哪儿？

newuniverse 发表于 2006-9-7 00:23

shafa

xinyuxf 发表于 2006-9-7 08:17

矩阵维数不一致？
f具体是什么样形式的矩阵？

newuniverse 发表于 2006-9-7 17:57

问题已经解决,好象是全局变量必须用之前要声明,我的main.m中忘记声明了

bainhome 发表于 2006-9-7 22:56

全局变量跟这个错误根本两码事儿，不是根本症结，除非循环调用某子函数中计算出的变量，在循环中自己维数发生了变化。

[ 本帖最后由 bainhome 于 2006-9-7 23:02 编辑 ]

AaronSpark 发表于 2006-9-8 00:56

不是这个问题造成的，我前面声不声名全局变量结果都是一致的，没有出现错误

newuniverse 发表于 2006-9-8 14:36

那是怎么回事?

wangyongq246 发表于 2006-9-8 14:48

看不出来!!学习为主!!

happy 发表于 2006-9-8 15:01

原帖由 newuniverse 于 2006-9-8 14:36 发表
那是怎么回事?

估计是搂主其他地方搞错了，和上述代码无关
页: [1]
查看完整版本: [求助]请达人看看我的问题在哪里