MATLAB错误：实现公式和接收/不能使用，因为矩阵维度不一致错误

Question

我有一个简单的MATLAB程序，但是我遇到了错误。我试图将TS1，TS2和TS3作为t的函数从0到40进行绘制。但是，我遇到了各种错误。

我得到的错误之一是"/ can't be used because matrix dimensions do not agree"。我还附上了我应该实施的实际公式的图片。任何帮助将不胜感激。

%Ground Surface Temperature Plot

clc;
clear;

var t
var e1
var e2
var e3
var Ts1
var Ts2
var Ts3


t = linspace(0,40,1000);
Ti = 285.93;
Twind1 = 244.2611111;
Twind2 = 238.7055556;
Twind3 = 233.15;
x = 0; 
h = 16.3987;
a = 13.78579*10^(-6);
k = 0.52;
L = 60.96;

e1 = erfc(x/(2*sqrt(a*t*86400)));
e2 = exp(h*x/k + h^2*a*t*86400/k^2);
e3 = erfc(x/(2*sqrt(a*t*86400))+(h*sqrt(a*t*86400)/k));

Ts1 = (e1 - e2*e3)*(Twind1 - Ti) + Ti;
Ts2 = (e1 - e2*e3)*(Twind2 - Ti) + Ti;
Ts3 = (e1 - e2*e3)*(Twind3 - Ti) + Ti;

figure
Plot(t,Ts1,t,Ts2,t,Ts3);

正在编码的公式如下。

Answer 1

此行发生错误。

e1 = erfc(x/(2*sqrt(a*t*86400)));

这些是你试图划分的东西的大小。

>> size(2*sqrt(a*t*86400))
ans =
           1        1000
>> size(x)
ans =
     1     1

我相信您正在尝试进行元素明智的划分，因此请使用./，如下所示。

e1 = erfc(x./(2*sqrt(a*t*86400)));

请注意，无论何时在数字和矩阵之间执行元素分割或乘法运算，都需要执行此操作，因此请相应地修改e2和e3。