Рекомендации к выполнению модернизированной лабораторной работы №1
(для развития навыков составления сценариев в среде MatLab
c целью оформления результатов в документе редактора Word)
Это файл, который берется в обработку с целью построения характеристик потокосцепления в зависимости от тока в обмотке и рабочего зазора электромагнита, а затем и получения характеристик его тяговых сил.
Данные получены в результате расчета энергии магнитного поля в конечно-разностной форме средствами Fortran(a) и последующего пересчета к потокосцеплению (это выполняется долго, трудно, поэтому и точек мало ).
6 % Количество точек по току
6 % Количество точек по зазору
3.000000000000000 % Зазор, мм
5.000000000000000
7.000000000000000
9.000000000000002
11.000000000000000
13.000000000000000
0.000000000000000E+000 % Ток, А
2.000000000000000E-001
4.000000000000000E-001
6.000000000000000E-001
8.000000000000000E-001
1.000000000000000
0.000000000000000E+000 % Потокосцепление, Вб (заданы по мере нарастания тока,
1.625374659634170E-001 % при переходе от минимального зазора к максимальному)
1.908930753149628E-001
1.916477140434351E-001
1.929942836417339E-001
1.948794733596508E-001
0.000000000000000E+000
1.360330138427829E-001
1.778294946315807E-001
1.924191242633246E-001
1.933254460428406E-001
1.948180015537299E-001
0.000000000000000E+000
1.110193061643157E-001
1.622445101564819E-001
1.833204160114363E-001
1.840943203141589E-001
1.854084720405954E-001
0.000000000000000E+000
9.429297291277157E-002
1.466409128641101E-001
1.761291910832546E-001
1.873534947397786E-001
1.881066357463121E-001
0.000000000000000E+000
8.257826389968255E-002
1.312141505896361E-001
1.645784204287423E-001
1.788499754697633E-001
1.795624482418855E-001
0.000000000000000E+000
7.360346534203498E-002
1.170809577019471E-001
1.527694072684852E-001
1.699037922511276E-001
1.705798567020393E-001
Далее – по сценарию на рабочем столе среды MatLab.
Формируем вектор-столбец значений тока.
Данные переносятся из приведенного выше файла копированием в буфер обмена.
i=[0.000000000000000E+000
2.000000000000000E-001
4.000000000000000E-001
6.000000000000000E-001
8.000000000000000E-001
1.000000000000000];
Захотели убедиться, пожалуйста, смотрите на рабочий стол среды MatLab. Данные в рабочем пространстве.
» i
i =
0
0.2000
0.4000
0.6000
0.8000
1.0000
»
Это самый короткий формат представления данных в среде MatLab (см. preferences в выпадающем меню).
Составляя протокол параллельно нужно даные перетащить и в таблицу помня о том, что в инженерной практике при оформлении ставят 3 знака после запятой.
Формируем вектор-столбец значений зазора и не забываем о таблице в протоколе.
z=[3.000000000000000
5.000000000000000
7.000000000000000
9.000000000000002
11.000000000000000
13.000000000000000];
Если желаете, то убедитесь, что вектор-столбец для зазора тоже существует в рабочем пространстве.
Теперь можно взять вектор значений потокосцепления
p=[0.000000000000000E+000
1.625374659634170E-001
1.908930753149628E-001
1.916477140434351E-001
1.929942836417339E-001
1.948794733596508E-001
0.000000000000000E+000
1.360330138427829E-001
1.778294946315807E-001
1.924191242633246E-001
1.933254460428406E-001
1.948180015537299E-001
0.000000000000000E+000
1.110193061643157E-001
1.622445101564819E-001
1.833204160114363E-001
1.840943203141589E-001
1.854084720405954E-001
0.000000000000000E+000
9.429297291277157E-002
1.466409128641101E-001
1.761291910832546E-001
1.873534947397786E-001
1.881066357463121E-001
0.000000000000000E+000
8.257826389968255E-002
1.312141505896361E-001
1.645784204287423E-001
1.788499754697633E-001
1.795624482418855E-001
0.000000000000000E+000
7.360346534203498E-002
1.170809577019471E-001
1.527694072684852E-001
1.699037922511276E-001
1.705798567020393E-001
];
Превратим столбец р в матрицу рр, указав предварительно число строк (по току) и число столбцов (по зазору).
ni=6; nz=6;
pp=reshape(p,ni,nz); % Это выполнено с применением функции MatLab. Подробности можете
% узнать, набрав help reshape. Посмотрите и в документ перетащите!
Размножили узлы для интерполяции (любопытство удовлетворите через help linspace)
» ii=linspace(i(1),i(end)); % по току 100 точек
» zz=linspace(z(1),z(end)); % по зазору 100 точек
Любопытство здесь и далее удовлетворяйте через help CSAPS.
» P=[1 1]; % исключили сглаживание м.н.к
» VALUES = CSAPS({i,z},pp,P,{ii,zz}); % выполнили двумерную интерполяцию
% кубическим сплайном, результат в VALUES
» figure(1); plot(ii,VALUES(:,1:10:end)); % выдали график по току (10 линий)
» figure(2); plot(zz,VALUES(1:10:end,:)); % выдали график по зазору (10 линий)
Оформить в документе согласно ГОСТ 7.32-81 ( Отчет о научно-исследовательской работе. Общие требования к правилам оформления).
Оформление графика выполнять по возможности полно находясь в MatLab(е) через меню фигуры.
Для переноса в протокол фигуру необходимо скопировать, а затем уж вставить, прибегнув к процедуре специальной вставки и работы с командами формата рисунка. В протоколе сделайте все правильно под ГОСТ.
Найдите различия в оформлении графиков а и б. Примите решение как будете оформлять (с легендой либо выносками).
а |
б |
Здесь разместить подрисуночную надпись.Это один рисунок, но в нем дла вида! |
Чтобы прочувствовать одну из проблем автоматизированной обработки результатов далее получите и оформите графики, соответствующие выражению
, где ;.
В приложение к отчету представить олаженную программу. Иметь на дискете.
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.