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

Как построить такой график

28.06.2016, 22:45. Показов 1373. Ответов 12
Метки нет (Все метки)

Доброго времени суток! Подскажите пожалуйста, каким образом мне вот такой график прописать?
Миниатюры
Как построить такой график  
__________________
Помощь в написании контрольных, курсовых и дипломных работ здесь
0
Лучшие ответы (1)
Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
28.06.2016, 22:45
Ответы с готовыми решениями:

Как построить такой график?
Подскажите пожалуйста с помощью каких команд можно строить такие графики:

Как построить такой график ?
Mathlab график

Подскажите, как построить график такой функции
Не получается построить график функции: E=0.06*t^2-0.01172*t-0.066276;

Как построить такой график функции в Matlab?
Совсем не понимаю, как эти промежутки задавать(

12
6024 / 4312 / 1734
Регистрация: 02.02.2014
Сообщений: 11,570
28.06.2016, 23:07 2
варьянт
Кликните здесь для просмотра всего текста
Matlab M
1
2
3
4
5
6
7
8
9
10
clc; clear all; clf;
 
R=10;
for fi=0:pi/12:pi/2
x=R*sin(fi)
y=R*cos(fi)
plot([0,x],[0,y])
hold on
end
axis equal
0
0 / 0 / 0
Регистрация: 09.05.2012
Сообщений: 21
30.06.2016, 19:51  [ТС] 3
Спасибо) Этот вариант подойдет)
Подскажите еще пожалуйста
Делаю два графика в одном и в итоге получается вот так вот
Хотя вертикальные линии должны располагаться так, как показано на первом рисунке
Миниатюры
Как построить такой график  
0
6024 / 4312 / 1734
Регистрация: 02.02.2014
Сообщений: 11,570
06.07.2016, 06:19 4
покажите свой код, что за "50%" на вашем графике?
1
0 / 0 / 0
Регистрация: 09.05.2012
Сообщений: 21
07.07.2016, 23:00  [ТС] 5
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
R=20;
for fi=0:pi/14:pi/2
x=R*sin(fi)
y=R*cos(fi)
plot([0,x],[0,y],'k')
hold on
R1=20;
for fi1=0:pi/114:pi/2
x1=R1*sin(fi1)
y1=R1*cos(fi1)
plot([0,x1],[0,y1],'k--')
hold on
end
end
 
grid on 
text(19.9,2.2,'50%');
text(18.95,6.5,'50%');
text(17.05,10.5,'50%');
text(14.35,13.95,'50%');
text(10.95,16.75,'50%');
text(7.05,18.75,'50%');
text(2.75,19.82,'50%');
hold on
Это то, что вы мне сделали, я лишь только преобразил его немного
Весь код
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
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
clc;clear all; clf;
 
N=4096 %общее число обрабатываемых отсчетов
Fn=35000 %Несущая частота сигнала
Um=1 %Амплитуда сигнала
Dev=10000 %Девиация частоты
Phi0=0 %Начальная фаза сигнала
N3=1 %Число отсчетов задержки сигнала
Pn1=1 %Мощность шума в первом канала
Pn2=1 %Мощность шума во втором канале 
FreqBandOfReceiver=26000 %Полоса пропускания РПУ
Qp=100 %Отношение сигнал/шум по мощности
Fd=165000 %Частота дискретизации
 
N1=4096 %общее число обрабатываемых отсчетов
Fn1=35000 %Несущая частота сигнала
Um1=1 %Амплитуда сигнала
Dev1=10000 %Девиация частоты
Phi0=0 %Начальная фаза сигнала
N31=1 %Число отсчетов задержки сигнала
Pn11=1 %Мощность шума в первом канала
Pn21=1 %Мощность шума во втором канале 
FreqBandOfReceiver=26000 %Полоса пропускания РПУ
Qp=100 %Отношение сигнал/шум по мощности
Fd1=165000 %Частота дискретизации
 
