在代码中一些改进:
函数[STR_SDV] = SIS_OP_SDV (K_SDV T1_SDV)
%变量K_SDV T1_SDV
LS_SDV = 3.94 e-6;
DC_SDV = 0.2;
B_SDV = 0.02;
MTTRSD_SDV = 8;
N_SDV = 4;
BS_SDV = B_SDV / 2;
DCS_SDV = LS_SDV。* DC_SDV;
LSU_SDV = (1-DCS_SDV)。* LS_SDV;
LSD_SDV = DCS_SDV。* LS_SDV;
C =(阶乘(abs (N_SDV))) /(阶乘(abs (K_SDV) 1)) * ((1-BS_SDV)。* LSU_SDV + (1-BS_SDV)。* LSD_SDV)。^ (K_SDV);
l = 1:(K_SDV-1);
MDT_SDV_IND = (LSU_SDV. / (LSU_SDV + LSD_SDV)) * ((T1_SDV (:)。/ (l + 1)) + MTTRSD_SDV) + (LSD_SDV. / (LSU_SDV + LSD_SDV)。* MTTRSD_SDV;
PMDT_SDV_IND =刺激(MDT_SDV_IND, 2);
STR_SDV_IND = c * PMDT_SDV_IND;
STR_SDV_CCF = BS_SDV。* LSU_SDV + BS_SDV。* LSD_SDV;
STR_SDV = STR_SDV_IND + STR_SDV_CCF;
STD_SDV =重塑(STR_SDV、大小(T1_SDV));
结束
这些变化,你还必须与标量K_SDV调用该函数,但你可以叫一种基本T1_SDV,结果将是相同的大小,T1_SDV。
你的代码必须使用
整数
K_SDV——你用阶乘()和使用1:K_SDV-1 K_SDV这两个期望整数。1:使用非标量K_SDV K_SDV-1阻碍一部分,因为你使用不同的数组大小不同的K_SDV值就有点讨厌尝试在K_SDV进行向量化。
我再次强调:你
必须
使用整数K_SDV。没有linspace(1、3)为例。
T1_SDV没有限制:您可以使用非整数的值,如果你想要的。linspace(4380、17520)会好参数。
把这些结合在一起,您可以创建一个向量T1_SDV值的第一,然后arrayfun在你
整数
K_SDV值。
这将是一个可怜的轮廓确实如你只会有3个点的一个轴。
你可能会考虑使用γ(abs (K_SDV) 1 + 1)允许非整数K_SDV,但你仍然有问题,l = 1: (K_SDV-1)这意味着你想要一个不同的条件刺激()根据K_SDV,还有没有可行的方法来把条款如果K_SDV非整数的分数。
提示:1:0是空所以当K_SDV是1,1:(K_SDV-1)是1:0空,导致T1_SDV除以空虚(因为空+ 1是空的),导致MDT_SDV_IND是空的,然后刺激(空的)是1。所以当K_SDV是1,PMDT_SDV_IND将为1。当K_SDV 2 1: (K_SDV-1)是1:1 - - > 1所以MDT_SDV_IND将为该产品只有一个词。刺激()不开始变得有意义,直到K_SDV = 3导致1:(3 - 1)——> 1:2
