|
||||
Меню:
Главная
Форум
Литература: Программирование и ремонт Импульсные блоки питания Неисправности и замена Радиоэлектронная аппаратура Микросхема в ТА Рубрикатор ТА Кабельные линии Обмотки и изоляция Радиоаппаратура Гибкие диски часть 2 часть 3 часть 4 часть 5 Ремонт компьютера часть 2 Аналитика: Монтаж Справочник Электроника Мощные высокочастотные транзисторы 200 микросхем Полупроводники ч.1 Часть 2 Алгоритмические проблемы 500 микросхем 500 микросхем Сортировка и поиск Монады Передача сигнала Электроника Прием сигнала Телевидиние Проектирование Эвм Оптимизация Автомобильная электроника Поляковтрансиверы Форт Тензодатчик Силовые полевые транзисторы Распределение частот Резисторные и термопарные Оберон Открытые системы шифрования Удк |
[104] (define dy (stream-map f y)) y) Эта процедура не работает, потому что вызов integral в первой строке solve требует, чтобы был определен входной поток dy, а это происходит только во второй строке процедуры solve. С другой стороны, замысл, заключенный в этом определении, вполне здрав, поскольку мы можем, в принципе, начать порождать поток y и не зная dy. Действительно, integral и многие другие операции над потоками обладают свойствами, подобными cons-stream, а именно, мы можем породить часть ответа, даже если нам дана только частичная информация об аргументах. В случае integral, первый элемент выходного потока есть указанное начальное значение initial-value. Таким образом, можно породить первый элемент выходного потока и не вычисляя интегрируемую величину dy. А раз мы знаем первый элемент y, то stream-map во второй строке solve может начать работать и породить первый элемент dy, а с его помощью мы получим второй элемент y, и так далее. Чтобы воспользоваться этой идеей, переопределим integral так, чтобы он ожидал интегрируемый поток в виде задержанного аргумента (delayed argument). Integral будет размораживать вычисление входного потока через force только тогда, когда ему нужно породить элементы входного потока помимо первого: (define (integral delayed-integrand initial-value dt) (define int (cons-stream initial-value (let ((integrand (force delayed-integrand))) (add-streams (scale-stream integrand dt) int)))) int) Теперь можно реализовать процедуру solve, задержав вычисление dy внутри определения y:71 (define (solve f y0 dt) (define y (integral (delay dy) y0 dt)) (define dy (stream-map f y)) y) Теперь при любом вызове integral необходимо задерживать интегрируемый аргумент. Можно показать, что процедура solve работает, аппроксимируя e « 2.718 вычислением в точке y = 1 решения дифференциального уравнения dy/dt = y с начальным условием y(0) = 1: (stream-ref (solve (lambda (y) y) 1 0.001) 1000) 2.716924 71 Не гарантируется, что эта процедура будет работать во всех реализациях Scheme, но для любой реализации должен найтись простой способ заставить подобную процедуру работать. Проблемы связаны с тонкими различиями в том, как реализации Scheme обрабатывают внутренние определения. (См. раздел 4.1.6.) dyo yo Рис. 3.35: Диаграмма потока сигналов для решения линейного дифференциального уравнения второго порядка. Упражнение 3.77. Вышеприведенная процедура integral была аналогична «непрямому» определению бесконечного потока натуральных чисел из раздела 3.5.2. В виде альтернативы можно дать определение integral, более похожее на integers-starting-from (также в разделе 3.5.2): (define (integral integrand initial-value dt) (cons-stream initial-value (if (stream-null? integrand) the-empty-stream (integral (stream-cdr integrand) (+ (* dt (stream-car integrand)) initial-value) dt)))) В системах с циклами эта реализациея порождает такие же проблемы, как и наша исходная версия integral. Модифицируйте процедуру так, чтобы она ожидала integrand как задержанный аргумент, а следовательно, могла быть использована в процедуре solve. Упражнение 3.78. Рассмотрим задачу проектирования системы обработки сигналов для решения гомогенных линейных дифференциальных уравнений второго порядка d2y dy Выходной поток, моделирующий у, порождается сетью, содержащей цикл. Этот цикл возникает потому, что значение d2y/dt2 зависит от значений у и dy/dt, а они оба получаются интегрированием d2y/dt2. Диаграмма, которую нам хотелось бы закодировать, показана на рис. 3.35. Напишите процедуру solve-2nd, которая в качестве аргументов берет константы a, b и dt и начальные значения y0 и dy0 для y и dy, и порождает поток последовательных значений y. iR + VR + Vc -aav R C iL VL Рис. 3.36: Последовательная RLC-цепь Упражнение 3.79. Обобщите процедуру solve-2nd из упражнения 3.78 так, чтобы с ее помощью можно было решать дифференциальные уравнения второго порядка общего вида d2y/dy2 = f (dy/dt,y). Упражнение 3.80. Последовательная RLC-цепь (series RLC circuit) состоит из резистора, конденсатора и катушки индуктивности, соединенных последовательно, как показано на рис. 3.36. Если сопротивление, индуктивность и емкость равны, соответственно, R, L и C, то отношения между напряжением v и током i на трех элементах описываются уравнениями vr = vl а цепь диктует соотношения ic iR = vc =L =C iRR diL ~dt dvc iL = -ic = vl + VR Сочетание этих условий показывает, что состояние цепи (характеризуемое через vc, напряжение на конденсаторе, и iL, ток через катушку) описывается парой дифференциальных уравнений dvc diL 1 dt L vc С R L iL Диаграмма потока сигналов, представляющая эту систему дифференциальных уравнений, показана на рисунке 3.37. Напишите процедуру RLC, которая в качестве аргументов берет параметры цепи R, L и C и точность по времени dt. Подобно процедуре RC из упражнения 3.73, RLC должна порождать процедуру, которая берет начальные значения переменных состояния vCo и iLo и порождает (через cons) пару потоков состояния vc и iL. С помощью RLC породите пару потоков, которая моделирует поведение RLC-цепи c K = 1 ом, C = 0.2 фарад, L =1 генри, dt = 0.1 секунды, и начальными значениями iLo =0 ампер и vco 10 вольт. i C + |
Среды: Smalltalk80 MicroCap Local bus Bios Pci 12С ML Микроконтроллеры: Atmel Intel Holtek AVR MSP430 Microchip Книги: Емкостный датчик 500 схем для радиолюбителей часть 2 (4) Структура компьютерных программ Автоматическая коммутация Кондиционирование и вентиляция Ошибки при монтаже Схемы звуковоспроизведения Дроссели для питания Блоки питания Детекторы перемещения Теория электропривода Адаптивное управление Измерение параметров Печатная плата pcad pcb Физика цвета Управлении софтверными проектами Математический аппарат Битовые строки Микроконтроллер nios Команды управления выполнением программы Перехода от ahdl к vhdl Холодный спай Усилители hi-fi Электронные часы Сердечники из распылённого железа Анализ алгоритмов 8-разрядные КМОП Классификация МПК История Устройства автоматики Системы и сети Частотность Справочник микросхем Вторичного электропитания Типы видеомониторов Радиобиблиотека Электронные системы Бесконтекстный язык Управление техническими системами Монтаж печатных плат Работа с коммуникациями Создание библиотечного компонента Нейрокомпьютерная техника Parser Пи-регулятор ч.1 ПИ-регулятор ч.2 Обработка списков Интегральные схемы Шина ISAВ Шина PCI Прикладная криптография Нетематическое: Взрывной автогидролиз Нечеткая логика Бытовые установки (укр) Автоматизация проектирования Сбор и защита Дискретная математика Kb радиостанция Энергетика Ретро: Прием в автомобиле Управление шаговым двигателем Магнитная запись Ремонт микроволновки Дискретные системы часть 2 | ||