%Вычисления
deltaT=1/Fd %Интервал дискретизации
T3=deltaT*N3 %Время задержки
F=zeros(1,N);
deltaF=1/(N*deltaT)% Вычисление массивов частот
for k=1:N
    F(k)=(k-1)*deltaF;
end
FHalf=zeros(1,N/2);
for k=1:N/2
    FHalf(k)=(k-1)*deltaF;
end
 
Fu=55e3 %Верхняя граничная частота 
Kot=Fd/Fu %Проверка выполнимости условия теоремы Котельникова
Tdescrete=(0:deltaT:(N-1)*deltaT);
%Модулирующая функция
PSP=randn(1,N); %Генерация ПСП
plot(Tdescrete,PSP,'r-') %Изображение модулирующего сигнала
 
%Частотно-модулированный сигнал
Fmsignal=zeros(1,N);
lambda=0;
for k=1:N
    lambda=lambda+PSP(k);
    FMsignal(k)=Um*cos(2*pi*Fn*(deltaT*k)+2*pi*Dev*lambda*deltaT);
end
 
%Изображиние модулированного сигнала
 
plot(Tdescrete, FMsignal, 'r-')
hold off
 
%Вычисления
deltaT1=1/Fd1 %Интервал дискретизации
T3=deltaT1*N31 %Время задержки
F1=zeros(1,N1);
deltaF1=1/(N1*deltaT1)% Вычисление массивов частот
for k1=1:N1
    F1(k1)=(k1-1)*deltaF1;
end
FHalf1=zeros(1,N1/2);
for k1=1:N1/2
    FHalf1(k1)=(k1-1)*deltaF1;
end
 
Fu1=55e3 %Верхняя граничная частота 
Kot=Fd1/Fu1 %Проверка выполнимости условия теоремы Котельникова
Tdescrete1=(0:deltaT1:(N1-1)*deltaT1);
%Модулирующая функция
PSP1=randn(1,N1); %Генерация ПСП
plot(Tdescrete1,PSP1,'r-') %Изображение модулирующего сигнала
 
%Частотно-модулированный сигнал
Fmsignal1=zeros(1,N1);
lambda=0;
for k1=1:N1
    lambda=lambda+PSP1(k1);
    FMsignal1(k1)=Um1*cos(2*pi*Fn1*(deltaT1*k1)+2*pi*Dev1*lambda*deltaT1);
end
 
%Изображиние модулированного сигнала
 
plot(Tdescrete1, FMsignal1, 'b-')
hold off
 
%Задержанный сигнал
FMsignal112=zeros(1,N);
for k=N3+1:N
    FMsignal2(k)=FMsignal(k-N3);
end
for k= 1:N3
    FMsignal2(k)=FMsignal(N-N3+k);
end
 
%Вычисление спектральной рлотности 1 сигнала
 
FMsignalSpectr=fft(FMsignal);
Pyy=FMsignalSpectr.*conj(FMsignalSpectr);
plot(F,Pyy,'r-')
axis([2e4 5e4 0 10e4])
 
%Вычисление спектральной рлотности 2 сигнала
FMsignalSpectr2=fft(FMsignal2);
P2yy=FMsignalSpectr2.*conj(FMsignalSpectr2);
plot(F,P2yy,'r-')
axis([2e4 5e4 0 10e4])
 
%Моделирование шума в первом канале
 
DispersiaNoise1=Pn1 %Дисперсия шума
SigmaNoise1=sqrt(DispersiaNoise1) %СКО шума в первом канале
WhiteNoise1=SigmaNoise1*rand(1,N) %Генерация белого шума в первом канале
 
%Моделирование шума во втором канале
 
DispersiaNoise2=Pn2 %Дисперсия шума
SigmaNoise2=sqrt(DispersiaNoise2) %СКО шума во втором канале
WhiteNoise2=SigmaNoise2*rand(1,N) %Генерация белого шума во втором канале
 
%Смесь сигнала и шума в первом пункте приема
 
