[Lorem Solis]

Добрый день. На практике пишем программы в матлабе. И не понятно, что не так с моими прогами. Буду очень признательна за помощь.
1. y = –3•exp(0,5•x), x меняется в диапазоне [–10, 10], шум эксперимента задается случайной величиной, меняющейся в диапазоне [–20,20].
Что я написала:
clear all;
figure
%x=rand(1,5)
n = 5;
x= -10+20*rand(1,n);
x = sort(x);
y=-3*exp(-0.5*x);
plot(x,y,'*r')
hold on
% disp( y );
z = -0.003+ 0.006*rand(1,n);

%z=-0.01+0.02*rand(1,n)
y=y+z;
plot(x,y,'b')
% hold on
u=log(-y);
Mx=0;
My=0;
Mxx=0;
Mxy=0;
%for i=1:n
% Mx=Mx+x(i);
% My=My+u(i);
% Mxx=Mxx+((x(i)).^2);
%Mxy=Mxy+(x(i)*u(i));
%end
Mx=sum(x);
My=sum(y);
Mxx=sum(x.*x);
Myy=sum(y.*y);

Mx=Mx/n;
My=My/n;
Mxx=Mxx/n;
Mxy=Mxy/n;

a=exp((Mxy-(My*Mx))/(Mxx-(Mx*Mx)));
b=-exp(((Mxx*My)-(Mx*Mxy))/((Mxx-(Mx*Mx))));

y=a*x+b;
plot(x,y,'g')
% hold on
% konec
% close all
Что с ней не выходит: график должен накладываться на график, а у меня этого не происходит.

2.5.Паромная переправа
Механизированный паром должен перевозить 1000 тонн груза в течение суток через реку шириной 3 км. Стоимость парома без двигателя пропорциональна его грузоподъемности L и равна 2000*L, а стоимость меха-нического двигателя пропорциональна произведению грузоподъемности на куб скорости его движения V и равна 5000*L*V3. Рассчитать оптимальные грузоподъемность двигателя парома обеспечивающего перевоз груза с минимальной стоимостью перевозки
Как я получала функцию, которую надо минимизировать:
1000 / x туда и еще (1000 / x) - 1 обратно
заметим, что количество раз * 3 это расстояние в километрах
обозначим это за S
теперь скорость y
S/y = 24
y = S / 24
2000x + 5000x * y^3
в итоге функция, которую нужно минимизировать, получается такой:
2000*x+5000*x*((1000/x-1)^3)*(3^3)/(2^9);


но она монотонна и поэтому прога выдает всегда верхнюю границу отрезка, на котором я ищу минимум (ищу методом золотого сечения).

сама прога

clear all
%f=500000/y+1250000*y^2;
%RetVal=2000*x+5000*x*((1000/x-1)^3)*(3^3)/(2^9);

a=0
b=5
k=0.618
e=0.001
while (b-a>e)
x1=a+(1-k)*(b-a);
x2=a+k*(b-a);
if (f(x1)>f(x2))
a=x1;
end
if (f(x1)<f(x2))
b=x2;
end
end
gruzopod=f(x1)
close all

Очень-очень выручите. Хотя бы со второй. На мой взгляд, она легче.
Очень прошу. Заранее благодарю.

@темы: Вопрос, C#

Comments

[All_ex]

1000 / x туда и еще (1000 / x) - 1 обратно - как-то смущает это утверждение ... вроде паром к началу следующих суток должен оказаться на стартовой стороне...

И вид функции у меня другой получился...
`T = 6/V` - время на один рейс...
`24/T = 4*V` - число рейсов за сутки... оно же равно `1000/L`, откуда, `V = 250/L` ...
Итого, функция стоимости имеет вид `2000*L + 5000*L*(250/L)^3`... и ответ легко руками вычисляется ...

По-моему, так...

[All_ex]

Я могу ошибаться... но в №1 смущает следующее
Mx=sum(x);
My=sum(y);
Mxx=sum(x.*x);
Myy=sum(y.*y); - тут зачем-то вычисляется сумма квадратов игреков, которая не используется далее...

Mx=Mx/n;
My=My/n;
Mxx=Mxx/n;
Mxy=Mxy/n; - а где выше вычислена сумма произведений?...

[Lorem Solis]

All_ex у меня теперь тоже такое уравнение :3 ага, все вручную делается. не знаю, зачем нам такие задания дают

1ую тоже поправила. но дело, похоже, было не в этом. кажется, даже работает :3

спасибо большое :3

[All_ex]

Lorem Solis, не знаю, зачем нам такие задания дают - ну, простое учебное задание... и проверить свои вычисления можно...
похоже, было не в этом. - ... ...

в общем, welcome...