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

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

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

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

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
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/p*log(1+p*x)-(w-1)*(y^3-1)/eps2)*exp(x/(1+p*x))-eps1*y*x*(1+p*x)));

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

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

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

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
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
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
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() от + и - бесконечности сильно отличаются, поэтому разрыв

Новые блоги и статьи Все статьи Все блоги /
Символьное дифференцирование igorrr37 13.02.2026 / * Логарифм записывается как: (x-2)log(x^2+2) - означает логарифм (x^2+2) по основанию (x-2). Унарный минус обозначается как ! */ #include <iostream> #include <stack> #include <cctype>. . .	Камера Toupcam IUA500KMA Eddy_Em 12.02.2026 Т. к. у всяких "хикроботов" слишком уж мелкий пиксель, для подсмотра в ESPriF они вообще плохо годятся: уже 14 величину можно рассмотреть еле-еле лишь на экспозициях под 3 секунды (а то и больше),. . .	И ясному Солнцу zbw 12.02.2026 И ясному Солнцу, и светлой Луне. В мире покоя нет и люди не могут жить в тишине. А жить им немного лет.	«Знание-Сила» zbw 12.02.2026 «Знание-Сила» «Время-Деньги» «Деньги -Пуля»
SDL3 для Web (WebAssembly): Подключение Box2D v3, физика и отрисовка коллайдеров 8Observer8 12.02.2026 Содержание блога Box2D - это библиотека для 2D физики для анимаций и игр. С её помощью можно определять были ли коллизии между конкретными объектами и вызывать обработчики событий столкновения. . . .	SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL_LoadPNG (без SDL3_image) 8Observer8 11.02.2026 Содержание блога Библиотека SDL3 содержит встроенные инструменты для базовой работы с изображениями - без использования библиотеки SDL3_image. Пошагово создадим проект для загрузки изображения. . .	SDL3 для Web (WebAssembly): Загрузка PNG с прозрачным фоном с помощью SDL3_image 8Observer8 10.02.2026 Содержание блога Библиотека SDL3_image содержит инструменты для расширенной работы с изображениями. Пошагово создадим проект для загрузки изображения формата PNG с альфа-каналом (с прозрачным. . .	Установка Qt-версии Lazarus IDE в Debian Trixie Xfce volvo 10.02.2026 В общем, достали меня глюки IDE Лазаруса, собранной с использованием набора виджетов Gtk2 (конкретно: если набирать текст в редакторе и вызвать подсказку через Ctrl+Space, то после закрытия окошка. . .

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

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