Qp1=10000; %Отношение с/ш по мощности
Pc1=Qp1*Pn1 %Мощность сигнала на входе приемника
Qd1=10*log10(Qp1) %отношение с/ш в дБ
Um1=sqrt(2*Pc1) %Множитель амплитуды сигнала на входе РПУ
FMsignalReceiver1=FMsignal.*Um1 %Сигнал на входе РПУ
FMsignalNoise1=FMsignalReceiver1+WhiteNoise1 %смесь сигнала с шумом
 
%Смесь сигнала и шума во втором пункте приема
 
Qp2=10000; %Отношение с/ш по мощности
Pc2=Qp2*Pn2 %Мощность сигнала на входе приемника
Qd2=10*log10(Qp2) %отношение с/ш в дБ
Um2=sqrt(2*Pc2) %Множитель амплитуды сигнала на входе РПУ
FMsignalReceiver2=FMsignal2.*Um2 %Сигнал на входе РПУ
FMsignalNoise2=FMsignalReceiver2+WhiteNoise2 %смесь сигнала с шумом
 
%Вычисление спектральной плотности 1 сигнала
 
FMignalReceiverSpectr1=fft(FMsignalNoise1);
Pyy1=FMignalReceiverSpectr1.*conj(FMignalReceiverSpectr1);
 
%Изображение спектральной плотности 1 сигнала
 
pyy1half=zeros(1,N/2);
for k =1:N/2
    Pyy1Half(k)=Pyy1(k);
end
plot(FHalf,Pyy1Half, 'r-')
 
FMignalReceiverSpectr2=fft(FMsignalNoise2);
Pyy2=FMignalReceiverSpectr2.*conj(FMignalReceiverSpectr2);
 
pyy2half=zeros(1,N/2);
for k =1:N/2
    Pyy2Half(k)=Pyy2(k);
end
plot(FHalf,Pyy2Half, 'r-')
 
%Вычисление взаимного фазового спектра сигналов
 
P12yy=angle(FMignalReceiverSpectr1.*conj(FMignalReceiverSpectr2))
 
P12yyHalf=zeros(1,N/2);
for k=1:N/2
    P12yyHalf(k)=P12yy(k)
end
plot(FHalf,P12yyHalf, 'r-')
hold off
 
Coeff=polyfit(FHalf, P12yyHalf, 1)
Linia=zeros(1,N/2)
for k=1:N/2
    Linia (k)=Coeff(1)*FHalf(k)+Coeff(2);
end
plot(FHalf,Linia, 'r-')
 
hold off
 
%Вычисление задержки
 
Tau=Coeff(1)/2/pi
 
%Задержанный сигнал
FMsignal1121=zeros(1,N1);
for k1=N31+1:N1
    FMsignal21(k1)=FMsignal1(k1-N31);
end
for k1= 1:N31
    FMsignal21(k1)=FMsignal1(N-N3+k1);
end
 
%Вычисление спектральной рлотности 1 сигнала
 
FMsignalSpectr1=fft(FMsignal1);
Pyy11=FMsignalSpectr1.*conj(FMsignalSpectr1);
plot(F1,Pyy11,'r-')
axis([2e4 5e4 0 10e4])
hold off
%Вычисление спектральной рлотности 2 сигнала
FMsignalSpectr21=fft(FMsignal21);
P2yy1=FMsignalSpectr21.*conj(FMsignalSpectr21);
plot(F1,P2yy1,'r-')
axis([2e4 5e4 0 10e4])
hold off
%Моделирование шума в первом канале
 
DispersiaNoise11=Pn11 %Дисперсия шума
SigmaNoise11=sqrt(DispersiaNoise11) %СКО шума в первом канале
WhiteNoise11=SigmaNoise11*rand(1,N1) %Генерация белого шума в первом канале
 
%Моделирование шума во втором канале
 
DispersiaNoise21=Pn21 %Дисперсия шума
SigmaNoise21=sqrt(DispersiaNoise21) %СКО шума во втором канале
WhiteNoise21=SigmaNoise21*rand(1,N1) %Генерация белого шума во втором канале
 
%Смесь сигнала и шума в первом пункте приема
 
