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。