How to fix matrix dimensions

10 次查看（过去 30 天）
Jonathon Klepatzki
Jonathon Klepatzki 2024-4-10
移动：Rik 2024-4-11
Ran in:
Given the specific values of y, sigmay, and sigmaz (see code), I must calculate them and plot each line on a graph
close all;
clear all;
clc
Q = 1;
y = [10, 20, 30, 40, 50, 65, 80, 100, 200, 300, 400, 500, 650, 800, 1000];
sigmay = [27, 62, 115, 165, 210];
sigmaz = [14, 38, 105, 250 450];
u = 2;
for j = 1:y
for k = 1:sigmay
for l = 1:sigmaz
C = (Q./pi.*sigmay.*sigmaz.*u).*exp(-0.5.*(y./sigmay).^2); %% this is where the error is at%%
end
end
end
Warning: Colon operands must be real scalars. This warning will become an error in a future release.
Warning: Colon operands must be real scalars. This warning will become an error in a future release.
Warning: Colon operands must be real scalars. This warning will become an error in a future release.
Arrays have incompatible sizes for this operation.
plot(y,C)
title('Excercise 1')
xlabel('distance (m)')
ylabel('concentrations')
grid on
As you can see, I tried using the ./ or .* but, it didn't work.
  4 个评论
显示 2更早的评论
Stephen23
Stephen23 2024-4-10
"Any suggestions?"
Do you have a reference for the formula?
Jonathon Klepatzki
Jonathon Klepatzki 2024-4-10
Just this, C = (Q/pi*sigmaz*sigmay*u)*exp(-1/2(y/sigmay)^2)

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

the cyclist
the cyclist 2024-4-10
移动：Rik 2024-4-11
Ran in:
In this expression
Q/pi*sigmaz*sigmay*u
only pi is going to be in the denominator.
Since your expression bears some resemblance to the PDF of a normal distribution, I expect that at least sigmaz and sigmay are also supposed to be in the denominator. Probably u as well.
I expect that you actually need
Q/(pi*sigmaz*sigmay*u)
which would give
Q = 1;
y = [10, 20, 30, 40, 50, 65, 80, 100, 200, 300, 400, 500, 650, 800, 1000];
sigmay = [27, 62, 115, 165, 210];
sigmaz = [14, 38, 105, 250 450];
u = 2;
C = (Q./(pi.*sigmay.*sigmaz.*u)).*exp(-0.5.*(y(:)./sigmay).^2)
C = 15x5
1.0e-03 * 0.3931 0.0667 0.0131 0.0039 0.0017 0.3200 0.0641 0.0130 0.0038 0.0017 0.2271 0.0601 0.0127 0.0038 0.0017 0.1405 0.0549 0.0124 0.0037 0.0017 0.0758 0.0488 0.0120 0.0037 0.0016 0.0232 0.0390 0.0112 0.0036 0.0016 0.0052 0.0294 0.0103 0.0034 0.0016 0.0004 0.0184 0.0090 0.0032 0.0015 0.0000 0.0004 0.0029 0.0019 0.0011 0.0000 0.0000 0.0004 0.0007 0.0006
<mw-icon class=""></mw-icon>
<mw-icon class=""></mw-icon>
This looks to be in the ballpark of what you expected.

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Startup and Shutdown 的更多信息

标签

产品


版本

R2020b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by