Qp11=10000; %Отношение с/ш по мощности
Pc11=Qp11*Pn11 %Мощность сигнала на входе приемника
Qd11=10*log10(Qp11) %отношение с/ш в дБ
Um11=sqrt(2*Pc11) %Множитель амплитуды сигнала на входе РПУ
FMsignalReceiver11=FMsignal1.*Um11 %Сигнал на входе РПУ
FMsignalNoise11=FMsignalReceiver11+WhiteNoise11 %смесь сигнала с шумом
 
%Смесь сигнала и шума во втором пункте приема
 
Qp21=10000; %Отношение с/ш по мощности
Pc21=Qp21*Pn21 %Мощность сигнала на входе приемника
Qd21=10*log10(Qp21) %отношение с/ш в дБ
Um21=sqrt(2*Pc21) %Множитель амплитуды сигнала на входе РПУ
FMsignalReceiver21=FMsignal21.*Um21 %Сигнал на входе РПУ
FMsignalNoise21=FMsignalReceiver21+WhiteNoise21 %смесь сигнала с шумом
 
%Вычисление спектральной плотности 1 сигнала
 
FMignalReceiverSpectr11=fft(FMsignalNoise11);
Pyy11=FMignalReceiverSpectr11.*conj(FMignalReceiverSpectr11);
 
%Изображение спектральной плотности 1 сигнала
 
pyy1half1=zeros(1,N1/2);
for k1 =1:N1/2
    Pyy1Half1(k1)=Pyy11(k1);
end
plot(FHalf1,Pyy1Half1, 'r-')
hold off
 
FMignalReceiverSpectr21=fft(FMsignalNoise21);
Pyy21=FMignalReceiverSpectr21.*conj(FMignalReceiverSpectr21);
 
pyy2half1=zeros(1,N1/2);
for k1 =1:N1/2
    Pyy2Half1(k1)=Pyy21(k1);
end
plot(FHalf1,Pyy2Half1, 'r-')
hold off
 
%Вычисление взаимного фазового спектра сигналов
 
P12yy1=angle(FMignalReceiverSpectr11.*conj(FMignalReceiverSpectr21))
 
P12yyHalf1=zeros(1,N1/2);
for k1=1:N1/2
    P12yyHalf1(k1)=P12yy1(k1)
end
plot(FHalf1,P12yyHalf1, 'r-')
 
Coeff1=polyfit(FHalf1, P12yyHalf1, 1)
Linia1=zeros(1,N1/2)
for k1=1:N/2
    Linia1 (k1)=Coeff1(1)*FHalf1(k1)+Coeff1(2);
end
plot(FHalf,Linia, 'r-',FHalf1,Linia1, 'b-')
hold on
R=20;
for fi=0:pi/14:pi/2
x=R*sin(fi)
y=R*cos(fi)
plot([0,x],[0,y],'k')
hold on
R1=20;
for fi1=0:pi/114:pi/2
x1=R1*sin(fi1)
y1=R1*cos(fi1)
plot([0,x1],[0,y1],'k--')
hold on
end
end
 
grid on 
text(19.9,2.2,'50%');
text(18.95,6.5,'50%');
text(17.05,10.5,'50%');
text(14.35,13.95,'50%');
te
 
xt(10.95,16.75,'50%');
text(7.05,18.75,'50%');
text(2.75,19.82,'50%');
hold on
 
%Вычисление задержки
Tau1=Coeff1(1)/2/pi
0
6024 / 4312 / 1734
Регистрация: 02.02.2014
Сообщений: 11,570
09.07.2016, 08:38 6
можно чуток проще
Кликните здесь для просмотра всего текста
Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
clc; clear all; clf;
 
R=20; R1=18;
for fi=0:pi/18:pi/2
    x=R*sin(fi);
    y=R*cos(fi);
    plot([0,x],[0,y],'k');
    txt=[num2str(fix(fi*57.3)),'%'];
    text(x,y,txt);
    hold on    
end
    for fi1=0:pi/(4*18):pi/2
        x1=R1*sin(fi1);
        y1=R1*cos(fi1);
        plot([0,x1],[0,y1],'k:');
        hold on
    end
