function [dat,model,na]=datagenerate(id,nois)

k=10;
t=1:30*k;
t1=1:10*k; 
t2=t1+10*k; 
t3=t2+10*k; 
X=[]; Xn=[];
Y=[]; Yn=[]; 

if id==11
X(t)=10+4*sin(0.08*t); 
Y(t)=10+4*cos(0.04*t+4);
Z(t)=2*X(t)+1*Y(t)+2;
model=[0 2 1 0 0];
na=['  True model:  Z = 2X +Y + 2.0'];

elseif id==12 
X(t1)=-0.1*t1+10; 
X(t2)=0.05*t1+0; 
X(t3)=0.05*t1+5;
Y(t)=3+2*sin(0.1*t-3.0);
Z(t)=4*Y(t)+1*X(t).^2-3*X(t).*Y(t)+10;
model=[0 0 4 0 1 0 -3];
na=['  True model:  Z = X^2 - 3XY + 4Y + 10'];

end

Xn= X + randn(size(X))*sqrt(nois);
Yn= Y + randn(size(Y))*sqrt(nois); 
Zn= Z + randn(size(Z))*sqrt(nois);
dat=[Zn;Xn;Yn]';