0 / 0 / 0
Регистрация: 25.11.2013
Сообщений: 7
1

Построение АЧХ и ФЧХ

25.11.2013, 02:55. Показов 7711. Ответов 7
Метки нет (Все метки)

Здраствуйте.. помогите пожалуйста с кодом.. пишет ошибку "Index exceeds matrix dimensions." но я не могу разобраться.. в Matlabе
Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
% Находим передаточные характеристики
maxK = 0;
nn = 1;
a1=0;
a=0;
dt = 10000;
for F = Fr-Fr/2:Fr/dt:Fr+Fr/2;
    w=2*pi*F;
 
a1_n=[Z11+R1+Z12,0,0;
    0,Z22+i*Wr*L1+Z21,-i*Wr*M;
    0,-i*Wr*M,i*Wr*L2+r+(1/i*Wr*C1)+Zn];
 
a_n=[Z11+R1+Z12,0,1;
    0,Z22+i*Wr*L1+Z21,0;
    0,-i*Wr*M,0];
 
K_n(nn)=Zn*det(a_n)/det(a1_n);
K0_n(nn)=abs(Zn*det(a_n)/det(a1_n));
faza_n(nn)=angle(Zn*det(a_n)/det(a1_n))/(pi)*180;
 if maxK<=K0_n(nn)
  maxK=K0_n(nn);
 end
 if maxK<=K0_n(nn)
  maxN=nn;
 end
 FF(nn)=F;
 nn=nn+1;
end
%--------------------------------------------------------------------------
%находим полосу пропускания
for NN = 1:1:maxN;
if K0_n(NN)<=0.707*maxK
f1=NN;
end 
end
NN=0;
for NN = maxN:1:nn-1;
if K0_n(NN)>=0.707*maxK
 f2=NN;
end
end
fn=f1;
fw=f2;
fmin07=FF(f1);
fmax07=FF(f2);
%--------------------------------------------------------------------------
% Строим Ачх и Фчх
figure, plot(FF(fn-50:fw+50),K0_n(fn-50:fw+50));grid;title('АЧХ');
figure, plot(FF(fn-200:fw+200),faza_n(fn-200:fw+200));grid;title('ФЧХ');
__________________
Помощь в написании контрольных, курсовых и дипломных работ здесь
0
Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
25.11.2013, 02:55
Ответы с готовыми решениями:

Построение АЧХ и ФЧХ сигнала
Здравствуйте! У меня есть следующий код: handle = fvtool(h); set(handle,'Analysis','freq')...

АЧХ и ФЧХ сигнала
Привет всем, хотел узнать, правильно ли я ищу АЧХ и ФЧХ заданного сигнала(задан тремя функциями). ...

семейство АЧХ,ФЧХ
Дано несколько передаточных функций h1= 6e06 ------------------------ ,...

Построить годограф, АЧХ и ФЧХ
Пожалуйста помогите построить графики в Matlab, очень надо. Весь интернет перерыла, но никаких...

7
5172 / 3510 / 361
Регистрация: 02.04.2012
Сообщений: 6,422
Записей в блоге: 17
26.11.2013, 13:30 2
А можешь скинуть недостающие данные (Fr, Wr, Rx, Zxx, Cx), чтобы можно было запустить и посмореть, что к чему
0
0 / 0 / 0
Регистрация: 25.11.2013
Сообщений: 7
26.11.2013, 18:13  [ТС] 3
Цитата Сообщение от Зосима Посмотреть сообщение
А можешь скинуть недостающие данные (Fr, Wr, Rx, Zxx, Cx), чтобы можно было запустить и посмореть, что к чему
да, вот
Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
%Наменалы элементов
R1=180*10^3;
L1=150*10^-6;
L2=800*10^-6;
C1=910*10^-12;
C3=51*10^-6;
C4=2.2*10^-12;
Rn=150*10^3;
Cn=30*10^-12;
M=50*10^-6;
Q=100;
Czi=10*10^-12;
Czs=5*10^-12;
Ri=180*10^3;
S=4*10^-3;
ru=200;
%--------------------------------------------------------------------------
%Резонансная частота
C1n=C1+Cn;
Wr=1/(sqrt(L2*C1n));
Fr=Wr/(2*pi)
%Граничная частота
Fs=(1+S*ru)/(2*pi*Czi*ru);
%--------------------------------------------------------------------------
%Расчёт Y параметров
Y11=j*Wr*(Czi+Czs);
Y12=-j*Wr*Czs;
Y21=S-j*Wr*Czs;
Y22=1/Ri+j*Wr*Czs;
Y=[Y11, Y12;
Y21, Y22];
%Расчёт Z параметров
Z11=Y22/(det(Y));
Z12=Y12/(det(Y));
Z21=Y21/(det(Y));
Z22=Y11/(det(Y));
Z=[Z11,Z12;
Z21,Z22];
%--------------------------------------------------------------------------
Xc3=1/j*Wr*C3;
Xc4=1/j*Wr*C4;
%коэффициент включения контура
p=L2/L1;
%определяем сопротивление потерь контура
ro=sqrt(L2/C1)*Q; %-характерестическое сопротивление
r=Wr*L2/Q ;%-сопротивление потерь
0
Эксперт по электронике
938 / 838 / 121
Регистрация: 23.11.2012
Сообщений: 2,489
26.11.2013, 20:33 4
Toshka123, а Вы не хотите привести схему и огласить, что, собственно, нужно сделать?
Глядишь проще получится. Или Вам нужно сделать так_как_сказал_преподаватель (с)?
0
0 / 0 / 0
Регистрация: 25.11.2013
Сообщений: 7
26.11.2013, 23:37  [ТС] 5
Цитата Сообщение от R2D2 Посмотреть сообщение
Toshka123, а Вы не хотите привести схему и огласить, что, собственно, нужно сделать?
Глядишь проще получится. Или Вам нужно сделать так_как_сказал_преподаватель (с)?
там правда вместо с3 с2
Вложения
Тип файла: docx курсач.docx (62.4 Кб, 32 просмотров)
0
5172 / 3510 / 361
Регистрация: 02.04.2012
Сообщений: 6,422
Записей в блоге: 17
27.11.2013, 18:16 6
Вот набросал, но что-то не то... определитель det(a_n) получается равным 0, поэтому так и получается Ищи, смотри.
Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
clear, clc
 
