создание функции в matlab

@ЭСАУ · Регистрация: 07.11.2012

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

Создал функцию, сохранили её.

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
function [y,z] = may_diff(x0,xk,dx,e,dell)
 
x=x0:dx:xk;
 
if e==1
Er=randn(size(x))/dell; % рандомный массив ошибок размерностью size(x)
else Er=0;
end;
 
y=sin(x) + Er;
N=length(x); % length(x) - длинна массива х
for i=1:N-1
m(i)=(y(i+1)-y(i))/dx;
end
m(i+1)=m(i);
 
figure; plot( x , m ,'-r' ); % график функции
figure; plot( x , y, '-k' ); % график продифференцированной функции
 
end

теперь пробую вызвать: may_diff(0,10,0.01,1,10)
строятся графики, всё нормально!

как сделать чтобы значения x,y,m,e - записывались в рабочей области ? чтобы в дальнейшем можно было к ним обратиться или при желании проводить какой-либо анализ.

P.S. функцию пишу первый, может второй раз.

Зосима · 21.05.2014, 23:59

тут разве что сделать эти переменные глобальными, тогда они появятся в рабочем пространстве

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
function [y,z] = may_diff(x0,xk,dx,e,dell)
global x y m e 
x=x0:dx:xk;
 
if e==1
Er=randn(size(x))/dell; % рандомный массив ошибок размерностью size(x)
else Er=0;
end;
 
y=sin(x) + Er;
N=length(x); % length(x) - длинна массива х
for i=1:N-1
m(i)=(y(i+1)-y(i))/dx;
end
m(i+1)=m(i);
 
figure; plot( x , m ,'-r' ); % график функции
figure; plot( x , y, '-k' ); % график продифференцированной функции
 
end

а потом вызываем эту красоту

Matlab M
1
2
>> global x y m e 
>> may_diff(0,10,0.01,1,10)

@ЭСАУ · 22.05.2014, 00:26 **[ТС]**

Можно ли, сложно ли сделать так чтобы при вызове функции у пользователя запрашивались данные эти данные.
к примеру: введи начальное значение (хО)

ждет ввода), введите шаг (dx): (ждет ввода)

что то подобное делалось в scilabe.

Зосима · 22.05.2014, 09:28

ты хочешь, чтобы данные вводились внутри функции?

вообщем это делается с помощью команды input (описание):

Matlab M
1
2
x0 = input('Введите начальное значение х0: ');
dx = input('Введите шаг dx: ');

@ЭСАУ · 22.05.2014, 13:42 **[ТС]**

Большое спасибо! всё заработало как хотелось!

Добавлено через 2 часа 33 минуты

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
function [ ] = may_diff2( )
clc;
global x y m er 
 
disp('Метод цифрового дифференцирования сигналов,');
disp('основанный на использовании конечных приращений');
fun = input('Введите функцию (пример sin(x)): ');
x0 = input('Введите начальное значение, х0: ');
xk = input('Введите конечно значение, xk: ');
dx = input('Введите шаг дискретизации, dx: ');
e = input('Для добавления ошибки введите 1, иначе 0: ');
 
x=x0:dx:xk;
 
y1 = fun;
if e == 1
disp('к основному сигналу будет добавлена ошибка');    
er = randn(size(x)); % рандомный массив ошибок размерностью size(x)
y = plus(y1,er);
elseif e == 0 
disp('сигнал без ошибок')
y = y1;
else
disp('некоректно')
end;  
 
N = length(x); % length(x) - длинна массива х
for i = 1:N-1
m(i) = (y(i+1)-y(i))/dx;
end
m(i+1) = m(i);
 
figure; plot( x , y, '-k' ); % график продифференцированной функции
figure; plot( x , m ,'-r' ); % график функции 
 
end

Часто выдает такие ошибки:

??? Attempted to access y(1002); index out of bounds because numel(y)=1001.

Error in ==> may_diff2 at 30
m(i) = (y(i+1)-y(i))/dx;

??? Error using ==> plus
Matrix dimensions must agree.

