写出下列实验代码，2. 利用Matlab或C语言编写交织和解交织程序。3. 利用Matlab或C语言编写和调用信道编码（前面三次实验任选一种）的编码和译码方法。4. 利用Matlab或C语言编写调制解调程序。

由于这是四个不同的实验，我会分别为您提供相关代码。

2. 交织和解交织程序(Matlab)

% 交织和解交织程序
clc; clear all; close all;

% 信号长度
N = 8;

% 交织矩阵
intMat = [1 5 2 6 3 7 4 8];

% 原始信号
x = randn(1,N);

% 交织后的信号
interleaved_x = x(intMat);

% 解交织后的信号
deinterleaved_x = interleaved_x(intMat);

% 显示结果
disp('Original Signal:');
disp(x);
disp('Interleaved Signal:');
disp(interleaved_x);
disp('De-Interleaved Signal:');
disp(deinterleaved_x);

3. 信道编码的编码和译码方法(Matlab)

% 信道编码的编码和译码方法
clc; clear all; close all;

% 生成码率为1/2的卷积码
trellis = poly2trellis(3, [7 5]);

% 随机生成一个信息序列
data = randi([0 1], 1, 10);

% 进行编码
coded_data = convenc(data, trellis);

% 添加噪声，模拟信道传输
received_data = awgn(coded_data, 10, 'measured');

% 进行译码
decoded_data = vitdec(received_data, trellis, 5, 'trunc', 'hard');

% 显示结果
disp('Original Data:');
disp(data);
disp('Coded Data:');
disp(coded_data);
disp('Received Data:');
disp(received_data);
disp('Decoded Data:');
disp(decoded_data);

4. 调制解调程序(Matlab)

% 调制解调程序
clc; clear all; close all;

% 信号长度
N = 1000;

% 频率为1kHz，幅度为1的正弦波
x = sin(2*pi*1000*(0:N-1)/8000);

% BPSK调制
modulated_x = pskmod(x, 2);

% AWGN信道模拟
received_x = awgn(modulated_x, 10, 'measured');

% BPSK解调
demodulated_x = pskdemod(received_x, 2);

% 显示结果
subplot(3,1,1);
plot(x);
title('Original Signal');

subplot(3,1,2);
plot(modulated_x);
title('Modulated Signal');

subplot(3,1,3);
plot(demodulated_x);
title('Demodulated Signal');

以上代码仅供参考，实际应用中可能需要根据具体情况进行修改。同时，C语言代码也可以根据类似的思路实现。