用MATLAB仿真实现双相码给定码元序列为：s=[1 1 1 1 0 0 1 0 1 1 1 0 0 1]的代码函数加脚本

以下是使用MATLAB实现给定码元序列为s=[1 1 1 1 0 0 1 0 1 1 1 0 0 1]的双相码（Biphase Mark Code）的代码和脚本示例：

代码函数：biphase_mark_code.m

function [code] = biphase_mark_code(data)
% 双相码生成函数
% data：输入的二进制数据
% code：生成的双相码序列

code = zeros(1, 2*length(data)); % 初始化双相码序列
prev_bit = 0; % 上一个比特位的值
for i = 1:length(data)
    cur_bit = data(i); % 当前比特位的值
    if cur_bit == 0 % 如果当前比特位为0
        if prev_bit == 0 % 如果上一个比特位也为0
            code(2*i-1) = 1; % 则当前码为正脉冲
            code(2*i) = -1; % 下一个码为负脉冲
        else % 如果上一个比特位为1
            code(2*i-1) = -1; % 则当前码为负脉冲
            code(2*i) = 1; % 下一个码为正脉冲
        end
    else % 如果当前比特位为1
        if prev_bit == 0 % 如果上一个比特位为0
            code(2*i-1) = -1; % 则当前码为负脉冲
            code(2*i) = 1; % 下一个码为正脉冲
        else % 如果上一个比特位也为1
            code(2*i-1) = 1; % 则当前码为正脉冲
            code(2*i) = -1; % 下一个码为负脉冲
        end
    end
    prev_bit = cur_bit; % 记录当前比特位为上一个比特位
end
end

脚本文件：biphase_mark_simulate.m

% 双相码仿真脚本
clc; clear; close all;

% 给定码元序列
s = [1 1 1 1 0 0 1 0 1 1 1 0 0 1];

% 生成双相码
code = biphase_mark_code(s);

% 绘制双相码波形图
t = 0:0.01:length(code)*0.01-0.01;
plot(t, repmat(code, 1, 100)); grid on;
ylim([-1.5, 1.5]); xlabel('时间'); ylabel('幅值'); title('双相码波形图');

使用方法：

1. 在MATLAB中新建一个文件夹，例如biphase_mark_simulation；
2. 将biphase_mark_code.m和biphase_mark_simulate.m文件复制到biphase_mark_simulation文件夹中；
3. 在MATLAB命令窗口中，使用`cd`命令切换到biphase_mark_simulation文件夹中；
4. 运行biphase_mark_simulate.m文件，例如：`biphase_mark_simulate`；
5. MATLAB将生成给定码元序列s的双相码波形图。