Error in ==> may_diff2 at 19
y = plus(y1,er);

??? Error using ==> plot
Vectors must be the same lengths.

Error in ==> may_diff2 at 33
figure; plot( x , y, '-k' ); % график продифференцированной функции

вроде бы размерности всех элементов заданы корректно, и никаких проблем быть не должно ( на мой взгляд)

Зосима · 22.05.2014, 13:56

у тебя вообще как-то странно все: y1 - текстовая строка, как вообще в y может появиться массив чисел из 1002 элементов...

в 28-й строке напиши for i = 1:length(y)-1 так точно с у проблем не будет

кстать, а почему ты не используешь команду diff, которая делает ( y(i+1)-y(i) ) вместо цикла?
и если уж ты все данные вводишь, то можно и функцию убрать

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
clear, clc;
 
disp('Метод цифрового дифференцирования сигналов,');
disp('основанный на использовании конечных приращений');
fun = input('Введите функцию (пример sin(x)): ');
x0 = input('Введите начальное значение, х0: ');
xk = input('Введите конечно значение, xk: ');
dx = input('Введите шаг дискретизации, dx: ');
e = input('Для добавления ошибки введите 1, иначе 0: ');
 
x=x0:dx:xk;
 
y1 = eval(fun,x); % считаем значения ф-ции
if e == 1
  disp('к основному сигналу будет добавлена ошибка');    
  er = randn(size(x)); % рандомный массив ошибок размерностью size(x)
  y = y1+er;
elseif e == 0 
  disp('сигнал без ошибок')
  y = y1;
else
  disp('некоректно')
end;  
 
m = diff(y)/dx;
m(end+1) = m(end);
 
plot( x , y, '-k' ); % график продифференцированной функции
figure; plot( x , m ,'-r' ); % график производной

@ЭСАУ · 22.05.2014, 16:22 **[ТС]**

при использование функции eval выдает такую ошибку.
??? Undefined function or method 'eval' for input arguments of type 'double'.

если задать функцию заранее, то всё работает.

Но, причем если использовать:

Matlab M
1
2
m = diff(y)/dx;
m(end+1) = m(end);

Если использовать:

Matlab M
1
2
3
4
for i = 1:length(y)-1
m(i) = (y(i+1)-y(i))/dx;
end
m(i+1) = m(i);

временами выдает ошибку о несоответствии размерностей

Зосима · 22.05.2014, 16:40

покажи пример входных данных, при которых у тебя вылазит ошибка

чтобы eval не давало ошибки, ввод функции надо сделать в строковом формате:

Matlab M
fun = input('Введите функцию (пример sin(x)): ', 's');

@ЭСАУ · 22.05.2014, 16:52 **[ТС]**

Введите функцию (пример sin(x)): sin(x)
Введите начальное значение, х0: 0
Введите конечно значение, xk: 13
Введите шаг дискретизации, dx: 0.05
Для добавления ошибки введите 1, иначе 0: 1
к основному сигналу будет добавлена ошибка
??? Error using ==> plot
Vectors must be the same lengths.

Error in ==> may_diff2 at 37
figure; plot( x , m ,'-r' ); % график производной

Код программы:

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
function [ ] = may_diff2( )
clc;
global x y m er
 
disp('Метод цифрового дифференцирования сигналов,');
disp('основанный на использовании конечных приращений');
fun = input('Введите функцию (пример sin(x)): ');
x0 = input('Введите начальное значение, х0: ');
xk = input('Введите конечно значение, xk: ');
dx = input('Введите шаг дискретизации, dx: ');
e = input('Для добавления ошибки введите 1, иначе 0: ');
 
x=x0:dx:xk;
 
y1 = sin(x);
 
if e == 1
disp('к основному сигналу будет добавлена ошибка');    
er = randn(size(x)); % рандомный массив ошибок размерностью size(x)
y = y1 + er;
elseif e == 0 
disp('сигнал без ошибок')
y = y1;
else
disp('некорректно')
end;  
 
for i = 1:length(y)-1
m(i) = (y(i+1)-y(i))/dx;
end
m(i+1) = m(i);
 
