Вычисление переменной с относительной точностью

@mat_edu · Регистрация: 14.10.2014

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

Я новичок в работе с системой Matlab и столкнулся с одной проблемой. Кто-нибудь мог бы помочь?

Задача такова:
Для квазиодномерного изоэнтропического течения идеального газа в канале переменного сечения есть соотношение между числом Маха М и площадью Потока А:
e=1./M.*(2/(g+1)*(1+(g-1)/2.*M.^2)).^(g+1)/(2*(g-1))
нужно вычислить М (с относительной точностью 0.001) для дозвукового и сверхзвукового течения при e=10 и g=1.4
Пробовал fzero, но функция не имеет корней.

@tvoretsmira · 15.10.2014, 09:30

Matlab M
1
2
3
4
5
6
7
8
9
10
e=10
g=1.4
M=-1:0.01:6;
zer=@(M)-e+1./M.*(2/(g+1)*(1+(g-1)/2.*M.^2)).^(g+1)/(2*(g-1));
plot(M,zer(M))
hold on
grid on
M1=fzero(zer,0.14)
M2=fzero(zer,4.5)
plot([M1,M2],[0,0],'ro')

Вычисление переменной с относительной точностью

@tvoretsmira · 15.10.2014, 09:48

По умолчанию fzero вычисляет результат с точностью eps, разделим эту точность на полученные корни, и вычислим относительную точность:

Matlab M
1
RelAcc=[eps/M1,eps/M2]

@mat_edu · 31.10.2014, 05:54 **[ТС]**

Большое спасибо за помощь, tvoretsmira! К сожалению, я допустил ошибку при написании уравнения, т.к. степень вычисляется в первую очередь, а выражение (g+1)/(2*(g-1)) у меня является степенью, то я должен был занести его в скобки e=1./M.*(2/(g+1)*(1+(g-1)/2.*M.^2)).^((g+1)/(2*(g-1))). Но это не важно, потому что ваше решение в любом случае верно. Еще раз спасибо!

Новые блоги и статьи Все статьи Все блоги /
Очистка реквизитов документа при копировании Maks 09.04.2026 Алгоритм из решения ниже применим как для типовых, так и для нетиповых документов на самых различных конфигурациях. Задача: при копировании документа очищать определенные реквизиты и табличную. . .	модель ЗдравоСохранения 8. Подготовка к разному выполнению заданий anaschu 08.04.2026 https:/ / github. com/ shumilovas/ med2. git main ветка * содержимое блока дэлэй из старой модели теперь внутри зайца новой модели 8ATzM_2aurI	Блокировка документа от изменений, если он открыт у другого пользователя Maks 08.04.2026 Алгоритм из решения ниже реализован на примере нетипового документа, разработанного в конфигурации КА2. Задача: запретить редактирование документа, если он открыт у другого пользователя. / / . . .	Система безопасности+живучести для сервера-слоя интернета (сети). Двойная привязка. Hrethgir 08.04.2026 Далее были размышления о системе безопасности. Сообщения с наклонным текстом - мои. А как нам будет можно проверить, что ссылка наша, а не подделана хулиганами, которая выбросит на другую ветку и. . .
Модель ЗдрввоСохранения 7: больше работников, больше ресурсов. anaschu 08.04.2026 работников и заданий может быть сколько угодно, но настроено всё так, что используется пока что только 20% kYBz3eJf3jQ	Дальние перспективы сервера - слоя сети с космологическим дизайном интефейса карты и логики. Hrethgir 07.04.2026 Дальнейшее ближайшее планирование вывело к размышлениям над дальними перспективами. И вот тут может быть даже будут нужны оценки специалистов, так как в дальних перспективах всё может очень сильно. . .	Горе от ума kumehtar 07.04.2026 Эта мне ментальная установка, что вот прямо сейчас, мол, мне для полного счастья не хватает (нужное вписать), и когда я этого достигну - тогда и полный кайф. Одна из самых сильных ловушек на пути. . . .	Использование значений реквизитов справочника в документе, с определенными условиями и правами Maks 07.04.2026 1. Контроль срока действия договора Алгоритм из решения ниже реализован на примере нетипового документа "ЗаявкаНаРаботу", разработанного в конфигурации КА2. Задача: уведомлять пользователя, если. . .

@mat_edu 0 / 0 / 0 Регистрация: 14.10.2014 Сообщений: 6

	Вычисление переменной с относительной точностью 15.10.2014, 00:12. Показов 1150. Ответов 3 Метки нет (Все метки) Я новичок в работе с системой Matlab и столкнулся с одной проблемой. Кто-нибудь мог бы помочь? Задача такова: Для квазиодномерного изоэнтропического течения идеального газа в канале переменного сечения есть соотношение между числом Маха М и площадью Потока А: e=1./M.(2/(g+1)(1+(g-1)/2.M.^2)).^(g+1)/(2(g-1)) нужно вычислить М (с относительной точностью 0.001) для дозвукового и сверхзвукового течения при e=10 и g=1.4 Пробовал fzero, но функция не имеет корней. 0

@mat_edu 0 / 0 / 0 Регистрация: 14.10.2014 Сообщений: 6
	31.10.2014, 05:54 [ТС]
	Большое спасибо за помощь, tvoretsmira! К сожалению, я допустил ошибку при написании уравнения, т.к. степень вычисляется в первую очередь, а выражение (g+1)/(2(g-1)) у меня является степенью, то я должен был занести его в скобки e=1./M.(2/(g+1)(1+(g-1)/2.M.^2)).^((g+1)/(2*(g-1))). Но это не важно, потому что ваше решение в любом случае верно. Еще раз спасибо! 0