通过Maple中的for循环定义一组函数

阿德里亚人

在Maple中,我想通过两个for循环定义一组函数:

printlevel:=2;

# Node coordinates.
N_x:=5;
N_y:=4;
N_elx:=N_x-1;
N_ely:=N_y-1;
h_x:=(x_e-x_s)/N_elx;
h_y:=(y_e-y_s)/N_ely;
x_n:=[seq(x_s+j*h_x,j=0..N_elx)];
y_n:=[seq(y_s+j*h_y,j=0..N_ely)];

# Partition of unity.
for j from 2 by 1 to N_x-1 do
    for k from 2 by 1 to N_y-1 do
        phi[j,k]:=(x,y)->(x-x_n[j-1])*(x-x_n[j+1])*(y-y_n[k-1])*(y-y_n[k+1])/((x_n[j]-x_n[j-1])*(x_n[j]-x_n[j+1])*(y_n[k]-y_n[k-1])*(y_n[k]-y_n[k+1]));
    od;
od;

然而,这使输出其中[j,k]不同,但x_n[j-1]x_n[j+1]y_n[j-1]y_n[j+1]不被评估。因此,功能定义不正确。例如,使用调用j=4,k=2函数phi[4,2](x,y);,我得到输出

phi[4, 2] called with arguments: x, y
#(phi[4,2],1): (x-x_n[j-1])*(x-x_n[j+1])*(y-y_n[k-1])*(y-y_n[k+1])/(x_n[j]-x_n[j-1])/(x_n[j]-x_n[j+1])/(y_n[k]-y_n[k-1])/(y_n[k]-y_n[k+1])
Error, (in phi[4, 2]) invalid subscript selector

而不是所需的输出

$ \ phi_ {4,2}:= \ dfrac {(x-x_3)(x-x_5)(y-y_1)(y-y_3)} {(x_4-x_3)(x_4-x_5)(y_2-y_1) (y_2-y_3)} = \ dfrac {(x- \ dfrac {1} {2})(x-1)(y-0)(y- \ dfrac {2} {3})} {(\ dfrac { 3} {4}-\ dfrac {1} {2})(\ dfrac {3} {4} -1)/(\ dfrac {1} {3} -0)(\ dfrac {1} {3}- \ dfrac {2} {3})} $ 在此处输入图片说明

我该如何解决?

阿德里亚人

来自宏cer的回答:

如果要使用运算符而不是表达式,请使用函数unapply

N_x:=5;
N_y:=4;
N_elx:=N_x-1;
N_ely:=N_y-1;
h_x:=(x_e-x_s)/N_elx;
h_y:=(y_e-y_s)/N_ely;
x_n:=[seq(x_s+j*h_x,j=0..N_elx)];
y_n:=[seq(y_s+j*h_y,j=0..N_ely)];

# Partition of unity.
for j from 2 by 1 to N_x-1 do
    for k from 2 by 1 to N_y-1 do
        phi[j,k]:=unapply( (x-x_n[j-1])*(x-x_n[j+1])
                           *(y-y_n[k-1])*(y-y_n[k+1])
                           /((x_n[j]-x_n[j-1])*(x_n[j]-x_n[j+1])
                            *(y_n[k]-y_n[k-1])*(y_n[k]-y_n[k+1])),
                           [x,y]);
    od;
od;

本文收集自互联网,转载请注明来源。

如有侵权,请联系[email protected] 删除。

编辑于
0

我来说两句

0条评论
登录后参与评论

相关文章

来自分类Dev

如何在for循环中为R中的预定义结构的一组现有变量应用函数?

来自分类Dev

如何使用用户定义的相等函数在Javascript中创建一组用户定义的对象?

来自分类Dev

通过一组链接填充下拉值循环

来自分类Dev

在Maple中编写一个简单的用户定义函数

来自分类Dev

在Maple中编写一个简单的用户定义函数

来自分类Dev

在Common Lisp中定义一个宏,该宏定义一组函数和函数调用

来自分类Dev

R:定义一个函数(和/或使用apply()或for循环)以重复执行一组过程

来自分类Dev

数组中具有一组值的循环迭代

来自分类Dev

如何在活动中循环一组视图?

来自分类Dev

尝试在PowerShell中循环一组命令

来自分类Dev

在R中循环一组数据帧

来自分类Dev

在R中循环一组生成数字的命令

来自分类Dev

如何在活动中循环一组视图?

来自分类Dev

数组中具有一组值的循环迭代

来自分类Dev

循环浏览谷歌表格中的一组值

来自分类Dev

循环ggscatter函数,用于数据框列中x,y变量的唯一组合

来自分类Dev

遍历D中的一组自定义对象

来自分类Dev

遍历D中的一组自定义对象

来自分类Dev

在数组中依次调用一组函数

来自分类Dev

在数组中依次调用一组函数

来自分类Dev

使用一组函数作为javascript中的原型

来自分类Dev

在Excel中,对一组Lookup值的Lookup函数实现

来自分类Dev

定义一组循环视图以在同一布局中使用多个时间

来自分类Dev

创建一组带有定义时间定义的参数的匿名函数

来自分类Dev

创建一组带有定义时间定义的参数的匿名函数

来自分类Dev

一组中的总和

来自分类Dev

熊猫在一组上循环

来自分类Dev

无限循环一组图像

来自分类Dev

Typescript 循环遍历一组值

Related 相关文章

热门标签

归档