Поиск с возвращением
Курсовой проект - Физика
Другие курсовые по предмету Физика
лгоритм ручного расчета законами Кирхгофа. Данные, полученные в схемном эмуляторе ненамного отличаются от результатов работы программы, т.к. в эмуляторе используются элементы, близкие к реальным, а в программном растете элементы идеальные.
Тем не менее, программная реализация расчета цепи значительно ускорит процесс расчета, повысит эффективность работы и упростит анализ результатов.
Список используемой литературы
1. Теоретические основы электротехники: Учебник для вузов. Том 1. - 4-е изд. /К. С. Демирчан, Л. Р. Нейман. СПб.:Питер,2004. - 463с.:ил.
. Джон Г. Мэтьюз, Куртис Д. Финк. Численные методы. Использование Matlab. - М.: Издательский дом Вильямс, 2001. - 720 с.
Приложение
Модуль данных - файл data.m:
%модуль данных
%данные цепи
%время=0.25^10-3
%источник тока=1
%резисторы=200=200=200=200=200=200
%циклическая частота=2*3.14*10^6
%емкость конденсатора=5*10^-6=5*10^-6=5*10^-6=5*10^-6
%индуктивность катушки=2*10^-7=2*10^-7
%напряжение=100=150=200
Модуль подготовки данных для расчета - файл rasch.m:
n=11=zeros(n,n);q=1:nw=1:n(q,w)=0.000000000001;=[1 -1 0.000000001 0.000000001 0.000000001 0.000000001 0.000000001 0.000000001 0.000000001 0.000000001 -1;%1
.0000000001 1 -1 -1 0.000000001 0.000000001 0.000000001 0.000000001 0.000000001 0.000000001 0.000000001;%2
.000000001 0.000000001 1 0.000000001 0.000000001 0.000000001 0.000000001 0.000000001 0.00001 1 1%7
.000000001 0.000000001 0.000000001 1 -1 -1 0.000000001 0.000000001 0.000000001 0.000000001 0.000000001;%3
.000000001 0.000000001 0.000000001 0.000000001 1 0.000000001 0.000000001 1 1 -1 0.000000001;%6
.000000001 0.000000001 0.000000001 0.000000001 0.000000001 1 1 0.000000001 -1 0.000000001 0.000000001;%4
0.000000001 0.000000001 0.000000001 0.000000001 0.000000001 1 1 0.000000001 0.000000001 0.000000001;%5
];k=1:nj=1:7(j,k)=AK1(j,k);qw=1:n(qw,1)=0.000000001;(4,1)=J/(2^(1/2));(7,1)=J/(2^(1/2));(9,5)=-(i/(am*C5)+R5);(9,6)=R6;(9,9)=am*L9*i;(10,4)=R4;(10,5)=i/(am*C5)+R5;(10,10)=R10;(10,1)=-E3/(2^(1/2));(8,7)=i/(am*C7);(8,9)=am*L9*i;(8,8)=-R8;(11,8)=R8;(11,10)=R10;(11,11)=-(i/(am*C11)+R11);(11,1)=E1/(2^(1/2)) ;
Модуль расчета детерминанта компактной схемой Гаусса - determinant.m:
function D=determinant(A,n)
%определение определителя
=zeros(n,n);Al=zeros(n,n);
(1,:)=A(1,:);k=2:n,
Al(k,1)=A(k,1)/Y(1,1);l=2:n,=0;j=l:n,=0;k=1:l-1,=sum+(Al(l,k)*Y(k,j));(l,j)=A(l,j)-sum;j=l+1:n=0;k=1:l-1=sum+(Al(j,k)*Y(k,l));(j,l)=1/Y(l,l)*(A(j,l)-sum);l=1:nj=l:n(l,j)=Y(l,j);l=2:nj=l+1:n(j,l)=Al(j,l);=AA(1,1);l=2:n,=D*AA(l,l);;
Главный модуль - файл main.m
data
%Проверим невырожденность системы(A);
%По правилу Крамера=zeros(1,n);=determinant(A,n);i=1:n= A;(:,i)=B;=determinant(A1,n);=D1/D;(1,i)=sqrt(imag(x1)^2+real(x1)^2);=atan(imag(x1)/real(x1));(i)=x(1,i)*sqrt(2)*sin(am*t+fi);