Организовать расчет плотности в цикле

24.12.2012, 21:41. **Показов** 899. **Ответов** 1

Author24 — интернет-сервис помощи студентам

День добрый! Нужна помощь в создании цикла, где необходимо найти плотность воды.
Есть массив данных (я перевел сплошной текст в матрицу, которая читаема)
В массиве 3 столбца: глубина, температура, соленость(по 200 строк)

И я наработал правильно(!) уравнение состояния, по которому находится эта плотность.
Вот оно:
(t и s - задаются; t1-t3 - постоянные )

Matlab M

t1=3.98
t2=283
t3=67.26
a1=t*(4.7867-0.098185*t+0.0010843*(t^2))*(10^(-3));
b1=t*(18.030-0.8164*t+0.01667*(t^2))*(10^(-6));
cl=s/1.80655;
plotnost0=-0.069+1.4708*cl-0.001570*(cl^2)+0.0000398*(cl^3);
dt=-(((t-t1)^2)/503.570)*((t+t2)/(t+t3));
P=dt+(plotnost0+0.1324)*(1-a1+b1*(plotnost0-0.1324));

P - окончательная плотность для задающейся t(температура) и s(соленость)

Вопрос вот в чем: как должен выглядеть цикл, который будет подставлять t и s в это уравнение, дабы получать искомую P?
Должен получиться массив данных(1 столбец в 200 значений плотности(P))

Заранее благодарю за советы!

Зосима · 25.12.2012, 10:24

Дружочек, если бы ты еще выложил данные (в файле например, то было бы вообще супер!)

Набросал на скорую руку и без циклов:

Matlab M

clear, clc
 
t = linspace(1,40,200); % 200 значений температуры
s = linspace(0.02,0.1,200); % 200 значений солености
% или M = load('data.txt');
% t = M(:,2);
% s = M(:,3);
 
t1=3.98; t2=283; t3=67.26;
a1 = t.*(4.7867-0.098185*t+0.0010843*t.^2)*1e-3;
b1 = t.*(18.030-0.8164*t+0.01667*t.^2)*1e-6;
cl = s/1.80655;
plotnost0 = -0.069+1.4708*cl-0.001570*cl.^2 + 0.0000398*cl.^3;
dt = -(((t-t1).^2)/503.570).*((t+t2)./(t+t3));
P = dt + (plotnost0+0.1324).*(1-a1+b1.*(plotnost0 - 0.1324));
disp([t', s', P'])

А если так хочешь в циклами, то

Matlab M

clear, clc
 
t = linspace(1,40,200); % 200 значений температуры
s = linspace(0.02,0.1,200); % 200 значений солености
% или M = load('data.txt');
% t = M(:,2);
% s = M(:,3);
 
t1=3.98; t2=283; t3=67.26;
 
for i = 1:length(t) % важно чтобы t и s были одной длинны!
    a1 = t(i)*(4.7867-0.098185*t(i)+0.0010843*t(i)^2)*1e-3;
    b1 = t(i)*(18.030-0.8164*t(i)+0.01667*t(i)^2)*1e-6;
    cl = s(i)/1.80655;
    plotnost0 = -0.069+1.4708*cl-0.001570*cl^2 + 0.0000398*cl^3;
    dt = -(((t(i)-t1)^2)/503.570)*((t(i)+t2)/(t(i)+t3));
    P(i) = dt + (plotnost0+0.1324).*(1-a1+b1*(plotnost0 - 0.1324));
end
disp([t', s', P'])

Значения t и s вводил от балды, поэтому и результат соответствующий.

Но по большому счету P(t,s) - это поверхность, т.е. матрица, а не массив! Но мы пока не будем влазить в дебри матричной магии

	25.12.2012, 10:24