ispgen.m

来自「压缩包中包括1插值与拟合 2规划问题3 绘图 4解方程 5数据分析」· M 代码 · 共 13 行

function [u,df]=ISPgen(e,N,w,lamda)
u=1+fix(unifrnd(0,N));
df=0;
for v=1:N
    if e(v)==1
        df=df+w(v,u);
    end
end
df=1-lamda*df;
if e(u)==1
    df=-df;
end