grid on
axis equal
0
0 / 0 / 0
Регистрация: 09.05.2012
Сообщений: 21
09.07.2016, 20:01  [ТС] 7
Теперь получается вот так) две линии лежат на горизонтали(синяя и красная), хотя так не должно быть
Миниатюры
Как построить такой график  
0
6024 / 4312 / 1734
Регистрация: 02.02.2014
Сообщений: 11,570
09.07.2016, 20:14 8
Цитата Сообщение от evilred93 Посмотреть сообщение
две линии лежат на горизонтали(синяя и красная), хотя так не должно быть
как написано, так и получилось.. новый код где?
еще линии могут слиться из-за малой величины по отношению к размеру осей..
0
0 / 0 / 0
Регистрация: 09.05.2012
Сообщений: 21
09.07.2016, 23:54  [ТС] 9
Эти линии должны быть вот такими
А вот код
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
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
clc;clear all; clf;
 
N=4096 %общее число обрабатываемых отсчетов
Fn=35000 %Несущая частота сигнала
Um=1 %Амплитуда сигнала
Dev=10000 %Девиация частоты
Phi0=0 %Начальная фаза сигнала
N3=1 %Число отсчетов задержки сигнала
Pn1=1 %Мощность шума в первом канала
Pn2=1 %Мощность шума во втором канале 
FreqBandOfReceiver=26000 %Полоса пропускания РПУ
Qp=100 %Отношение сигнал/шум по мощности
Fd=165000 %Частота дискретизации
 
N1=4096 %общее число обрабатываемых отсчетов
Fn1=50000 %Несущая частота сигнала
Um1=1 %Амплитуда сигнала
Dev1=10000 %Девиация частоты
Phi0=0 %Начальная фаза сигнала
N31=1 %Число отсчетов задержки сигнала
Pn11=1 %Мощность шума в первом канала
Pn21=1 %Мощность шума во втором канале 
FreqBandOfReceiver=26000 %Полоса пропускания РПУ
Qp=100 %Отношение сигнал/шум по мощности
Fd1=195000 %Частота дискретизации
 
 
%Вычисления
deltaT=1/Fd %Интервал дискретизации
T3=deltaT*N3 %Время задержки
F=zeros(1,N);
deltaF=1/(N*deltaT)% Вычисление массивов частот
for k=1:N
    F(k)=(k-1)*deltaF;
end
FHalf=zeros(1,N/2);
for k=1:N/2
    FHalf(k)=(k-1)*deltaF;
end
 
Fu=55000 %Верхняя граничная частота 
Kot=Fd/Fu %Проверка выполнимости условия теоремы Котельникова
Tdescrete=(0:deltaT:(N-1)*deltaT);
%Модулирующая функция
PSP=randn(1,N); %Генерация ПСП
plot(Tdescrete,PSP,'r-') %Изображение модулирующего сигнала
 
%Частотно-модулированный сигнал
Fmsignal=zeros(1,N);
lambda=0;
for k=1:N
    lambda=lambda+PSP(k);
    FMsignal(k)=Um*cos(2*pi*Fn*(deltaT*k)+2*pi*Dev*lambda*deltaT);
end
 
%Изображиние модулированного сигнала
 
plot(Tdescrete, FMsignal, 'r-')
hold off
 
 
%Вычисления
deltaT1=1/Fd1 %Интервал дискретизации
T3=deltaT1*N31 %Время задержки
F1=zeros(1,N1);
deltaF1=1/(N1*deltaT1)% Вычисление массивов частот
for k1=1:N1
    F1(k1)=(k1-1)*deltaF1;
end
FHalf1=zeros(1,N1/2);
for k1=1:N1/2
    FHalf1(k1)=(k1-1)*deltaF1;
end
 
Fu1=65e3 %Верхняя граничная частота 
Kot=Fd1/Fu1 %Проверка выполнимости условия теоремы Котельникова
Tdescrete1=(0:deltaT1:(N1-1)*deltaT1);
%Модулирующая функция
PSP1=randn(1,N1); %Генерация ПСП
plot(Tdescrete1,PSP1,'r-') %Изображение модулирующего сигнала
 
