Ситуация, когда прибыль превышает, издержки
%Задаем
начальные условия
t=0:.01:.5;
a1=cos(3.*t+5);
a2=t.^2+3.*t+4;
N0=5;
p=5;
%Прибыль от одного акта продажи
s=13;
%Издержки на рекламные действия
%%%%%%%%%%%%%%%%%%%%%%
[t,y]=
ode45(@fun1,[0 3],0,[],N0);
P=p*N0*(1/3)*sin(3*t);
S=s*(1/3)*sin(3*t);
plot(t,P,'r-');
xlabel('Time');
ylabel('Benefit');
hold on;
axis([0 .5
0 10]);
grid;
pause;
plot(t,S,'b-');
ylabel('Benefit,
Expenses');
legend('Benefit','Expenses',2);
hold off;
%%%%%%%%%%%%%%%%%%%%%%%%%%
//fun1.m
function
dN=fun1(t,N,N0);
%
a1=cos(3.*t+5);
a2=t.^2+3.*t+4;
dN=(a1+a2*N)*(N0-N);
//
Приложение
2:
Ситуация,
когда рекламная кампания убыточная.
…
s=40;
…
Приложение
3:
Ситуация
с постоянными коэффициентами 1
и 2:
a1=.1;
a2=.05;
p=5;
s=13;
N0=200;
t=4;
[t,N]=ode45(@fun2,[0
t],0,[],a1,a2,N0);
plot(t,N,'r-')
grid;
xlabel('time')
ylabel(
'Quantity of clients')
pause;
//
Fun2.m
function
dN=rec1(t,N,a1,a2,N0)
dN=(a1+a2*N)*(N0-N);
//
n0=(a1/a2)+N0;
k=a1/a2;
%Constants
[t
m]=ode45(@fun3,[0 6],k,[],a1,a2,n0);
n=a2*m.*(n0-m);
plot(t,n,'r-');
xlabel('Time');
ylabel('Quantity
of clients');
legend('For
a long period of time');
grid;
//
Fun3.m
function
du=rec2(t,u,a1,a2,n0)
du=a2*u*(n0-u);
//
Приложение
4:
Конкурентная
борьба двух компаний за покупателя.
a1=0.68;
a2=0.28;
b1=0.98;
b2=0.53;
k1=
cos(3.*p+5)+p.^2+3.*p;
k2=.8*sin(6.*p+4)+exp(-p);
N0=20;
N1=13;
N2=N0-N1;
[p,y]=
ode45(@fun2,[0 10],[N1 N2],[],N0);
plot(p,y(:,1),p,y(:,2),'g');
grid;
legend('Company_1','Company_2');
xlabel('Price');
ylabel('Company
1,2');
pause;
plot(y(:,1),y(:,2),'r');
xlabel('Company
1');
ylabel('Company
2');
grid;
//
Fun4.m
function
fun=fun1(p,y,a1,a2,b1,b2,k1,k2,N0)
%
fun=zeros(2,1);
a1=0.68;
a2=0.28;
b1=0.98;
b2=0.53;
N0=20;
%
fun(1)=(a1+a2.*y(1))*(N0-y(1))*(cos(3.*p+5)+p.^2+3.16*p);
fun(2)=(b1+b2*y(2))*(N0-y(1)-y(1))*(.8*sin(6.*p+4)+exp(-p));
//
16