plot( x , y, '-k' ); % график продифференцированной функции
figure; plot( x , m ,'-r' ); % график производной
 
end

если заменить 28-31 строки на:

Matlab M
1
2
m = diff(y)/dx;
m(end+1) = m(end);

тогда при этих же входных данных всё работает!
использую Matlab 7.6.0 (R2008a)

Зосима · 22.05.2014, 17:14

ага! вроде разобрался

смотри, m у тебя глобальная. При запуске с одним шагом будет одно кол-во элементов, например 1000, в другой раз при запуске m не стирается, а так и остается длинны 1000! при другом запуске, когда поменяли данные и точек стало меньше, например 100, в цикле пересчитаются первые 100 элементов, но общая длинна m так и будет 1000, при этом длинна х будет 100 - вот и несоответствие!

Команды diff и zeros полностью перезаписывают содержимое переменной m, в то время, как цикл изменяет лишь некоторые значения!
чтобы это "вылечить", я добавил инициализацию массива m (строка 27)

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
clc;
global x y m er
 
disp('Метод цифрового дифференцирования сигналов,');
disp('основанный на использовании конечных приращений');
fun = input('Введите функцию (пример sin(x)): ', 's');
x0 = input('Введите начальное значение, х0: ');
xk = input('Введите конечно значение, xk: ');
dx = input('Введите шаг дискретизации, dx: ');
e = input('Для добавления ошибки введите 1, иначе 0: ');
 
x=x0:dx:xk;
 
y1 = eval(fun,x);
 
if e == 1
    disp('к основному сигналу будет добавлена ошибка');
    er = randn(size(x)); % рандомный массив ошибок размерностью size(x)
    y = y1 + er;
elseif e == 0
    disp('сигнал без ошибок')
    y = y1;
else
    disp('некорректно')
end;
 
m = zeros(size(y));
for i = 1:length(y)-1
    m(i) = (y(i+1)-y(i))/dx;
end
m(i+1) = m(i);
 
plot( x , y, '-k' ); % график продифференцированной функции
figure; plot( x , m ,'-r' ); % график производной

@ЭСАУ · 22.05.2014, 17:27 **[ТС]**

перед вызовом функции каждый раз очищал рабочее пространство, командой

Matlab M
1
clear

возможно в какой то момент упустил.

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Desktop (MinGW): Рисуем цветные прямоугольники с помощью рисовальщика SDL3 на Си и C++ 8Observer8 17.03.2026 Содержание блога Финальные проекты на Си и на C++: finish-rectangles-sdl3-c. zip finish-rectangles-sdl3-cpp. zip	Символические и жёсткие ссылки в Linux. algri14 15.03.2026 Существует два типа ссылок — символические и жёсткие. Ссылка в Linux — это запись в каталоге, которая может указывать либо на inode «файла-ИСТОЧНИКА», тогда это будет «жёсткая ссылка» (hard link),. . .	[Owen Logic] Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ФедосеевПавел 14.03.2026 Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ВВЕДЕНИЕ Выполняя задание на управление насосной группой заполнения резервуара,. . .	делаю науч статью по влиянию грибов на сукцессию anaschu 13.03.2026 прикрепляю статью
SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:	Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .	Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .	Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .

@ЭСАУ 1 / 1 / 0 Регистрация: 07.11.2012 Сообщений: 42
	22.05.2014, 00:26 [ТС]
	Можно ли, сложно ли сделать так чтобы при вызове функции у пользователя запрашивались данные эти данные. к примеру: введи начальное значение (хО)ждет ввода), введите шаг (dx): (ждет ввода) что то подобное делалось в scilabe. 0

Зосима 5245 / 3573 / 379 Регистрация: 02.04.2012 Сообщений: 6,477 Записей в блоге: 18
	22.05.2014, 16:40
	покажи пример входных данных, при которых у тебя вылазит ошибка чтобы eval не давало ошибки, ввод функции надо сделать в строковом формате: Matlab M fun = input('Введите функцию (пример sin(x)): ', 's'); 0