質問編集履歴
3
誤字
test
CHANGED
|
@@ -1 +1 @@
|
|
|
1
|
-
MATLABで
|
|
1
|
+
MATLABで復調するコードのエラー
|
test
CHANGED
|
File without changes
|
2
誤字
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -2,15 +2,9 @@
|
|
|
2
2
|
|
|
3
3
|
MATLABでフレーム単位でオーディオデータを取得し、部分的に復調しようとしています。
|
|
4
4
|
|
|
5
|
-
長いオーディオデータを処理するために、並行してデータの取得と復調を行うことが目的です。
|
|
5
|
+
長いオーディオデータを処理するために、並行してデータの取得と復調を行うことが目的です。
|
|
6
6
|
|
|
7
7
|
|
|
8
|
-
|
|
9
|
-
[参考サイト1](https://jp.mathworks.com/help/matlab/matlab_prog/system-objects-vs-matlab-functions.html)
|
|
10
|
-
|
|
11
|
-
[参考サイト2](https://jp.mathworks.com/help/matlab/matlab_prog/system-design-in-matlab-using-system-objects.html)
|
|
12
|
-
|
|
13
|
-
参考サイト2では、「コンポーネントの作成と構成の同時実行」と「システムでのコンポーネントの接続」部分を参考にしています。
|
|
14
8
|
|
|
15
9
|
|
|
16
10
|
|
|
@@ -26,20 +20,6 @@
|
|
|
26
20
|
|
|
27
21
|
audioRead = dsp.AudioFileReader(...
|
|
28
22
|
|
|
29
|
-
'sample.wav',...
|
|
30
|
-
|
|
31
|
-
'SamplesPerFrame',framesize);
|
|
32
|
-
|
|
33
|
-
deviceWriter = audioDeviceWriter(...
|
|
34
|
-
|
|
35
|
-
'SampleRate',Fsa);
|
|
36
|
-
|
|
37
|
-
|
|
38
|
-
|
|
39
|
-
while ~isDone(fileReader)
|
|
40
|
-
|
|
41
|
-
signal = audioRead();
|
|
42
|
-
|
|
43
23
|
|
|
44
24
|
|
|
45
25
|
deviceWriter(s);
|
|
@@ -50,62 +30,6 @@
|
|
|
50
30
|
|
|
51
31
|
|
|
52
32
|
|
|
53
|
-
###該当のソースコード
|
|
54
|
-
|
|
55
|
-
オーディオデータを全て読み込んでから復調するコード（実装済）
|
|
56
|
-
|
|
57
|
-
```Matlab
|
|
58
|
-
|
|
59
|
-
%% 復調
|
|
60
|
-
|
|
61
|
-
Fs = 20000; %audio out sampling
|
|
62
|
-
|
|
63
|
-
Fc = 5000; %搬送波周波数
|
|
64
|
-
|
|
65
|
-
[yy, Fss] = audioread('sample.wav');%ファイルを読む
|
|
66
|
-
|
|
67
|
-
[n,d] = rat(Fss/Fs);%再標本化レート
|
|
68
|
-
|
|
69
|
-
[r,c] = size(yy);
|
|
70
|
-
|
|
71
|
-
|
|
72
|
-
|
|
73
|
-
t = (0:1/Fss:r/Fss-1/Fss);%0~2秒
|
|
74
|
-
|
|
75
|
-
f = (Fss/r:Fss/r:Fss);%周波数成分
|
|
76
|
-
|
|
77
|
-
yc = cos(2*pi*Fc*t);
|
|
78
|
-
|
|
79
|
-
|
|
80
|
-
|
|
81
|
-
z2c = yy' .* yc; %復調
|
|
82
|
-
|
|
83
|
-
Rp = 0.00057565; % 0.01 dB リップル
|
|
84
|
-
|
|
85
|
-
Rst = 1e-6; % 80 dB 阻止域の減衰
|
|
86
|
-
|
|
87
|
-
eqnum = firceqrip(20,Fc/(Fss/2),[Rp Rst],'passedge');
|
|
88
|
-
|
|
89
|
-
|
|
90
|
-
|
|
91
|
-
fvtool(eqnum,'Fs',Fss,'Color', 'White') %フィルターの可視化
|
|
92
|
-
|
|
93
|
-
|
|
94
|
-
|
|
95
|
-
lowpassFIR = dsp.FIRFilter('Numerator', eqnum); %ローパスフィルター（低域通過濾波器）
|
|
96
|
-
|
|
97
|
-
z = lowpassFIR(z2c');
|
|
98
|
-
|
|
99
|
-
fz = resample(z, d, n); %サンプリングレートの変更
|
|
100
|
-
|
|
101
|
-
|
|
102
|
-
|
|
103
|
-
sound(fz, Fs);
|
|
104
|
-
|
|
105
|
-
```
|
|
106
|
-
|
|
107
|
-
|
|
108
|
-
|
|
109
33
|
###補足情報(言語/FW/ツール等のバージョンなど)
|
|
110
34
|
|
|
111
35
|
MatlabR2017b
|
1
誤字脱字の修正
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -10,7 +10,7 @@
|
|
|
10
10
|
|
|
11
11
|
[参考サイト2](https://jp.mathworks.com/help/matlab/matlab_prog/system-design-in-matlab-using-system-objects.html)
|
|
12
12
|
|
|
13
|
-
「コンポーネントの作成と構成の同時実行」と「システムでのコンポーネントの接続」部分を参考にしています。
|
|
13
|
+
参考サイト2では、「コンポーネントの作成と構成の同時実行」と「システムでのコンポーネントの接続」部分を参考にしています。
|
|
14
14
|
|
|
15
15
|
|
|
16
16
|
|