%Частотно-модулированный сигнал
Fmsignal1=zeros(1,N1);
lambda=0;
for k1=1:N1
    lambda=lambda+PSP1(k1);
    FMsignal1(k1)=Um1*cos(2*pi*Fn1*(deltaT1*k1)+2*pi*Dev1*lambda*deltaT1);
end
 
%Изображение модулированного сигнала
 
plot(Tdescrete1, FMsignal1, 'b-')
hold off
 
%Задержанный сигнал
FMsignal112=zeros(1,N);
for k=N3+1:N
    FMsignal2(k)=FMsignal(k-N3);
end
for k= 1:N3
    FMsignal2(k)=FMsignal(N-N3+k);
end
 
%Вычисление спектральной рлотности 1 сигнала
 
FMsignalSpectr=fft(FMsignal);
Pyy=FMsignalSpectr.*conj(FMsignalSpectr);
plot(F,Pyy,'r-')
axis([2e4 5e4 0 10e4])
 
%Вычисление спектральной рлотности 2 сигнала
FMsignalSpectr2=fft(FMsignal2);
P2yy=FMsignalSpectr2.*conj(FMsignalSpectr2);
plot(F,P2yy,'r-')
axis([2e4 5e4 0 10e4])
 
%Моделирование шума в первом канале
 
DispersiaNoise1=Pn1 %Дисперсия шума
SigmaNoise1=sqrt(DispersiaNoise1) %СКО шума в первом канале
WhiteNoise1=SigmaNoise1*rand(1,N) %Генерация белого шума в первом канале
 
%Моделирование шума во втором канале
 
DispersiaNoise2=Pn2 %Дисперсия шума
SigmaNoise2=sqrt(DispersiaNoise2) %СКО шума во втором канале
WhiteNoise2=SigmaNoise2*rand(1,N) %Генерация белого шума во втором канале
 
%Смесь сигнала и шума в первом пункте приема
 
Qp1=10000; %Отношение с/ш по мощности
Pc1=Qp1*Pn1 %Мощность сигнала на входе приемника
Qd1=10*log10(Qp1) %отношение с/ш в дБ
Um1=sqrt(2*Pc1) %Множитель амплитуды сигнала на входе РПУ
FMsignalReceiver1=FMsignal.*Um1 %Сигнал на входе РПУ
FMsignalNoise1=FMsignalReceiver1+WhiteNoise1 %смесь сигнала с шумом
 
%Смесь сигнала и шума во втором пункте приема
 
Qp2=10000; %Отношение с/ш по мощности
Pc2=Qp2*Pn2 %Мощность сигнала на входе приемника
Qd2=10*log10(Qp2) %отношение с/ш в дБ
Um2=sqrt(2*Pc2) %Множитель амплитуды сигнала на входе РПУ
FMsignalReceiver2=FMsignal2.*Um2 %Сигнал на входе РПУ
FMsignalNoise2=FMsignalReceiver2+WhiteNoise2 %смесь сигнала с шумом
 
%Вычисление спектральной плотности 1 сигнала
 
FMignalReceiverSpectr1=fft(FMsignalNoise1);
Pyy1=FMignalReceiverSpectr1.*conj(FMignalReceiverSpectr1);
 
%Изображение спектральной плотности 1 сигнала
 
pyy1half=zeros(1,N/2);
for k =1:N/2
    Pyy1Half(k)=Pyy1(k);
end
plot(FHalf,Pyy1Half, 'r-')
 
FMignalReceiverSpectr2=fft(FMsignalNoise2);
Pyy2=FMignalReceiverSpectr2.*conj(FMignalReceiverSpectr2);
 
pyy2half=zeros(1,N/2);
for k =1:N/2
    Pyy2Half(k)=Pyy2(k);
end
plot(FHalf,Pyy2Half, 'r-')
 
%Вычисление взаимного фазового спектра сигналов
 
