代码搜索结果

function y = cirshftt(x,m,N) % 长度为 N 的x序列: (时域)作m采样点圆周移位 % ------------------------------------------------------------------- % [y] = cirshftt(x,m,N) % y = 包含圆周移位的输出序列 % x = 长度

function [y] = ovrlpsav(x,h,N) % 用混叠相加法作块卷积 % ---------------------------------------- % [y] = ovrlpsav(x,h,N) % y = 输出序列 % x = 输入序列 % h = 脉冲响应 % N = 块长 % Lenx = length(x); M = length(h); M1

% Q.m % address function function [address_i,in_ani]=Q(in_ani) in_amax=0.7; cap=512; if in_ani>=0.7 in_ani=0.7; address_i=512; elseif in_ani==0 in_ani=0.0001; address_i=1;

function y = up(x,M) % y = up(x,M) % M-fold up-sampling of a 1-D signal [r,c] = size(x); if r > c y = zeros(M*r,1); else y = zeros(1,M*c); end y(1:M:end) = x;

function y = up(x,M) % y = up(x,M) % M-fold up-sampling of a 1-D signal [r,c] = size(x); if r > c y = zeros(M*r,1); else y = zeros(1,M*c); end y(1:M:end) = x;

function y = up(x,M) % y = up(x,M) % M-fold up-sampling of a 1-D signal [r,c] = size(x); if r > c y = zeros(M*r,1); else y = zeros(1,M*c); end y(1:M:end) = x;

%编写惩罚函数 pen.m function y=pen(x) s=0; for i=1:1:length(x) s=s+x(i) end y=min(0,min(x)).^2+(1-s).^2; %编写交易费用函数h.m function y=h(x,u,p) if (x>u) y=x*p else if (x>0) y=u*p

%写一个计算目标函数值M-file myfun.m function F = myfun(x) k = 1:10; F = 2 + 2*k-exp(k*x(1))-exp(k*x(2)); %调用函数lsqnonlin来解此等式 x0 = [0.3 0.4] %初始点 [x,resnorm] = lsqnonlin(@myfun,x0) % 使用优化函数

function rdot=e31a(t,r) c=0.01; g=9.81; m=1; %空气阻力系数及重力加速度（m/s^2） vm=sqrt(r(3)^2+r(4)^2); %速度大小 rdot=[r(3);r(4);-c*vm*r(3)/m;-c*v

function [M,m,df]=fftseq(m,ts,df) % [M,m,df]=fftseq(m,ts,df) % [M,m,df]=fftseq(m,ts) %FFTSEQ generates M, the FFT of the sequence m. % The sequence is zero padded to meet th