Вопрос по MATLAB

[Krypto]

Помогите, я не халявщик и действительно пытаюсь понять.У меня возникла проблема. Задача была изобразить зависимость давления 1 моля Ван-дер-Ваальсовского газа от обьема (2D график) и от объема и температуры (3D график). Уравнение Ван-дер-Ваальса p=RT/(V-b)-a/V^2; параметры даны: R=0.082; a=3.67; b=0.013; T=298 (я опустил размерности). Ну вроде бы как всё. Задал значения переменных в интервалах и с шагом попытался написать само уравнение. Тут и возникли трудности. Как организовать поэлементные операции? Ведь перед значком операции нужно ставить точку. Но MATLABу это не нравится - ругается. В итоге я поубирал некоторые точки и ввёл такое: p=R*T./(V-0.013)-3.67./V.^2. Вообщем получилось у меня что-то невразумительное с отрицательным давлением. Как правильно ввести формулу? И главное - как ввести формулу для трёхмерного графика? 2D изотерма должна получиться такой (одна из них): http://www.chem.msu.su/rus/teaching/realgases/pic1.4.gif
Вместо этого получается какой-то бред.

[photon]

Очень трудно понять, что и как Вы задавали.

Вы не могли бы процитировать текст вашей программы, используя тег code. Типа

R=0.82;
a=3.67;
...

plot(V,p)

[Krypto]

>> R=0.082;
>> T=298;
>> a=3.67;
>> b=0.013;
>> V=[0:.1:10];
>> p=R*T./(V-0.013)-3.67./V.^2;
Warning: Divide by zero.
>> plot(V,p)

Какая-то изотерма непонятная. Чего-то снизу идёт.
Я не знаю как правильно ввести формулу для р (ну чтобы вычисления проводились с каждым элементом массива).

[photon]

Во-первых, исходя из физического смысла, $V>b$, что у Вас не выполняется.

Код Ваш работоспособный, только непонятно, зачем Вы определяли $a$, $b$, если потом внесли их в формулу в виде чисел.

[Krypto]

Сам не знаю почему. Меня тут все дергают. Отвлёкся.
У меня не получается постороить трёхмерный график (от объёма и температуры). Я пытаюсь, но он выдаёт мне какую-то плоскую линию.

[photon]

Krypto писал(а):

Сам не знаю почему. Меня тут все дергают. Отвлёкся.
У меня не получается постороить трёхмерный график (от объёма и температуры). Я пытаюсь, но он выдаёт мне какую-то плоскую линию.

посмотрите примеры с использованием команд meshgrid(), mesh(), surf()

[Krypto]

Большое всем спасибо! Я уже разобрался. Ещё вопрос.
Вот M-file("liqe.m"):

function y=liqe(k,p);
y=k(1)*k(2).*p./(1+k(1).*p)+(1-k(2))*k(3).*p./(1+k(3).*p);

А вот сам код:

>> teta=[.222 .224 .157 .128 .409 .482 .569 .669 .781 .806];
>> p=[2.9 3.3 7.7 12.3 23 44.9 93 180.6 360.1 722.2];
>> plot(p,teta,'o')
>> hold on
>> xlabel('Axis p')
>> ylabel('Axis teta')
>> k0=[1 1 1];
>> lsqcurvefit('liqe',k0,p,teta)
>> x=[0:800];
>> y=ans(1)*ans(2).*x./(1+ans(1).*x)+(1-ans(2))*ans(3).*x./(1+ans(3).*x);
>> plot(x,y)

Здесь всё очень сильно зависит от начального приближения. Что делать? Как сделать так, чтобы верхняя правая часть графика проходила ближе к точке. Как узнать эти начальные приближения?

[Kirill.Netreba]

Люди знающие, напишите, пожалуйста, как с математической точки зрения объяснить, что делает функция smooth? Или с физической...

[photon]

Kirill.Netreba в сообщении #144294 писал(а):

что делает функция smooth?

усредняет по нескольким точкам.

Как именно усредняет, зависит от алгоритма усреднения, который будет выбран:

Цитата:

Z = SMOOTH(Y) uses the moving average method with span 5 and X=1:length(Y).

Z = SMOOTH(Y,7) uses the moving average method with span 7 and X=1:length(Y).

Z = SMOOTH(Y,'sgolay') uses the Savitzky-Golay method with DEGREE=2, SPAN = 5, X = 1:length(Y).

Z = SMOOTH(X,Y,'lowess') uses the lowess method with SPAN=5.

Z = SMOOTH(X,Y,SPAN,'rloess') uses the robust loess method.

Z = SMOOTH(X,Y) where X is unevenly distributed uses the 'lowess' method with span 5.

Z = SMOOTH(X,Y,8,'sgolay') uses the Savitzky-Golay method with span 7 (8 is reduced by 1 to make it odd).

Z = SMOOTH(X,Y,0.3,'loess') uses the loess method where span is 30% of the data, i.e. span = ceil(0.3*length(Y)).