MATLAB中算出各点曲率后如何返回到COMSOL？

hgniang

我在泡泡上种了好多点，然后在MATLAB中读取了点的坐标之后，算出每个点的曲率，
con=geominfo(PT1,'out','mp');
r=con(1,;）
z=con(2,;）
h=pi/2/49;
dr=gradient(r,h);
dz=gradient(z,h);
ddr=2*2*del2(r,h);
ddz=2*2*del2(z,h);
CHAFEN6=abs(dr.*ddz-dz.*ddr)./(dr.^2+dz.^2).^1.5;
得出的CHAFEN6这样一个数组如何加载到COMSOL中呢？

kyolee

可以尝试将曲率写成一个function，然后再comsol里面进行调用这个function，最后在边界上运用，注意矩阵的维度问题。等待你的结果!

hgniang

可以尝试将曲率写成一个function，然后再comsol里面进行调用这个function，最后在边界上运用，注意矩阵的维度问题。等待你的结果!
kyolee 发表于 2010-9-10 08:09

我有点疑惑呀，我是把这整个的四分之一个圆上所有的点的曲率算出来放在了矩阵中，那调用的时候该咋个写法呢？是每一条边都调用这个function ,还是只是赋矩阵中相应的那个维呢？

还有，我试着写成function，但是，有个geominfo()函数在这里，我要得到这些点的信息呀。

function chafen6=chafenu(PT1)
con=geominfo(PT1,'out','mp');
r=con(1,;)
z=con(2,;)
h=pi/2/49;
dr=gradient(r,h);
dz=gradient(z,h);
ddr=2*2*del2(r,h);
ddz=2*2*del2(z,h);
chafenu=abs(dr.*ddz-dz.*ddr)./(dr.^2+dz.^2).^1.5;

运行function之后，

??? Input argument " PT1" is undefined.
Error in ==> CHAFEN6 at 2
con=geominfo(PT1,'out','mp');

这个PT1该咋个定义呢？不好意思，我对MATLAB不熟。请指教哈

kyolee

我几乎不用调用，记忆中调用的时候，是写成一个function和数个argument的. 你的问题是否可以把argument写成是r和z呢，因为调用的时候是需要调用argument的。
写成这样 funciton out=chafen6=chafenu(r,z)
****
chafenu = ***

然后使用matlab function： 注意需要写入你的argument， 然后去边界上应用这个funciton。

你试试，看看可否实现！

期待你的结果！