拟合方程不是次数越高越好。
创新互联专注于遂昌企业网站建设,响应式网站设计,成都做商城网站。遂昌网站建设公司,为遂昌等地区提供建站服务。全流程按需设计,专业设计,全程项目跟踪,创新互联专业和态度为您提供的服务
用分段样条函数也可以。
必要时要删去一些坏的点子。
时序离散数据可以用移动窗修匀。
S域不好编写,你必须离散化传递函数,用Z变换转换为对应的传递函数,然后用c写就可以了。
既然是零阶保持输入,那手算的时候应该也有零阶保持器的吧。
既然离散系统,那么递推关系式也许可以吧。但是自己用C语言写,估计很麻烦,既然有Matlab就用Matlab算啊。
对于含有零阶保持器的离散系统的传递函数你可以上网搜,也可以参考自动控制原理的书籍啊! 至于Matlab使用的方面,实现ABCD可调节也是可以的,你可以搜索Matlab里关于syms或者input函数的使用。 也可以尝试学习GUI的编程。
百度一下
有C++的解题思路和代码。
思路可以借鉴。
离散化:将a数组的备份temp[]排序,然后把不重复的值都弄到X数组中,接下来开始挨着求出a中的每一个元素在X中的位置,用index记录。
这样每次到a[i],index[i]中记录的就是a[i]在X[]中的位置。
尺取法:
通过观察发现,所求序列的第一个一定是在序列中只出现1次的,不然就可以直接把这个舍去了。
设置s,e分别为所求序列的起始和结束。
e每次都++,然后当序列中元素个数==非重复元素个数len时,要用minn记录此时序列长度。然后再s++(直到X[index[s]]==1)。
最后到e不小于n然后结束。
#includecstdio
#include iostream
#includemap
#includeset
#includealgorithm
#include cstring
#define MAXN 1000005
using namespace std;
int n;
int cou;
int a[MAXN];//所有元素
int X[MAXN];//不重复元素
int temp[MAXN];//临时
int inde[MAXN];//存储a[]中每一个元素在X中的下标
int bin_search(int cou,int aa){
int s=0,e=cou-1;
int mid;
while(s=e){
mid=(s+e)1;
if(X[mid]==aa){
return mid;
}else{
if(aaX[mid]){
e=mid-1;
}else{
s=mid+1;
}
}
}
}
void discrete(){
cou=1;
sort(temp,temp+n);
X[0]=temp[0];
for(int i=1;in;i++){
if(temp[i]!=temp[i-1]){
X[cou++]=temp[i];
}
}
for(int i=0;in;i++){
inde[i]=bin_search(cou,a[i]);
}
}
int main()
{
int t;
scanf("%d",t);
while(t--){
int minn=1000005;
scanf("%d",n);
for(int i=0;in;i++){
scanf("%d",a[i]);
temp[i]=a[i];
}
discrete();
memset(X,0,sizeof(X));
int len=cou;
int s=0,e=0;
int number=0;
while(en){
if(X[inde[e]]==0) number++;
X[inde[e]]++;
while(X[inde[s]]=2){
X[inde[s]]--;
s++;
}
if(number==len){
minn=min(minn,e-s+1);
/*if(X[inde[s]]==1)*/ number--;
X[inde[s]]--;
s++;
}
e++;
}
printf("%d\n",minn);
}
return 0;
}
在循环中使用rand函数产生随机数,若想每次运行时产生不同的随机数,在程序开始使用srand函数置一下随机数种子。
传递函数G(s)要离散化成G(z),然后变成差分方程,然后就可以用c语言编程了~
比如Y(S)/U(S)=G(s)=1/(a*s+1);
用零阶保持器,你的采样周期是Ts(比如你的控制函数是在中断服务函数里实现的,那Ts就是你的中断周期),得到Y(z)/U(z)=G(z)=(1-z^(-1))*Z[G(s)/s]=[1-e^(-Ts/a)]/[z-e^(-Ts/a)]={[1-e^(-Ts/a)]*z^(-1)}/[1-e^(-Ts/a)*z^(-1)];
所以得到差分方程:y(k)=e^(-Ts/a)*y(k-1)+ (1-e^(-Ts/a))*u(k-1),k表示当前时刻。