Не выводится график

@Anna111 · Регистрация: 17.12.2015

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

Добрый вечер!
Помогите пожалуйста решить проблему, график не выводится( Я что-то не так делаю?

Matlab M

clear, clc
h = 0.1; % шаг
a = -5; % нижний предел
b = 5; % верхний предел
p = 2.5;
eps1 = 2.71;
eps2 = 0.8;
q = 1.5;
w = 0.1;
% функция системы:
dsys = @(x,y)[1/q((1-q/p*ln(1+p*x)-(w-1)*(y^3-1)/eps2)*exp(x/(1+p*x))-eps1*y*x*(1+p*x))];
x = a:h:b;
 
y = zeros(size(x));
for j=1:length(x)-1
    k1=h*dsys(x(j), y(j));
    k2=h*dsys(X(j)+h/2,Y(j)+k1(1)/2);
    k3=h*dsys(X(j)+h/2,Y(j)+k2(1)/2);
    k4=h*dsys(X(j)+h,Y(j)+k3(1));
    dU = (k1+2*k2+2*k3+k4)/6;
    y(j+1)=y(j) + dU;
end 
 
plot(x,y);
legend('y(x)');

@Krasme · 07.03.2016, 21:31

скобки исправить и ln заменить на log
а также перепроверить формулу, т.к. после первого q неизвестно, какой должен быть знак, я поставила умножение (наобум) ..

Matlab M

1	dsys = @(x,y)(1/q((1-q/plog(1+px)-(w-1)(y^3-1)/eps2)exp(x/(1+px))-eps1yx(1+px)));

а также надо учитывать, что матлаб различает регистры, т.е. x и X - разные переменные

@Anna111 · 08.03.2016, 10:19 **[ТС]**

Помогите пожалуйста исправить ошибки и как сделать так, чтобы 2 графика в одном окне графическом отображались?

Matlab M

clear, clc
h = 0.1; % шаг
a = -1; % нижний предел
b = 1.5; % верхний предел
p = 2.5;%beta
eps1 = 2.71;
eps2 = 0.8;
q = 1.5;%gama
w = 0.1;%psi
% функция системы:
dsys = @(x,y)(((1-(w-1)*(y^3-1)/eps2)*exp(x)-eps1*y*x));
x = a:h:b;
 
y = zeros(size(x));
for j=1:length(x)-1
    k1=h*dsys(x(j), y(j));
    k2=h*dsys(x(j)+h/2,y(j)+k1(1)/2);
    k3=h*dsys(x(j)+h/2,y(j)+k2(1)/2);
    k4=h*dsys(x(j)+h,y(j)+k3(1));
    dU = (k1+2*k2+2*k3+k4)/6;
    y(j+1)=y(j) + dU;
end 
 
function z = myfun(X, Y)
p = 2.5;
eps1 = 2.71;
eps2 = 0.8;
q = 1.5;
w = 0.1;
z = (1 + (w - 1.0) * (1 - Y^3) / eps2) * exp(X / (1 + p * X))...
    - eps1 * Y * X * (1 + p * X);
end
 
ezplot(@(x,y)myfun(x,y), [-1,1.5,-1,1]);
plot(x,y);
legend('y(x)');

@Dmitriu_L · 09.03.2016, 13:23

Если все в одном файле тогда так:

Matlab M

function main
clear, clc
h = 0.1; % шаг
a = -1; % нижний предел
b = 1.5; % верхний предел
p = 2.5;%beta
eps1 = 2.71;
eps2 = 0.8;
q = 1.5;%gama
w = 0.1;%psi
% функция системы:
dsys = @(x,y)(((1-(w-1)*(y^3-1)/eps2)*exp(x)-eps1*y*x));
x = a:h:b;
 
y = zeros(size(x));
for j=1:length(x)-1
    k1=h*dsys(x(j), y(j));
    k2=h*dsys(x(j)+h/2,y(j)+k1(1)/2);
    k3=h*dsys(x(j)+h/2,y(j)+k2(1)/2);
    k4=h*dsys(x(j)+h,y(j)+k3(1));
    dU = (k1+2*k2+2*k3+k4)/6;
    y(j+1)=y(j) + dU;
end 
figure
subplot(2,1,1);
ezplot(@(x,y)myfun(x,y), [-1,1.5,-1,1]);
subplot(2,1,2);
plot(x,y);
legend('y(x)');
 
function z = myfun(X,Y)
p = 2.5;
eps1 = 2.71;
eps2 = 0.8;
q = 1.5;
w = 0.1;
z = (1 + (w - 1.0) * (1 - Y.^3) / eps2).* exp(X ./ (1 + p * X))...
    - eps1*Y.* X.* (1 + p * X);

Ознакомьтесь с темой во избежание дальнейших ошибок при работе с векторами СЕКРЕТНО! Зловещая точка претыкания

@Anna111 · 15.03.2016, 19:16 **[ТС]**

Спасибо, но не подскажете как в решении myfun увеличить точность, чтобы не было разрыва.

@SSC · 16.03.2016, 13:19

при х=-0.4 знаменатель равен 0, а exp() от + и - бесконечности сильно отличаются, поэтому разрыв

@Anna111 0 / 0 / 0 Регистрация: 17.12.2015 Сообщений: 22
	15.03.2016, 19:16 [ТС]	5
	Спасибо, но не подскажете как в решении myfun увеличить точность, чтобы не было разрыва. 0

@SSC $Эксперт по математике/физике$ 3390 / 1913 / 571 Регистрация: 09.04.2015 Сообщений: 5,365
	16.03.2016, 13:19	6
	при х=-0.4 знаменатель равен 0, а exp() от + и - бесконечности сильно отличаются, поэтому разрыв 1