我想把一个符号矩阵数值。我能怎么做? ?

4视图(30天)
清晰的所有
clc
信谊x y
一个= [1 x ^ 2 + y ^ 2; x - y x ^ 2 + y ^ 2);
我想数字集成所有元素的x限制(清廉)和ylimits (0-15)。
谢谢你的时间和帮助。

接受的答案

沃尔特·罗伯森
沃尔特·罗伯森 2023年4月28日
使用嵌套vpaintegral()调用。内部调用将首先执行,找出它无法整合,因为有一个额外的自由变量,所以它将返回一个vpaintegral()的数据形式。外vpaintegral将识别vpaintegral数据形式的内在和能够进行2 d集成。
2的评论

登录置评。

更多的答案(2)

Torsten
Torsten 2023年4月27日
信谊x y
一个= [1 x ^ 2 + y ^ 2; x - y x ^ 2 + y ^ 2);
IntA = int (int (A、x 0, 10), y, 0, 15)
IntA =
1评论
不过朋友
不过朋友 2023年4月28日
我想做这种集成数字. .我有一个大的方程组在矩阵的每个元素的形成是不同的变量的函数。某些情况下出现的价值整合是恒定的。“int”操作很长对于我的情况。谢谢你不管怎样. .

登录置评。


保罗
保罗 2023年4月28日
我不认为有一个函数的数值积分一个arrayvalued两个变量的函数。这是一个循环的方法来整合单独使用的每个元素 integral2
信谊x y
一个= [1 x ^ 2 + y ^ 2; x - y x ^ 2 + y ^ 2);
2 = 1:2,
jj = 1:2
Afun = matlabFunction ((ii, jj),“var”,{' x ' ' y '});
%添加0 * x,以确保适当的尺寸当集成(1,1)
结果(ii, jj) = integral2 (@ (x, y) Afun x (x, y) + 0 * 0, 10 0, 15);
结束
结束
格式短e
结果
结果= 2×2
1.0 e + 00 * 1.5000 1.6250 e + e + 02年04 -3.7500 1.6250 e + e + 02年04
2的评论
保罗
保罗 2023年4月29日
如果我很惊讶 vpaintegral 是速度比 integral2 相同的公差。让我们试一试。
信谊x y
一个= [1 x ^ 2 + y ^ 2; x - y x ^ 2 + y ^ 2);
抽搐
kk = 1:10 0
结果=南(2,2);
2 = 1:2,
jj = 1:2
Afun = matlabFunction ((ii, jj),“var”,{' x ' ' y '});
%添加0 * x,以确保适当的尺寸当集成(1,1)
结果(ii, jj) = integral2 (@ (x, y) Afun x (x, y) + 0 * 0, 10 0, 15);
结束
结束
结束
toc
运行时间是4.600858秒。
抽搐
kk = 1:10 0
结果= vpaintegral (vpaintegral (A、x 0, 10), y, 0, 15);
结束
toc
运行时间是20.426502秒。

登录置评。

类别

找到更多的在数值积分和微分帮助中心文件交换

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!