P12yy=angle(FMignalReceiverSpectr1.*conj(FMignalReceiverSpectr2))
 
P12yyHalf=zeros(1,N/2);
for k=1:N/2
    P12yyHalf(k)=P12yy(k);
end
plot(FHalf,P12yyHalf, 'r-')
hold off
 
 
Coeff=polyfit(FHalf, P12yyHalf, 1)
Linia=zeros(1,N/2)
for k=1:N/2
    Linia (k)=Coeff(1)*FHalf(k)+Coeff(2);
end
 
%Вычисление задержки
 
Tau=Coeff(1)/2/pi
 
 
%Задержанный сигнал2
FMsignal1121=zeros(1,N1);
for k1=N31+1:N1
    FMsignal21(k1)=FMsignal1(k1-N31);
end
for k1= 1:N31
    FMsignal21(k1)=FMsignal1(N-N3+k1);
end
 
%Вычисление спектральной рлотности 1 сигнала2
 
FMsignalSpectr1=fft(FMsignal1);
Pyy11=FMsignalSpectr1.*conj(FMsignalSpectr1);
plot(F1,Pyy11,'r-')
axis([2e4 5e4 0 10e4])
hold off
%Вычисление спектральной рлотности 2 сигнала2
FMsignalSpectr21=fft(FMsignal21);
P2yy1=FMsignalSpectr21.*conj(FMsignalSpectr21);
plot(F1,P2yy1,'r-')
axis([2e4 5e4 0 10e4])
hold off
%Моделирование шума в первом канале2
 
DispersiaNoise11=Pn11 %Дисперсия шума
SigmaNoise11=sqrt(DispersiaNoise11) %СКО шума в первом канале2
WhiteNoise11=SigmaNoise11*rand(1,N1) %Генерация белого шума в первом канале2
 
%Моделирование шума во втором канале2
 
DispersiaNoise21=Pn21 %Дисперсия шума
SigmaNoise21=sqrt(DispersiaNoise21) %СКО шума во втором канале2
WhiteNoise21=SigmaNoise21*rand(1,N1) %Генерация белого шума во втором канале2
 
%Смесь сигнала и шума в первом пункте приема2
 
Qp11=10000; %Отношение с/ш по мощности2
Pc11=Qp11*Pn11 %Мощность сигнала на входе приемника2
Qd11=10*log10(Qp11) %отношение с/ш в дБ2
Um11=sqrt(2*Pc11) %Множитель амплитуды сигнала на входе РПУ2
FMsignalReceiver11=FMsignal1.*Um11 %Сигнал на входе РПУ2
FMsignalNoise11=FMsignalReceiver11+WhiteNoise11 %смесь сигнала с шумом2
 
%Смесь сигнала и шума во втором пункте приема2
 
Qp21=10000; %Отношение с/ш по мощности2
Pc21=Qp21*Pn21 %Мощность сигнала на входе приемника2
Qd21=10*log10(Qp21) %отношение с/ш в дБ2
Um21=sqrt(2*Pc21) %Множитель амплитуды сигнала на входе РПУ2
FMsignalReceiver21=FMsignal21.*Um21 %Сигнал на входе РПУ2
FMsignalNoise21=FMsignalReceiver21+WhiteNoise21 %смесь сигнала с шумом2
 
%Вычисление спектральной плотности 1 сигнала2
 
FMignalReceiverSpectr11=fft(FMsignalNoise11);
Pyy11=FMignalReceiverSpectr11.*conj(FMignalReceiverSpectr11);
 
%Изображение спектральной плотности 1 сигнала2
 
pyy1half1=zeros(1,N1/2);
for k1 =1:N1/2
    Pyy1Half1(k1)=Pyy11(k1);
end
plot(FHalf1,Pyy1Half1, 'r-')
hold off
 
FMignalReceiverSpectr21=fft(FMsignalNoise21);
Pyy21=FMignalReceiverSpectr21.*conj(FMignalReceiverSpectr21);
 
pyy2half1=zeros(1,N1/2);
for k1 =1:N1/2
    Pyy2Half1(k1)=Pyy21(k1);
