Необходимо рассчитать интеграл и построить график

@Валентина17 · Регистрация: 02.06.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
26
%s=0.07;
v0=3;
x=0.1;
vz=(v0*0.265)/x;
r=1;
z=[0.001:0.1:10];
c0=0.082;
w=r/c0;
v=vz*exp(-0.5*(w*z).^2);
syms z
y=v;
v1=diff(y);
A2=2;
A3=1;
j=A3*v1;
f0=40000;
c=331;
f=(v/c)*2*f0;
t=[0.01:0.1:10];
%t=10;
a=0;
b=10;
tol=1.e-6;
S=A2.*v1.*exp(-j).*cos(6.28.*(f0+f).*t);
I=guad(S,a,b,tol);
plot (S,z);

Выдает такую ошибку ??? Error using ==> times
Matrix dimensions must agree.

Error in ==> diplom at 24
S=A2.*v1.*exp(-j).*cos(6.28.*(f0+f).*t);

Считаю интеграл по упрощенной формуле.

@Centurio · 28.03.2016, 01:04

Сообщение от Валентина17

Выдает такую ошибку ??? Error using ==> times
Matrix dimensions must agree.
Error in ==> diplom at 24
S=A2.*v1.*exp(-j).*cos(6.28.*(f0+f).*t);

В сообщении говорится, что размеры матриц не совпадают. У вас вектор v1 состоит из 99 элементов, а вектор t - из 100 элементов, вот и ошибку выдаёт.

@SSC · 28.03.2016, 07:46

Функция diff (стр.12) возвращает массив на единицу меньше чем исходный для диффиренцирования

@Валентина17 · 28.03.2016, 12:44 **[ТС]**

То есть можно увеличить z на одну единицу и тогда будет работать?

@SSC · 28.03.2016, 13:38

Сообщение от Валентина17

и тогда будет работать

Работать то может и будет, а будет ли считать то что нужно?
Нам Ваша задача неизвестна.
Запустив Ваш код, мы выявили ошибку и причину ее вызывающую (что и было в вопросе).
Зачем нужно дифферинцирование для

Сообщение от Валентина17

Считаю интеграл по упрощенной формуле

мне совершенно непонятно.
Как будет вести себя оператор

Matlab M
1
plot (S,z);

при выше определенной

Matlab M
1
syms z

мне неведомо, скорее всего даст ошибку.
Про функцию

Matlab M
1
I=guad(S,a,b,tol);

я тоже ничего не знаю и help ничего не выдает (но может у меня нет нужного приложения).
А может это функция quad? Но тогда первым параметром должна быть функция, а не массив чисел S.
Про остальное мне просто неведомо ...

Не по теме:

А в остальном, прекрасная маркиза, ...........

@Валентина17 · 28.03.2016, 14:09 **[ТС]**

Есть два графика и исходные данные необходимо построить S, я применяла diff для того чтобы вычислить производную от скорости v, и guad для нахождения интеграла по переменной z. Прикрепила фото с моими данными что есть, а так все беру со статьи http://cyberleninka.ru/аrticle... ogo-potoka

@bobah16 · 28.03.2016, 14:15

Сообщение от Валентина17

guad для нахождения интеграла по переменной z

еще раз, внимательнее, может quad?

@Валентина17 · 28.03.2016, 14:19 **[ТС]**

Ой, да. Точно!!! Спасибо!!

@Centurio · 28.03.2016, 14:30

Сообщение от Валентина17

я применяла diff для того чтобы вычислить производную от скорости v

Функция diff для массива вычисляет не производную, а разность между соседними элементами. Чтобы найти производную, функцию следует применять к символьному выражению.

@bobah16 · 28.03.2016, 14:40

Валентина17, да, Centurio верно заметил. Нужно еще на шаг поделить, чтобы производную получить.

@Валентина17 · 28.03.2016, 20:04 **[ТС]**

у меня все равно выдает ошибку. не можете исправленный код показать?

@bobah16 · 28.03.2016, 21:35

Сообщение от Валентина17

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

Всю ту же? Если да, то вам уже сказали в чем проблема.

@Валентина17 · 01.04.2016, 18:44 **[ТС]**

здравствуйте. ошибку я исправила,но не могу посчитать интегрaл от S,т.к. получается массив. Так как возможно посчитать тогда интеграл от S? У меня и v1,j,f зависят от z

@Валентина17 · 09.05.2016, 14:30 **[ТС]**

