matlab制作双声道正弦波wav文件

matlab制作ch=2，bitspersample=24，sample=48k， time=1s 的正弦波wav文件

fs：采样率

sinewave1_fs，正弦波1的频率， 可以自己改

sinewave2_fs，正弦波1的频率， 可以自己改

time，总时长，可以改。

clear; 
close all; 
clc;

Df=1;     %频率间隔
fs = 48000; %采样频率
T = 1/fs;  %采样周期
N=fs/Df    %序列点数
time = (N-1)*T %声音片段的总时长
t=0:T:time;
sinewave1_fs = 10
sinewave2_fs = 20
vol = 0; %声音片段的音量 0db
db2mag_val = (db2mag(vol))
y1 = sin(2*pi*sinewave1_fs*t)*(db2mag(vol)); %生成left ch声音片段，注意需要用db2mag()函数把dB转换成magnitude。
y2 = sin(2*pi*sinewave2_fs*t + pi/2)*(db2mag(vol));%生成right ch声音片段
y = [y1',y2'];
sound(y,fs) %可以播放声音的函数 sound()
filename = ('stereo_sinwave.wav'); %给文件取名
subplot(4,1,1);
plot(t,y1)
subplot(4,1,2);
plot(t,y2)
audiowrite(filename,y,fs,'BitsPerSample',24) %存储.wav音频文件，在这里文件名为sinwave.wav
info = audioinfo(filename)
%%%%%%%%%%%%%%%%%%%%%%%%
[x,FS]=audioread(filename); % 将 WAV 文件转换成变量
x1 = x(:,1); % 抽取第 1 声道
x2 = x(:,2); % 抽取第 2 声道
subplot(4,1,3);
plot(t,x1)
subplot(4,1,4);
plot(t,x2)
%z = dec2hex(x)