Занятия по MATLAB в компьютерном классе
Занятия по MATLAB в компьютерном классе (6 занятий по 90 мин.) Темы занятий 1. Общие сведения о MATLAB'е. 2. Переменные. 3. Элементы xy-графики. 4. Простые примеры, иллюстрирующие эффективность MATLAB'а. 5. Графический способ решения уравнений. 5. Полиномы. 7. Итерации. 8. Системы линейных алгебраических уравнений. Литература. Тема 1. Общие сведения о MATLAB'е MATLAB - матричная лаборатория - наиболее развитая система программирования для научно-технических расчетов, дополненная к настоящему времени несколькими десятками более частных приложений, относящихся к вычислительной математике, обработке информации, конструированию электронных приборов, экономике и ряду других разделов прикладной науки. Изучение MATLAB'а по фирменной документации, которая теперь прилагается на инсталляционном компакт-диске, занимает у начинающих пользователей слишком много времени не только из-за необходимости читать ее на английском языке со специфическим сленгом, но, главным образом, ввиду неизбежного для таких руководств последовательного и достаточно формального изложения большого объема информации, а имеющиеся на русском языке пособия в основном следуют этому стереотипу. Даже для опытного специалиста по расчетам на компьютере такое изучение сопряжено с неоправданно большими затратами труда. MATLAB предназначен прежде всего для программирования численных алгоритмов. Он разрабатывается уже более 15 лет и возник на основе более ранних прикладных пакетов LINPACK и EIGPACK, созданных в 1970-е гг. в США, и в свою очередь повлиял на появление таких систем, как MathCad, MAPLE и Mathematica. Совершенствование системы MATLAB происходило как в связи с достижениями в вычислительной математике, так и в связи с изменениями в архитектуре персональных компьютеров и развитием общесистемных средств. Со временем MATLAB был дополнен целым рядом уже упоминавшихся приложений (toolboxes), далеко раздвинувших границы его применимости. Далее речь пойдет лишь о ядре MATLAB'а, которое мы будем называть системой, и конкретно о ее версии 5.2, выпущенной фирмой MathWorks в январе 1998 г. MATLAB - система программирования высокого уровня, работающая как интерпретатор и включающая большой набор инструкций (команд) для выполнения самых разнообразных вычислений, задания структур данных и графического представления информации. Команды эти разбиты на тематические группы, расположенные в различных директориях системы. Теперь в системе около 800 команд, и примерно половина из них вполне доступна начинающему пользователю. Команды с большим возможным объемом вычислений написаны на С, но много и таких команд, которые представлены в терминах этих первых. Поэтому система оказывается почти открытой для пользователя. Имеются большие возможности для вывода двумерной и трехмерной графики и средства управления ею. Пользователь может без особых затруднений добавлять свои команды и писать программы в терминах уже существующих команд; несколько сложнее делать это в рамках Фортрана и С. Можно обмениваться данными с программами на этих языках, а из них обращаться к системе. Краткость и наглядность программирования и исключительные возможности визуализации результатов делают систему очень эффективной при поисках и апробации новых алгоритмов, при проведении разовых расчетов и в учебном процессе, поскольку ее можно осваивать без предварительного знакомства с основами программирования и выполнять такие сложные примеры, которые невозможно делать с использованием других систем. Документация по системе и ее приложениям содержит много тысяч страниц, и поэтому естественно встает вопрос о том, как ее осваивать. Работа с системой требует определенной математической подготовки, так что обучение можно начинать на втором курсе вуза. Основные сведения о системе изложены в руководствах /1/ - /2/: /1/ - это учебник с описанием вычислительных возможностей и архитектуры системы, /2/ - описание ее графических возможностей. Конечно, можно читать подряд /1/, /2/ и при необходимости обращаться за уточнениями к команде help или справочнику /3/, в котором описаны почти все команды. Но гораздо более эффективным, на наш взгляд, является изложение основных вычислительных процедур с помощью наиболее употребительных команд системы. Именно так мы и познакомимся с MATLAB'ом, а точнее примерно с 30-40 его командами. После этих занятий пользоваться документацией /1/ и /2/ будет гораздо легче. Приложений к системе мы здесь не касаемся, а изучать их можно только после предварительного ознакомления с ней, а также с тем разделом науки, которому посвящено конкретное приложение. Отметим только, что большинство приложений означают для пользователя просто расширение списка доступных ему команд. Очень удобно то, что вся документация по системе и приложениям находится на компакт-диске, с которого происходит их установка, и при желании она может быть размещена также и на винчестере. Для работы с системой достаточно иметь компьютер PC 486 с оперативной памятью хотя бы 16 Mb и с установленными на нем системами Windows 95 и MATLAB 5.2. В действительности MATLAB может работать и с другими операционными системами, такими, например, как Macintosh, Unix и OS/2. За рубежом вышло уже достаточно много учебных пособий по системе, но на русский язык ни одно из них пока не переводилось и даже в центральных библиотеках их теперь нет из-за сокращения финансирования. Изданные у нас пособия (например, /4/ - /12/) в основном следуют руководствам /1/ - /3/, тогда как нам представляется полезным дать менее формальное введение в предмет, опираясь прежде всего на интуицию слушателя. MATLAB - высокоуровневая система программирования, позволяющая резко сократить затраты труда при проверке алгоритмов и проведении прикидочных расчетов. Возможность проведения больших расчетов на MATLAB'е определяется в основном теми затратами времени, на которые может пойти пользователь: здесь приходится выбирать между легкостью и наглядностью программирования и представления результатов, с одной стороны, и затратами времени на счет - с другой. Система очень удобна для освоения и апробации численных методов, что мы и хотим показать здесь прежде всего. Именно поэтому она рекомендуется как одна из основных для физиков и многих других естественно-научных специальностей в ведущих американских университетах. Детальное освоение любой большой программной системы - это достаточно длительный процесс, основу которого составляют индивидуальная работа, и наши занятия призваны дать лишь первоначальный импульс этому процессу в отношении MATLAB'а. Темы 2 - 4 представляют сравнительно элементарное введение, а в остальных рассматриваются более сложные примеры, показывающие, как можно использовать программные и графические возможности системы для исследования численных алгоритмов. В некоторых местах текста вставлены упражнения, пронумерованные сквозным образом. Тема 2. ПеременныеПеременные могут быть числовыми, текстовыми и других типов. У нас будут только числовые (это во всех деталях) и текстовые (совсем немного). Название переменной начинается с латинской буквы, далее могут быть буквы и цифры (не более 31 символа). Строчные и прописные буквы здесь различаются. 1. Числовые переменные. Это числа, векторы, матрицы и многомерные массивы. В компьютере все числа представлены примерно с 16 десятичными знаками, под каждое вещественное число отводится 8 байтов, под комплексное - 16. 1.1. Ввод чисел. Целые числа. В системе они не выделяются явно. Наберем и выполним отдельно каждую команду:a=2 a=2.0 a=2; a=1:6 b=1:20 c=10:-2:5 Командное окно. Командная строка. Редактирование командной строки. Буфер исполненных команд. Как выбирать информацию из командного окна и из буфера исполненных командных строк. Нельзя допускать совпадения имени переменной с именем какой-либо команды. Вещественные числа. Выполним по отдельности следующие команды:d=0.5:0.3:2.5 d=.5:.3:2.5 d=.5+1:.3-.1:2.5*2 length(d) d(end) d(end-2) d(1) d(0) d(2:7) d(7:-1:2) d(150) f=linspace(1.5,30,143); length(f)Индексы всегда начинаются со значения 1. Команды набираются на малом латинском регистре. Возможна многопараметричность команд. Диапазон вещественных чисел: realmax realminДругие константы MATLAB'а: pi i j eps Их не следует портить. Комплексные числа: q=1+2*i q=1+2i real(q) imag(q) abs(q) conj(q) s=angle(q) (здесь -<s<=).q=1+2*i;r=3; fi=0:.01:pi; z=q+r*exp(i*fi); plot(z) Это верхняя полуокружность. 1.2. Ввод векторовВекторы-строки: a=1:6 linspace(1,6,10)Векторы-столбцы: a=(1:6)' linspace(1,6,10)' Операторы .' и ' :y1=linspace(1,6,4)'; y2=y1; y=y1+i*y2; y.' y' Команды linspace и : применимы для задания только вещественных векторов. 1.3. Ввод матриц. A(i,j) - элемент из i-й строки и j-го столбца. A(k) - k-й элемент таблицы, вытянутой в столбец. A=[1,2;3,4] A=[1;2,3;4] A(2,2) A(3) A(5) size(A) A(3,4)=10 size(A) A(5)=6 size(A) A(22)=3 A=A(:) A(22)=3 size(A) [m,n]=size(A)A=reshape(1:24,4,6) size(A) A([1,end],:)=[] A(:,[1,end])=[] size(A) 1.4. Некоторые специальные матрицы: m=3;n=4; eye(m,n) eye(m) eye(n) ones(m,n) ones(m) ones(n) zeros(m,n) rand(m,n) rand(m,n) rand('state',0) rand(m,n) rand(m) Это равномерное распределение на интервале (0, 1).randn(m,n) randn('state',100) Это нормальное распределение, у него мат.ожидание=0, дисперсия=1v1=1:4 v2=7:12 toeplitz(v1,v2) toeplitz(v1) 1.5.Некоторые простые команды:A=reshape(1:24,4,6) triu(A) triu(A,0) triu(A,2) triu(A,-1) tril(A) v=1:5 diag(v) diag(v,2) diag(v,-1) diag(A) diag(A,2) diag(A,-1) A=reshape(1:24,4,6) rot90(A) rot90(A,2) 1.6. Выдачи на экран. Команда format с различными опциями. В обычном формате (forrmat short) выдается 5 знаков, для целых чисел 9 знаков, порядки изменяются от -308 до +308. В полном формате (format long e) 16 знаков.a=2 a=.001 a=1e-3 a=1e-5 a=123456789 a=1234567891 a=1+3*i format long e, 2^.5, format short Опция format short e позволяет получать ровные столбцы. 2. Текстовые переменные. Они берутся в кавычки (на букве э на латинском регистре), символ занимает 2 байта. Используются для задания заголовков в числовых выдачах и на графиках, для задания формул и т.д. Можно переводить текстовые переменные в числовые и наоборот. Выполним в командной строкеt='Moscow - столица России' (после дефиса нужно перейти на русский шрифт и затем не забыть снова вернуться на латинский). 3. Контроль за переменными:who whos 4. Другие типы переменных - ячейки и структуры.5. Система help. help выдает список директорий системы;help <имя директории> выдает список команд директории;help <имя команды> выдает описание команды. type <имя команды> выдает текст команды или программы пользователя, если он составлен в терминах MATLAB'а. Упражнения1. Из заданной матрицы A выбрать вектор, компоненты которого есть все окаймляющие A элементы, взятые в порядке A(1,1), A(2,1), … , A(1,2).2. Из заданной матрицы A размера m*n построить матрицу B с m строками, у которой диагонали с номерами 0, 1,…, n-1 были бы столбцами A с номерами 1:n, а все остальные элементы равнялись бы нулю. 3. Построить последовательность из 1000 целых случайных чисел, которые равновероятно принимали бы значения 11:18, и выдать, сколько из них оказалось равным 11, 12, …, 18. Используйте команду round(k) округления k до ближайшего целого.Тема 3. Элементы xy-графики1.Как открывать графическое окно: figure whitebg zoom onТеперь построим график функципи y=sin(2x), 0<=x<=5, выполнив строкуx=0:1e-3:5; y=sin(2*pi*x); plot(y) plot(x,y) ,grid
Страницы: 1, 2, 3, 4, 5
|