end
plot(FHalf1,Pyy2Half1, 'r-')
hold off
 
 
 
%Вычисление взаимного фазового спектра сигналов
 
P12yy1=angle(FMignalReceiverSpectr11.*conj(FMignalReceiverSpectr21))
 
P12yyHalf1=zeros(1,N1/2);
for k1=1:N1/2
    P12yyHalf1(k1)=P12yy1(k1);
end
plot(FHalf1,P12yyHalf1, 'r-')
 
 
 
Coeff1=polyfit(FHalf1, P12yyHalf1, 1)
Linia1=zeros(1,N1/2)
for k1=1:N/2
    Linia1 (k1)=Coeff1(1)*FHalf1(k1)+Coeff1(2);
end
 
 plot(FHalf,Linia, 'r-',FHalf1,Linia1, 'b-')
 axis ([0 3e3 0 25e3]) 
hold on
R=20e3; R1=18e3;
for fi=0:pi/18:pi/2
    x=R*sin(fi);
    y=R*cos(fi);
    plot([0,x],[0,y],'k');
    txt=[num2str(fix(fi*57.3)),'%'];
    text(x,y,txt);
    hold on    
end
    for fi1=0:pi/(4*18):pi/2
        x1=R1*sin(fi1);
        y1=R1*cos(fi1);
        plot([0,x1],[0,y1],'k:');
        hold on
    end
grid on
axis equal
 
%Вычисление задержки
 
Tau1=Coeff1(1)/2/pi
Миниатюры
Как построить такой график  
0
6024 / 4312 / 1734
Регистрация: 02.02.2014
Сообщений: 11,570
10.07.2016, 08:28 10
еще линии могут слиться из-за малой величины по отношению к размеру осей..
про это я и говорю.. шкала оси до 30000, а величина линий всего 3...

Добавлено через 7 минут
в workspace можно посмотреть диапазон величин, у Linia и Linia1 максимальное значение всего 3.05, что вы хотите увидеть на шкале в 25000?
0
0 / 0 / 0
Регистрация: 09.05.2012
Сообщений: 21
10.07.2016, 23:41  [ТС] 11
Да то, что они сливаются, это не важно. Важно, чтобы они не просто на горизонтали были, а как бы по диагонали, как показано на графике сверху. Задача совместить вместе эти графики
0
6024 / 4312 / 1734
Регистрация: 02.02.2014
Сообщений: 11,570
11.07.2016, 06:37 12
Лучший ответ Сообщение было отмечено evilred93 как решение

Решение

внатяжку....
Кликните здесь для просмотра всего текста
Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
 plot(FHalf,Linia, 'r.',FHalf1,Linia1, 'b.')
%  axis ([0 3e3 0 25e3]) 
hold on
R=20; R1=18;
for fi=0:pi/18:pi/2
    x=R*sin(fi)*5e3;
    y=R*cos(fi);
    plot([0,x],[0,y],'k');
    txt=[num2str(fix(fi*57.3)),'%'];
    text(x,y,txt);
    hold on    
end
    for fi1=0:pi/(4*18):pi/2
        x1=R1*sin(fi1)*5e3;
        y1=R1*cos(fi1);
        plot([0,x1],[0,y1],'k:');
        hold on
    end
grid on
%axis equal
1
0 / 0 / 0
Регистрация: 09.05.2012
Сообщений: 21
12.07.2016, 10:25  [ТС] 13
Да, согласен, я туговат здесь)
Спасибо Вам большое)
0
IT_Exp
Эксперт
87844 / 49110 / 22898
Регистрация: 17.06.2006
Сообщений: 92,604
12.07.2016, 10:25

Как создать такой график в MATLAB?

Как построить такой рисунок
Может кто-нибудь сможет помочь в построении рисунка

Как построить график
Здравствуйте, у меня есть массив данных X,Y,Z, по которым нужно построить график, где X,Y, -...

Как построить график u от t
Подскажите пожалуйста, есть цикл который рассчитывает значение V, и соответсвующее значение u(1)...


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

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

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