function yi=New_Int(x,y,xi)
n=length(x);
m=length(y);
if n~=m
    error('The lengths of X and Y must be equal');
    return;
end
Y=zeros(n);
Y(:,1)=y';
for k=1:n-1]
    if abs(x(i+k)-x(i))<eps
        error('the DATA is error');
......................

function kj1_KeyPressFcn(hObject, eventdata, handles)
hm=uicontrol(gcf,'style','popupmenu',...
    'string','上下翻转|左右翻转|任意角度翻转',...
    'position',[250 20 50 20],'KeyPressFcn');
set(hm,'value',1)
v=get(hm,'string');
switch v,
        case 1
        axes(handles.axes2);
        x=(handles.img)
        if isrgb(handles.img)
            for k=1:3
......................

z1=80;
z2=80;
z3=60;
for f=1:3
d1=pi/4*2./f;
d2=pi/4*2./f;
d3=pi/3*2./f;
Z0=50;
Z=-j*z3*cot(d3);
Y=1./Z;
Y0=1/z1;
Y1=1/z2;
......................

matlab实现变论域模糊控制出错，三个输入，三个输出
function[sys,x0,str,ts]=bly(t,x,u,flag)

switch flag,
    case 0
        [sys,x0,str,ts]=mdlInitializeSizes; 
    case 3
        sys=mdlOutputs(t,x,u);
    case {1,2,4,9}
    sys=[];
    otherwise
    error(['unhandled flag=',num2str(flag)]);
......................

使用s函数实现变论域模糊控制，以e de n为输入，三个输出，仿真一半时出错
function[sys,x0,str,ts]=bly(t,x,u,flag)

switch flag,
    case 0
        [sys,x0,str,ts]=mdlInitializeSizes; 
    case 3
        sys=mdlOutputs(t,x,u);
    case {1,2,4,9}
    sys=[];
    otherwise
    error(['unhandled flag=',num2str(flag)]);
......................

C = [52 59 43.5 47.8 50 74 80 68 180 200;
        18 180 10 8 4 140 35 47 7 4;
        3 5.9 11 6 14.95 9.35 4.53 18 0.817 0.5;
        0.2 8.25 0.05 0.01 0 9 0.342 0 0 0;
        0.1 22 0.001 0.001 0 20 26.38 4 0 0;
        -5.2 -2.5 0 1.09 -2 29 50.7 5.24 100 100;
        1.806 2.5 5.084 3.142 5.055 2.6 1.7 9.5 1.5 2;
        2.65 2.5 2.71 2.851 2.98 2.67 2.51 2.87 1.03 0.9]
   d = [58.491;96.779;2.004;3.324;7.738;7.934;3.771;2.688]
A=ones(1,10);
A(1,10)=0;A(1,9)=0;
b=1
......................

求各位大神指导MATLAB程序的加密方法

The Merton Jump-Diffusion Model
Merton (1976) [1] jump-diffusion model is an extension to the Geometric Brownian Motion model, with the underlying asset exhibit jumps in addition to have continuous diffusion paths. The asset price evolves as 
dS t =μS t dt + σS t  dW t  + (η−1)dq  
where μ is the drift rate, σ the volatility of S  and dW t   the increment of a Wiener process. The independent Poisson process dq  has a value of 1  with probability λdt  and 0 otherwise. 

The following is a ThetaML implementation of the Merton Jump-Diffusion Model. 
Merton jump diffusion Process simulation
Implementation as ThetaML
model JumpDiffusionProcess
%Jump diffusion model for stock price movement
%Example input parameters:
%S0 = 100, r = 0.03, sigma = 0.2, lambda = 1
......................

求帮忙求一下x1和h的幅频特性曲线，初始值随便取，谢谢啊...
ms=200; 
ks=30000;
mu=40;
Fm=200;
kt=181000;
c0=1000;
B=100;
h=1;
%FM=(Fm*(1-exp(-B*(x(2)-x(4))))/(1+exp(-B*(x(2)-x(4)))));
xr=0.01*sin(2*pi*h*t);

......................

function [L1,XY1,L2,XY2]=JQRLJGH(XX,YY)
%% 基于Dijkstra和遗传算法的机器人路径规划
%  GreenSim团队——专业级算法设计&代写程序
%  欢迎访问GreenSim团队主页→http://blog.sina.com.cn/greensim
%输入参数在函数体内部定义
%输出参数为
%  L1    由Dijkstra算法得出的最短路径长度
%  XY1   由Dijkstra算法得出的最短路径经过节点的坐标
%  L2    由遗传算法得出的最短路径长度
%  XY2   由遗传算法得出的最短路径经过节点的坐标
%程序输出的图片有
%  Fig1  环境地图（包括：边界、障碍物、障碍物顶点之间的连线、Dijkstra的网络图结构）
......................