matlab 计算效率,Mathematica和Matlab相比的计算效率问题

冷不防 2023-01-14 04:59 145阅读 0赞

该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

之前一直用Matlab，最近想学一些Mathematica。听说MMA的数值计算效率和Matlab不相上下，但实际测试感觉差别太大，不知该怎么优化。如下一个简单的例子，用Jacobi迭代法解Laplace方程，速度可相差20倍左右。用MMA的NDSolve解速度会快一些，但仍大幅度慢于Matlab，请教大家有什么好方法？

Mathematica

n = 32; (\*grid size\*)

h = 1/(n - 1); (\*grid spacing\*)

u =

Table\[0., \{n\}, \{n\}\];(\*init to zero\*)

u\[\[1, All\]\] = 100;(\*B.C.\*)

Do\[

u\[\[2 ;; n - 1, 2 ;; n - 1\]\] =

Table\[1/4 (u\[\[i - 1, j\]\] + u\[\[i + 1, j\]\] + u\[\[i, j - 1\]\] +

u\[\[i, j + 1\]\]), \{i, 2, n - 1\}, \{j, 2, n - 1\}\],

100\]

ListPlot3D\[u, PlotRange -> All, ImagePadding -> 20, Mesh -> \{n, n\},

ImageSize -> 300, PlotLabel -> "solution of Laplace PDE on 2D"\]

Matlab

clear all; close all;

n = 32;

h = 1/(n-1); %grid spacing

u = zeros(n,n);

u(1,:) = 100; %B.C.

% coordinates

\[X,Y\] = meshgrid(0:h:1,0:h:1);

for k=1:100

tmp = u;

for i=2:n-1

for j=2:n-1

tmp(i,j)=...

(1/4)\*(u(i-1,j)+u(i+1,j)+...

u(i,j-1)+u(i,j+1));

end

end

u=tmp;

end

mesh(X,Y,u);

title('solution of Laplace PDE on 2D');

set(gcf,'Position',\[10,10,320,320\]);

发表评论取消回复

表情：

评论列表（有 0 条评论，145人围观）

还没有评论，来说两句吧...

相关阅读

相关使用MATLAB计算平均值

使用MATLAB计算平均值计算数学中最基本的操作之一就是计算一组数字的平均值。在MATLAB中，计算平均值是非常简单的，只需使用mean函数即可。下面是一个简单的例子展示如

矫情吗；*/ 2024年03月17日 17:20/ 0 赞/ 42 阅读

相关 MATLAB：计算数组长度

MATLAB：计算数组长度在MATLAB中，计算数组的长度是一项常见的操作。数组长度是指数组中元素的个数。在本文中，我们将讨论如何使用MATLAB计算数组的长度，并提供相应

港控/mmm°/ 2024年03月04日 03:56/ 0 赞/ 115 阅读

相关 Matlab 方位角计算

Matlab 方位角计算方位角是描述一个点相对于参考方向的角度。在Matlab中，可以使用一些数学函数和向量运算来计算方位角。本文将介绍如何使用Matlab计算方位角，并提

淡淡的烟草味﹌/ 2024年03月03日 08:44/ 0 赞/ 113 阅读

相关 MATLAB-向量相关计算

MATLAB 向量类型：行向量列向量 MATLAB 行向量：创建行向量括在方括号中的元素的集合，用空格或逗号分隔的元素。 r =

朱雀/ 2023年01月23日 03:56/ 0 赞/ 151 阅读

相关 matlab 集群计算,Matlab 集群计算使用心得

相信很多研究人员都和我一样使用Matlab分析和研究各种算法。最近无意之中发现matlab提供的并行工具箱，通过简单地修改(仅仅把for改为parfor，存在依赖关系的变量用新

偏执的太偏执、/ 2023年01月15日 12:30/ 0 赞/ 181 阅读

相关 matlab 计算效率,Mathematica和Matlab相比的计算效率问题

该楼层疑似违规已被系统折叠隐藏此楼查看此楼之前一直用Matlab，最近想学一些Mathematica。听说MMA的数值计算效率和Matlab不相上下，但实际测试感觉差别太

冷不防/ 2023年01月14日 04:59/ 0 赞/ 146 阅读

相关 MATLAB07:数值计算

转自：[https://blog.csdn.net/ncepu\_Chen/article/details/103165266][https_blog.csdn.net_nce

╰+攻爆jí腚メ/ 2022年12月04日 02:22/ 0 赞/ 302 阅读

相关 Matlab--符号计算

<table> <thead> <tr> <th align="left">函数</th> <th align="left">函数功能</th>

短命女/ 2022年08月03日 05:14/ 0 赞/ 242 阅读

相关 MatLab矩阵计算

取名MATLAB即Matrix Laboratory 矩阵实验室的意思。例如：在MATLAB命令窗口输入命令： a=\[1,1.5,2,9,7; 0,3.6,0.5,

系统管理员/ 2022年04月10日 08:22/ 0 赞/ 341 阅读

相关如何分析matlab程序的主要效率问题

利用profile on 在需要分析效率的程序段前后加入 profile on profile off 然后，在common line中输入profile viewer

刺骨的言语ヽ痛彻心扉/ 2021年07月16日 23:52/ 0 赞/ 294 阅读