查看单个帖子
旧 2019-12-10, 16:49   #1
poster
高级会员
 
注册日期: 2019-11-21
帖子: 3,006
声望力: 66
poster 正向着好的方向发展
帖子 如何在不使用任何循环的情况下生成两个变量的函数?

假设我有一个函数y(t,x)= exp(-t)* sin(x)

在Matlab中,我定义

t = [0: 0.5: 5]; x = [0: 0.1: 10*2*pi]; y = zeros(length(t), length(x)); % empty matrix init 现在,如何在不使用任何循环的情况下定义矩阵y,以使每个元素y(i,j)在(t(i), x(j))处包含所需函数y的值?以下是我使用for循环的方法。

for i = 1:length(t) y(i,:) = exp(-t(i)) .* sin(x); end
回答:
您的输入向量x是1xN , t是1xM ,输出矩阵y是MxN 。为了向量化代码, x和t必须与y具有相同的尺寸。

[x_,t_] = meshgrid(x,t); y_ = exp(-t_) .* sin(x_); 您的示例是一个简单的2D案例。函数meshgrid()也可用于3D。有时您无法避免循环,在这种情况下,当循环可以达到1:N或1:M时,请选择最短的一个。我用于准备向量化方程式(向量x矩阵乘法)的向量的另一个函数是diag() 。



更多&回答...
poster 当前离线   回复时引用此帖