Форум программистов, компьютерный форум, киберфорум
Наши страницы
Matlab
Войти
Регистрация
Восстановить пароль
 
Kimmi
0 / 0 / 0
Регистрация: 25.11.2017
Сообщений: 3
1

Построение семейства резонансных кривых

04.12.2017, 11:44. Просмотров 416. Ответов 2
Метки нет (Все метки)

Задание выглядит так: построить семейство резонансных кривых нелинейного осциллятора http://www.cyberforum.ru/cgi-bin/latex.cgi?\ddot{x} + \alpha \dot{x} + x - \mu  {x}^{3} = sin(2\pi ft) при разных α и μ.

Написал функцию:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
function y = ResKf(t,x)
 
global a;
global n;
global w;
 
%x(1) = x(t)
%x(2) = dx/dt
%y(1,1) = dy/dt
%y(2,1) = d2y/dt2
 
y(1,1)= x(2);
y(2,1)= sin(w*t)-a*x(2)-x(1)+n*((x(1))^3);
 
end
Скрипт:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
clear;
 
global a;
global n;
global w;
n = 0;
tk = [0 40];
y0 = [0.1 0.01];
for a = 0.3:0.1:0.8
    i=1;
    for n = 0.00:0.01:0.04
        for w = 0.3:0.01:2
[t,u] = ode45('ResKf',tk,y0);
u1 = u(1:end,1);
A(i)=(max(u1(t>10))-min(u1(t>10)))/2;
i=i+1;
        end
    end
plot(0.3:0.01:2,A) 
hold on
end
и выдает ошибку:

Matlab M
1
2
3
4
5
Error using plot
Vectors must be the same length.
 
Error in ResK (line 19)
plot(0.3:0.01:4,A)
При этом, если n (\mu) я меняю вручную, то все строится замечательно.

Прошу помочь, как все построить на одном графике и "без рук"?)
0
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
04.12.2017, 11:44
Ответы с готовыми решениями:

Построение семейства кривых
Доброго времени суток! Прошу при возможности помочь мне разобраться в том, как построить семейство...

Суммирование семейства кривых в Матлаб
Пусть требуется вывести семейство кривых для x€ , которое задано функцией, зависящей от параметра ...

Построение семейства графиков
Здравствуйте, форумчане! Пришел к Вам за помощью. Знаю, что можно и самому, но чего-то не...

Построение фокусировочных кривых
Доброго времени суток:) Занимаюсь автофокусировкой опт. микроскопа. Мне надо сравнить некоторые...

Построение кривых второго порядка.Гипербола и парабола
Доброго времени суток! Помогите, пожалуйста, построить гиперболу и параболу. Я не могу построить т....

2
SSC
Эксперт по математике/физике
3138 / 1534 / 465
Регистрация: 09.04.2015
Сообщений: 4,275
04.12.2017, 16:09 2
Цитата Сообщение от Kimmi Посмотреть сообщение
все строится замечательно.
Не знаю что у Вас строится замечательно, но в функции plot 1 вектор (по оси х ) имеет 171 элемент, а вектор A имеет 855 элементов. В А Вы по сути сохранили последовательно точки для 5 графиков, естественно ничего строится не будет.
Может выходом является построение двумерного массива A, и потом построение последовательно 5 графиков для различных n
1
Kimmi
0 / 0 / 0
Регистрация: 25.11.2017
Сообщений: 3
05.12.2017, 16:20  [ТС] 3
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
clear;
 
global a;
global n;
global w;
n = 0;
tk = [0 40];
y0 = [0.1 0.01];
 
i=0;
 
for a = 0.3:0.1:0.8  
    for n = 0.00:0.01:0.02
        i = i+1;
        j=1;
        for w = 0.3:0.01:2
            
            [t,u] = ode45('ResKf',tk,y0);
            u1 = u(1:end,1);
            A(i,j)=(max(u1(t>10))-min(u1(t>10)))/2;
            
            j=j+1;
        end
        plot(0.3:0.01:2,A(i,:))
        hold on
        title('Cемейство резонансных кривых нелинейного осциллятора'); 
    end
end
0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
05.12.2017, 16:20

Построение нескольких кривых одной функцией plot
Играюсь с одной анимацией в которой из-за большого количества обьектов начинает подтормаживать. В...

Построение кривых второго порядка (окружность, эллипс, парабола, гипербола)
общее задание: нужно построить гиперболу с заданными "параметрами" a и b, а также оси координат,...

Построение семейства резонансных характеристик фильтра в 3D
Нужно построить 3D график семейства резонансных характеристик фильтра, используя эти формулы из...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
3
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Рейтинг@Mail.ru