R1=180*10^3;
L1=150*10^-6;
L2=800*10^-6;
C1=910*10^-12;
C3=51*10^-6;
C4=2.2*10^-12;
Rn=150*10^3;
Cn=30*10^-12;
M=50*10^-6;
Q=100;
Czi=10*10^-12;
Czs=5*10^-12;
Ri=180*10^3;
S=4*10^-3;
ru=200;
%--------------------------------------------------------------------------
%Резонансная частота
C1n=C1+Cn;
Wr=1/(sqrt(L2*C1n));
Fr=Wr/(2*pi);
%Граничная частота
Fs=(1+S*ru)/(2*pi*Czi*ru);
j = sqrt(-1);
i = j;
%--------------------------------------------------------------------------
%Расчёт Y параметров
Y11=j*Wr*(Czi+Czs);
Y12=-j*Wr*Czs;
Y21=S-j*Wr*Czs;
Y22=1/Ri+j*Wr*Czs;
Y=[Y11, Y12;
    Y21, Y22];
%Расчёт Z параметров
Z11=Y22/(det(Y));
Z12=Y12/(det(Y));
Z21=Y21/(det(Y));
Z22=Y11/(det(Y));
Z=[Z11,Z12;
    Z21,Z22];
%--------------------------------------------------------------------------
Xc3=1/(j*Wr*C3);
Xc4=1/(j*Wr*C4);
%коэффициент включения контура
p=L2/L1;
%определяем сопротивление потерь контура
ro=sqrt(L2/C1)*Q; %-характерестическое сопротивление
r=Wr*L2/Q ;%-сопротивление потерь
 
% Находим передаточные характеристики
maxK = 0;
nn = 1;
a1=0;
a=0;
dt = 10000;
i = 1i;
for F = Fr-Fr/2:Fr/dt:Fr+Fr/2;
    w=2*pi*F;
    
    Zn = 1/(1/Rn + i*w*r*C1n);
    
    a1_n=[Z11+R1+Z12,0,0;
        0,Z22+i*Wr*L1+Z21,-i*Wr*M;
        0,-i*Wr*M,i*Wr*L2+r+(1/(i*Wr*C1))+Zn];
    
    a_n=[Z11+R1+Z12,0,1;
        0,Z22+i*Wr*L1+Z21,0;
        0,-i*Wr*M,0];
    
    K_n(nn)=Zn*det(a_n)/det(a1_n);
    K0_n(nn)=abs(Zn*det(a_n)/det(a1_n));
    faza_n(nn)=angle(Zn*det(a_n)/det(a1_n))/(pi)*180;
    if maxK<=K0_n(nn)
        maxK=K0_n(nn);
    end
    if maxK<=K0_n(nn)
        maxN=nn;
    end
    FF(nn)=F;
    nn=nn+1;
end
%--------------------------------------------------------------------------
%находим полосу пропускания
for NN = 1:1:maxN;
    if K0_n(NN)<=0.707*maxK
        f1=NN;
    end
end
NN=0;
for NN = maxN:1:nn-1;
    if K0_n(NN)>=0.707*maxK
        f2=NN;
    end
end
fn=f1;
fw=f2;
fmin07=FF(f1);
fmax07=FF(f2);
%--------------------------------------------------------------------------
% Строим Ачх и Фчх
figure, plot(FF,K0_n);grid;title('АЧХ');
figure, plot(FF,faza_n);grid;title('ФЧХ');
0
0 / 0 / 0
Регистрация: 25.11.2013
Сообщений: 7
27.11.2013, 18:47  [ТС] 7
спасибо хоть на этом)
0
0 / 0 / 0
Регистрация: 25.11.2013
Сообщений: 7
29.11.2013, 00:05  [ТС] 8
а вот у меня ещё 1 вопрос появился
можете посмотреть как преобразовать схему можно и матрицу сопротивлений? а то я пробовал пробовал и нифига не получается уже месяц рис " курсач" по схема экв. транзистора
и да какие элементы надо убрать из схемы? я думаю что С4 и С3(тоесть это С2)
буду очень благодарен очень очень)
Миниатюры
Построение АЧХ и ФЧХ   Построение АЧХ и ФЧХ  
0
IT_Exp
Эксперт
87844 / 49110 / 22898
Регистрация: 17.06.2006
Сообщений: 92,604
29.11.2013, 00:05

Построить АЧХ и ФЧХ фильтра
Построить в Matlab АЧХ и ФЧХ фильтра низких частот с частотой дискретизации равной 1000Гц, частотой...

Нахождение функции АЧХ и ФЧХ
Доброго времени суток! Хотелось бы узнать, можно ли в матлабе вывести функции А(w) и φ(w) для...

Необходимость постройки АЧХ и ФЧХ
Добрый день, помогите пожалуйста с постройкой ачх и фчх и вах Схема и формула во вложении

КИХ-фильтр по заданным АЧХ и ФЧХ
Здравствуйте. Мне нужно сделать КИХ-фильтр нижних частот по заданным АЧХ и ФЧХ аналогового фильтра....


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
8
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin
Copyright ©2000 - 2022, CyberForum.ru