Здравствуйте. Необходимо рассчитать интеграл от функции S(z,t) dz. Пытаюсь рассчитать,но он не выдает мне значения интеграла. и Построить график зависимости интеграла от времени t, I(t). Изначально брала время постоянной. помогите решить эту задачу. [удалено]
Matlab MВыделить код

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
%s=0.07;
v0=3;
x=0.1;
vz=(v0*0.265)/x;
r=1;
c0=0.082;
w=r/c0;
A2=2;
A3=1;
f0=40000;
c=331;
v=0;
v1=0;
f=0;
S=0;
for z=1:0.1:11
    syms z
    v=v+vz*exp(-0.5*(w*z).^2);
    y=v;
    v1=v1+diff(y)/0.1;
    j=A3*v1;
    f=(v/c)*2*f0;
t=1;% дожен меняться от 0 до 11
    S=S.'+A2.*v1.'*exp(-j)*cos(6.28*(f0+f.')*t); 
    I=@(z)(quadl('S+A2.*v1*exp(-j)*cos(6.28*(f0+f)*t)',1,5));
end
 
%plot(I,t);

@Валентина17 · 09.05.2016, 16:04 **[ТС]**

нашла функция, но не могу реализовать нахождение интеграла от этой функции по z. посоветуйте как можно сделать?
- (39750000*z*exp((19875000*z)/(1681*exp((125000*z^2)/1681))))/(1681*exp((125000*z^2)/1681)) - (119250000*z*exp((59625000*z)/(1681*exp((125000*z^2)/1681))))/(1681*exp((125000*z^2)/1681)) - (39750000*z*exp((19875000*z)/(1681*exp((125000*z^2)/1681)))*cos(3994080/(331*exp((125000*z^2)/1681)) + 251200))/(1681*exp((125000*z^2)/1681)) - (119250000*z*exp((59625000*z)/(1681*exp((125000*z^2)/1681)))*cos(7988160/(331*exp((125000*z^2)/1681)) + 251200))/(1681*exp((125000*z^2)/1681))

Новые блоги и статьи Все статьи Все блоги /
Символьное дифференцирование 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, то после закрытия окошка. . .

@SSC $Эксперт по математике/физике$ 3390 / 1913 / 571 Регистрация: 09.04.2015 Сообщений: 5,365
	28.03.2016, 07:46
	Функция diff (стр.12) возвращает массив на единицу меньше чем исходный для диффиренцирования 0

@Валентина17 0 / 0 / 0 Регистрация: 02.06.2015 Сообщений: 11
	28.03.2016, 12:44 [ТС]
	То есть можно увеличить z на одну единицу и тогда будет работать? 0

@Валентина17 0 / 0 / 0 Регистрация: 02.06.2015 Сообщений: 11
	28.03.2016, 14:09 [ТС]
	Есть два графика и исходные данные необходимо построить S, я применяла diff для того чтобы вычислить производную от скорости v, и guad для нахождения интеграла по переменной z. Прикрепила фото с моими данными что есть, а так все беру со статьи http://cyberleninka.ru/аrticle... ogo-potoka 0

@Валентина17 0 / 0 / 0 Регистрация: 02.06.2015 Сообщений: 11
	28.03.2016, 14:19 [ТС]
	Ой, да. Точно!!! Спасибо!! 0

@bobah16 373 / 343 / 42 Регистрация: 14.07.2015 Сообщений: 2,890
	28.03.2016, 14:40
	Валентина17, да, Centurio верно заметил. Нужно еще на шаг поделить, чтобы производную получить. 0

@Валентина17 0 / 0 / 0 Регистрация: 02.06.2015 Сообщений: 11
	28.03.2016, 20:04 [ТС]
	у меня все равно выдает ошибку. не можете исправленный код показать? 0

@Валентина17 0 / 0 / 0 Регистрация: 02.06.2015 Сообщений: 11
	01.04.2016, 18:44 [ТС]
	здравствуйте. ошибку я исправила,но не могу посчитать интегрaл от S,т.к. получается массив. Так как возможно посчитать тогда интеграл от S? У меня и v1,j,f зависят от z 0

@Валентина17 0 / 0 / 0 Регистрация: 02.06.2015 Сообщений: 11
	09.05.2016, 16:04 [ТС]
	нашла функция, но не могу реализовать нахождение интеграла от этой функции по z. посоветуйте как можно сделать? - (39750000zexp((19875000z)/(1681exp((125000z^2)/1681))))/(1681exp((125000z^2)/1681)) - (119250000zexp((59625000z)/(1681exp((125000z^2)/1681))))/(1681exp((125000z^2)/1681)) - (39750000zexp((19875000z)/(1681exp((125000z^2)/1681)))cos(3994080/(331exp((125000z^2)/1681)) + 251200))/(1681exp((125000z^2)/1681)) - (119250000zexp((59625000z)/(1681exp((125000z^2)/1681)))cos(7988160/(331exp((125000z^2)/1681)) + 251200))/(1681exp((125000z^2)/1681)) 0