mfunction.m
function[P] = mfunction(v,M,R,T)
P=4*pi*(M/(2*pi*R*T)).^(3/2)*v.^2*exp((-M*(v.^2))/(2*R*T));
end
我想制作图表。 x为v,范围为1: 1200
,
M = 0.032,R = 8.31, T= 300
我想绘制y=mfunction(x)
并弹出错误。我该如何绘制图表?
答案 0 :(得分:2)
问题在于中间乘法
P=4*pi*(M/(2*pi*R*T)).^(3/2)*v.^2 * exp((-M*(v.^2))/(2*R*T));
^
你在做什么与做
类似[1,2] * [1,2]
,它会给出错误“mtimes inner matrix dimensions must agree”,因为它将mx1矩阵乘以mx1矩阵。
根据您的需要,您可以执行以下操作之一:
>> [1,2] * [1,2]' %inner product
ans =
5
>> [1,2]' * [1,2]
ans =
1 2
2 4
>> [1,2] .* [1,2] %element-wise product
ans =
1 4