给蛋蛋看的

已于 2022-01-25 21:18:19 修改
阅读量129 收藏 1
点赞数 1
文章标签： matlab
于 2022-01-12 17:04:44 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Jack_Swan/article/details/122457756
版权

                        

clear all;yalmip('clear');
 
A{1}=[0.1 0.4;-0.8 0.9];B{1}=[1;1];
A{2}=[0.85 1;-0.3 0.9];B{2}=[1;1];
A{3}=[1 0;-0.3 1.3];B{3}=[1;1];
 
[~,l]=size(A);s=2;%被控对象模糊规则数，控制器模糊规则数
[n_x,n_u]=size(B{1});
 
mu_1=[1;0]; 
E=[eye(n_x) zeros(n_x,n_u);zeros(n_u,n_x) zeros(n_u,n_u)];
H=[zeros(n_x,n_u);mu_1(1)*eye(n_u);zeros(n_x,n_u);mu_1(2)*eye(n_u)];
 
P=sdpvar(n_x+n_u,n_x+n_u);
XX1=sdpvar(2*(n_x+n_u),n_x,'full'); 
for j=1:1:s
    G{j}=sdpvar(n_u,n_u);
    GK{j}=sdpvar(n_u,n_x,'full');
end
 
Sigma_d=mdiag(P,-E'*P*E);
for i=1:1:l
    for j=1:1:s
        Xi{i}{j}=Sigma_d...
        +[-XX1 -H*G{j} [XX1*A{i}+H*GK{j} XX1*B{i}-H*G{j}]]...
        +[-XX1 -H*G{j} [XX1*A{i}+H*GK{j} XX1*B{i}-H*G{j}]]';
    end
end
 
 
LMI_S=[lmi(P>=0):'P'];  
for i=1:1:l
    for j=1:1:s
        LMI_S=LMI_S+[lmi(Xi{i}{j}<=0):'LMI_S'];
    end
end
 
options = sdpsettings('solver','mosek','shift','0.00');  
solution = optimize(LMI_S, [],options)
check(LMI_S)
 
for j=1:1:s
    GG{j}=value(G{j})
    GGKK{j}=value(GK{j})
    K{j}=inv(GG{j})*GGKK{j}
end
 
for i=1:1:l
    A{i}
    eig(A{i})
    for j=1:1:s
        AA{i}{j}=A{i}+B{i}*K{j}
        eig(AA{i}{j})
    end
end
 
num=100;x0=[-1;2];
 
XX1=zeros(2,10);
x=x0;
for i=1:1:num
    XX1(:,i)=x;
    %计算隶属函数
    w1_down=(x(1)+2)/10;
    w1_up=(x(1)+3)/10;
    w2_down=(-x(1)+2)/10;
    w2_up=(-x(1)+4)/10;
    w1=(1-(sin(x(1)))^2)*w1_down+((sin(x(1)))^2)*w1_up;
    w2=0.5*w2_down+0.5*w2_up;
    w3=1-w1-w2;
    w{1}=w1;w{2}=w2;w{3}=w3;
    h{1}=w1/(w1+w2);h{2}=w2/(w1+w2);
    K_global=0;
    for j=1:1:s
        K_global=K_global+K{j}*w{j};
%         K_global=K_global+K{j}*h{j};%这应该才是对的，为什么最后不稳定在0处？
    end
    u=K_global*x;
    A_global=0;B_global=0;
    for j=1:1:l
        A_global=A_global+A{j}*w{j};
        B_global=B_global+B{j}*w{j};
    end
    x=A_global*x+B_global*u;
end
XX2=zeros(2,10);
x=x0;
for i=1:1:num
    XX2(:,i)=x;
    %计算隶属函数
    w1_down=(x(1)+2)/10;
    w1_up=(x(1)+3)/10;
    w2_down=(-x(1)+2)/10;
    w2_up=(-x(1)+4)/10;
    w1=(1-(sin(x(1)))^2)*w1_down+((sin(x(1)))^2)*w1_up;
    w2=0.5*w2_down+0.5*w2_up;
    w3=1-w1-w2;
    w{1}=w1;w{2}=w2;w{3}=w3;
    A_global=0;B_global=0;
    for j=1:1:l
        A_global=A_global+A{j}*w{j};
        B_global=B_global+B{j}*w{j};
    end
    x=A_global*x;
end
time=zeros(1,10);
for i=1:1:num
    time(i)=i;
end
 
figure(1)
plot(time, XX1(1,:),'-r', time,XX2(1,:),'-.b',  'linewidth',2);
xlabel('Time in samples');ylabel('x1');
h=legend('$$System\_withController(t)$$','$$System\_withoutController(t)$$');
set(h,'interpreter','latex');
 
figure(2)
plot(time, XX1(2,:),'-r', time,XX2(2,:),'-.b',  'linewidth',2);
xlabel('Time in samples');ylabel('x2');
h=legend('$$System\_withController(t)$$','$$System\_withoutController(t)$$');
set(h,'interpreter','latex');

仅考虑上一时刻的有记忆

clear all;yalmip('clear');

A{1}=[0.1 0.4;-0.8 0.9];B{1}=[1;1];
A{2}=[0.85 1;-0.3 0.9];B{2}=[1;1];
A{3}=[1 0;-0.3 1.3];B{3}=[1;1];

[~,l]=size(A);s=2;%被控对象模糊规则数，控制器模糊规则数
[n_x,n_u]=size(B{1});

mu_1=[1;0]; 
E=[eye(2*n_x) zeros(2*n_x,n_u);zeros(n_u,2*n_x) zeros(n_u,n_u)];
H=[zeros(n_x,n_u);zeros(n_x,n_u);mu_1(1)*eye(n_u);zeros(n_x,n_u);zeros(n_x,n_u);mu_1(2)*eye(n_u)];

P=sdpvar(2*n_x+n_u,2*n_x+n_u);
X1=sdpvar(4*n_x+2*n_u,n_x,'full');
X2=sdpvar(4*n_x+2*n_u,n_x,'full');
for j=1:1:s
    G{j}=sdpvar(n_u,n_u);
    GK1{j}=sdpvar(n_u,n_x,'full');
    GK0{j}=sdpvar(n_u,n_x,'full');
end

a=0.1;
Sigma_d=mdiag(P,(exp(-a)-1)*E'*P*E);

for i=1:1:l
    for j=1:1:s
        Xi{i}{j}=Sigma_d...
        +[-X1 -X2 -H*G{j} [H*GK1{j} X1+X2*A{i}+H*GK0{j} X2*B{i}-H*G{j}]]...
        +[-X1 -X2 -H*G{j} [H*GK1{j} X1+X2*A{i}+H*GK0{j} X2*B{i}-H*G{j}]]';
    end
end


LMI_S=[lmi(P>=0):'P'];  
for i=1:1:l
    for j=1:1:s
        LMI_S=LMI_S+[lmi(Xi{i}{j}<=0):'LMI_S'];
    end
end
 
options = sdpsettings('solver','mosek','shift','0.00');  
solution = optimize(LMI_S, [],options)
check(LMI_S)

for j=1:1:s
    GG{j}=value(G{j})
    GGKK1{j}=value(GK1{j})
    GGKK0{j}=value(GK0{j})
    K1{j}=inv(GG{j})*GGKK1{j}
    K0{j}=inv(GG{j})*GGKK0{j}
end

for i=1:1:l
    A{i}
    eig(A{i})
    for j=1:1:s
        AA{i}{j}=A{i}+B{i}*K1{j}+B{i}*K0{j}
        eig(AA{i}{j})
    end
end

num=100;x0=[-1;2];

XX1=zeros(2,10);
x=x0;
for i=1:1:num
    XX1(:,i)=x;
    %计算隶属函数
    w1_down=(x(1)+2)/10;
    w1_up=(x(1)+3)/10;
    w2_down=(-x(1)+2)/10;
    w2_up=(-x(1)+4)/10;
    w1=(1-(sin(x(1)))^2)*w1_down+((sin(x(1)))^2)*w1_up;
    w2=0.5*w2_down+0.5*w2_up;
    w3=1-w1-w2;
    w{1}=w1;w{2}=w2;w{3}=w3;
    K1_global=0;K0_global=0;
    for j=1:1:s
        K1_global=K1_global+K1{j}*w{j};
        K0_global=K0_global+K0{j}*w{j};
    end
    if i==1
        u=K0_global*x;
    else
        u=K1_global*XX1(:,i-1)+K0_global*x;
    end
    A_global=0;B_global=0;
    for j=1:1:l
        A_global=A_global+A{j}*w{j};
        B_global=B_global+B{j}*w{j};
    end
    x=A_global*x+B_global*u;
end
XX2=zeros(2,10);
x=x0;
for i=1:1:num
    XX2(:,i)=x;
    %计算隶属函数
    w1_down=(x(1)+2)/10;
    w1_up=(x(1)+3)/10;
    w2_down=(-x(1)+2)/10;
    w2_up=(-x(1)+4)/10;
    w1=(1-(sin(x(1)))^2)*w1_down+((sin(x(1)))^2)*w1_up;
    w2=0.5*w2_down+0.5*w2_up;
    w3=1-w1-w2;
    w{1}=w1;w{2}=w2;w{3}=w3;
    A_global=0;B_global=0;
    for j=1:1:l
        A_global=A_global+A{j}*w{j};
        B_global=B_global+B{j}*w{j};
    end
    x=A_global*x;
end
time=zeros(1,10);
for i=1:1:num
    time(i)=i;
end

figure(1)
plot(time, XX1(1,:),'-r', time,XX2(1,:),'-.b',  'linewidth',2);
xlabel('Time in samples');ylabel('x1');
h=legend('$$System\_withController(t)$$','$$System\_withoutController(t)$$');
set(h,'interpreter','latex');

figure(2)
plot(time, XX1(2,:),'-r', time,XX2(2,:),'-.b',  'linewidth',2);
xlabel('Time in samples');ylabel('x2');
h=legend('$$System\_withController(t)$$','$$System\_withoutController(t)$$');
set(h,'interpreter','latex');

仅考虑上一时刻的有记忆输出反馈

clear all;yalmip('clear');

A{1}=[0.1 0.4;-0.8 0.9];B{1}=[1;1];C1{1}=[1 1];
A{2}=[0.85 1;-0.3 0.9];B{2}=[1;1];C1{2}=[1 1];
A{3}=[1 0;-0.3 1.3];B{3}=[1;1];C1{3}=[1 1];

[~,l]=size(A);s=2;%被控对象模糊规则数，控制器模糊规则数
[n_x,n_u]=size(B{1});
[n_y,~]=size(C1{1});

mu_1=[1;0]; 
E=[eye(n_y+n_x) zeros(n_y+n_x,n_u);zeros(n_u,n_y+n_x) zeros(n_u,n_u)];
H=[zeros(n_y,n_u);zeros(n_x,n_u);mu_1(1)*eye(n_u);zeros(n_y,n_u);zeros(n_x,n_u);mu_1(2)*eye(n_u)];

P=sdpvar(n_y+n_x+n_u,n_y+n_x+n_u);
X1=sdpvar(2*(n_y+n_x+n_u),n_y,'full');
X2=sdpvar(2*(n_y+n_x+n_u),n_x,'full');
for j=1:1:s
    G{j}=sdpvar(n_u,n_u);
    GK1{j}=sdpvar(n_u,n_y,'full');
    GK0{j}=sdpvar(n_u,n_y,'full');
end

a=0.1;
Sigma_d=mdiag(P,(exp(-a)-1)*E'*P*E);

for i=1:1:l
    for j=1:1:s
        Xi{i}{j}=Sigma_d...
        +[-X1 -X2 -H*G{j} [H*GK1{j} X1*C1{i}+X2*A{i}+H*GK0{j}*C1{i} X2*B{i}-H*G{j}]]...
        +[-X1 -X2 -H*G{j} [H*GK1{j} X1*C1{i}+X2*A{i}+H*GK0{j}*C1{i} X2*B{i}-H*G{j}]]';
    end
end


LMI_S=[lmi(P>=0):'P'];  
for i=1:1:l
    for j=1:1:s
        LMI_S=LMI_S+[lmi(Xi{i}{j}<=0):'LMI_S'];
    end
end
 
options = sdpsettings('solver','mosek','shift','0.00');  
solution = optimize(LMI_S, [],options)
check(LMI_S)

for j=1:1:s
    GG{j}=value(G{j})
    GGKK1{j}=value(GK1{j})
    GGKK0{j}=value(GK0{j})
    K1{j}=inv(GG{j})*GGKK1{j}
    K0{j}=inv(GG{j})*GGKK0{j}
end

for i=1:1:l
    A{i}
    eig(A{i})
    for j=1:1:s
        AA{i}{j}=A{i}+B{i}*K1{j}*C1{i}+B{i}*K0{j}*C1{i}
        eig(AA{i}{j})
    end
end

num=100;x0=[-1;2];
u=0;
U=zeros(1,10);
XX1=zeros(2,10);
x=x0;
for i=1:1:num
    XX1(:,i)=x;
    %计算隶属函数
    w1_down=(x(1)+2)/10;
    w1_up=(x(1)+3)/10;
    w2_down=(-x(1)+2)/10;
    w2_up=(-x(1)+4)/10;
    w1=(1-(sin(x(1)))^2)*w1_down+((sin(x(1)))^2)*w1_up;
    w2=0.5*w2_down+0.5*w2_up;
    w3=1-w1-w2;
    w{1}=w1;w{2}=w2;w{3}=w3;
    K1_global=0;K0_global=0;
    for j=1:1:s
        for k=1:1:l
            K1_global=K1_global+K1{j}*C1{k}*w{j};
            K0_global=K0_global+K0{j}*C1{k}*w{j};
        end
    end
    if i==1
        u=K0_global*x;
    else
        u=K1_global*XX1(:,i-1)+K0_global*x;
    end
    A_global=0;B_global=0;
    for j=1:1:l
        A_global=A_global+A{j}*w{j};
        B_global=B_global+B{j}*w{j};
    end
    x=A_global*x+B_global*u;
    U(:,i)=u;
end
XX2=zeros(2,10);
x=x0;
for i=1:1:num
    XX2(:,i)=x;
    %计算隶属函数
    w1_down=(x(1)+2)/10;
    w1_up=(x(1)+3)/10;
    w2_down=(-x(1)+2)/10;
    w2_up=(-x(1)+4)/10;
    w1=(1-(sin(x(1)))^2)*w1_down+((sin(x(1)))^2)*w1_up;
    w2=0.5*w2_down+0.5*w2_up;
    w3=1-w1-w2;
    w{1}=w1;w{2}=w2;w{3}=w3;
    A_global=0;B_global=0;
    for j=1:1:l
        A_global=A_global+A{j}*w{j};
        B_global=B_global+B{j}*w{j};
    end
    x=A_global*x;
end
time=zeros(1,10);
for i=1:1:num
    time(i)=i;
end

figure()
plot(time, U, '-r' , 'linewidth',2);
xlabel('Time in samples');ylabel('u');

figure()
plot(time, XX1(1,:),'-r', time,XX2(1,:),'-.b',  'linewidth',2);
xlabel('Time in samples');ylabel('x1');
h=legend('$$System\_withController(t)$$','$$System\_withoutController(t)$$');
set(h,'interpreter','latex');

figure()
plot(time, XX1(2,:),'-r', time,XX2(2,:),'-.b',  'linewidth',2);
xlabel('Time in samples');ylabel('x2');
h=legend('$$System\_withController(t)$$','$$System\_withoutController(t)$$');
set(h,'interpreter','latex');

多层历史数据的，存在问题在改。

clear all;yalmip('clear');

A{1}=[0.1 0.4;-0.8 0.9];B{1}=[1;1];C1{1}=[1 1];
A{2}=[0.85 1;-0.3 0.9];B{2}=[1;1];C1{2}=[1 1];
A{3}=[1 0;-0.3 1.3];B{3}=[1;1];C1{3}=[1 1];

[~,l]=size(A);s=2;%被控对象模糊规则数，控制器模糊规则数
[n_x,n_u]=size(B{1});
[n_y,~]=size(C1{1});

lambda=1;%历史数据层数

mu_1=[1;0]; 
E=[eye(lambda*n_y+n_x) zeros(lambda*n_y+n_x,n_u);zeros(n_u,lambda*n_y+n_x) zeros(n_u,n_u)];
H=[zeros(lambda*n_y,n_u);zeros(n_x,n_u);mu_1(1)*eye(n_u);zeros(lambda*n_y,n_u);zeros(n_x,n_u);mu_1(2)*eye(n_u)];

P=sdpvar(lambda*n_y+n_x+n_u,lambda*n_y+n_x+n_u);
X1=sdpvar(2*(lambda*n_y+n_x+n_u),lambda*n_y,'full');
X2=sdpvar(2*(lambda*n_y+n_x+n_u),n_x,'full');
for j=1:1:s
    G{j}=sdpvar(n_u,n_u);
    for n=1:1:lambda+1
        GK{n}{j}=sdpvar(n_u,n_y,'full');
    end
%     GK1{j}=sdpvar(n_u,n_y,'full');
%     GK0{j}=sdpvar(n_u,n_y,'full');
end

a=0.1;
Sigma_d=mdiag(P,(exp(-a)-1)*E'*P*E);

for i=1:1:l
    for j=1:1:s
        GK_M=[];
        for n=1:1:lambda
            GK_M=[GK_M;GK{n}{j}];
        end
        Xi{i}{j}=Sigma_d...
        +[-X1 -X2 -H*G{j} [H*GK_M X1*C1{i}+X2*A{i}+H*GK{lambda+1}{j}*C1{i} X2*B{i}-H*G{j}]]...
        +[-X1 -X2 -H*G{j} [H*GK_M X1*C1{i}+X2*A{i}+H*GK{lambda+1}{j}*C1{i} X2*B{i}-H*G{j}]]';
    end
end


LMI_S=[lmi(P>=0):'P'];  
for i=1:1:l
    for j=1:1:s
        LMI_S=LMI_S+[lmi(Xi{i}{j}<=0):'LMI_S'];
    end
end
 
options = sdpsettings('solver','mosek','shift','0.00');  
solution = optimize(LMI_S, [],options)
check(LMI_S)

for j=1:1:s
    GG{j}=value(G{j})
    for n=1:1:lambda+1
        GGKK{n}{j}=value(GK{n}{j});
        K{n}{j}=inv(GG{j})*GGKK{n}{j}
    end
end

for i=1:1:l
    A{i}
    eig(A{i})
    for j=1:1:s
        for n=1:1:lambda+1
            AA{i}{j}=A{i}+B{i}*K{n}{j}*C1{i}
        end
        eig(AA{i}{j})
    end
end

num=100;x0=[-1;2];
u=0;
U=zeros(1,10);
XX1=zeros(2,10);
x=x0;
for i=1:1:num
    XX1(:,i)=x;
    %计算隶属函数
    w1_down=(x(1)+2)/10;
    w1_up=(x(1)+3)/10;
    w2_down=(-x(1)+2)/10;
    w2_up=(-x(1)+4)/10;
    w1=(1-(sin(x(1)))^2)*w1_down+((sin(x(1)))^2)*w1_up;
    w2=0.5*w2_down+0.5*w2_up;
    w3=1-w1-w2;
    w{1}=w1;w{2}=w2;w{3}=w3;
%     K1_global=0;K0_global=0;
    for n=1:1:lambda+1
        K_global{n}=0;
        k_global=0;
        for j=1:1:s
            for k=1:1:l
                k_global=k_global+K{n}{j}*C1{k}*w{j};%这里隶属函数有问题，但是算了，后面再说。
            end
        end
        K_global{n}=k_global;
    end
    if i<=lambda
        u=0;
        for n=lambda+1:1:lambda+2-i
            u=u+K_global{n}*XX1(:,i-(lambda+1-n));
        end
    else
        u=0;
        for n=lambda+1:1:1
            u=u+K_global{n}*XX1(:,i-(lambda+1-n));
        end
    end
    A_global=0;B_global=0;
    for j=1:1:l
        A_global=A_global+A{j}*w{j};
        B_global=B_global+B{j}*w{j};
    end
    x=A_global*x+B_global*u;
    U(:,i)=u;
end
XX2=zeros(2,10);
x=x0;
for i=1:1:num
    XX2(:,i)=x;
    %计算隶属函数
    w1_down=(x(1)+2)/10;
    w1_up=(x(1)+3)/10;
    w2_down=(-x(1)+2)/10;
    w2_up=(-x(1)+4)/10;
    w1=(1-(sin(x(1)))^2)*w1_down+((sin(x(1)))^2)*w1_up;
    w2=0.5*w2_down+0.5*w2_up;
    w3=1-w1-w2;
    w{1}=w1;w{2}=w2;w{3}=w3;
    A_global=0;B_global=0;
    for j=1:1:l
        A_global=A_global+A{j}*w{j};
        B_global=B_global+B{j}*w{j};
    end
    x=A_global*x;
end
time=zeros(1,10);
for i=1:1:num
    time(i)=i;
end

figure()
plot(time, U, '-r' , 'linewidth',2);
xlabel('Time in samples');ylabel('u');

figure()
plot(time, XX1(1,:),'-r', time,XX2(1,:),'-.b',  'linewidth',2);
xlabel('Time in samples');ylabel('x1');
h=legend('$$System\_withController(t)$$','$$System\_withoutController(t)$$');
set(h,'interpreter','latex');

figure()
plot(time, XX1(2,:),'-r', time,XX2(2,:),'-.b',  'linewidth',2);
xlabel('Time in samples');ylabel('x2');
h=legend('$$System\_withController(t)$$','$$System\_withoutController(t)$$');
set(h,'interpreter','latex');

信道衰落的一次历史数据有记忆输出反馈控制

clear all;yalmip('clear');

A{1}=[0.1 0.4;-0.8 0.9];B{1}=[1;1];C1{1}=[1 1];
A{2}=[0.85 1;-0.3 0.9];B{2}=[1;1];C1{2}=[1 1];
A{3}=[1 0;-0.3 1.3];B{3}=[1;1];C1{3}=[1 1];

[~,l]    =size(A);s=2;%被控对象模糊规则数，控制器模糊规则数
[n_x,n_u]=size(B{1});
[n_y,~]  =size(C1{1});

mu1=0.8;mu2=0.05;

mu_1=[1;1;0];mu_2=[1;1;0]; 
E   =[eye(n_y+n_x) zeros(n_y+n_x,n_u);zeros(n_u,n_y+n_x) zeros(n_u,n_u)];
H1  =[zeros(n_y,n_u);zeros(n_x,n_u);mu_1(1)*eye(n_u);zeros(n_y,n_u);zeros(n_x,n_u);mu_1(2)*eye(n_u);zeros(n_y,n_u);zeros(n_x,n_u);mu_1(3)*eye(n_u)];
H2  =[zeros(n_y,n_u);zeros(n_x,n_u);mu_2(1)*eye(n_u);zeros(n_y,n_u);zeros(n_x,n_u);mu_2(2)*eye(n_u);zeros(n_y,n_u);zeros(n_x,n_u);mu_2(3)*eye(n_u)];

P  =sdpvar(n_y+n_x+n_u,n_y+n_x+n_u);
X11=sdpvar(3*(n_y+n_x+n_u),n_y,'full');
X12=sdpvar(3*(n_y+n_x+n_u),n_x,'full');
X21=sdpvar(3*(n_y+n_x+n_u),n_y,'full');
X22=sdpvar(3*(n_y+n_x+n_u),n_x,'full');
for j=1:1:s
    G{j}  =sdpvar(n_u,n_u);
    GK1{j}=sdpvar(n_u,n_y,'full');
    GK0{j}=sdpvar(n_u,n_y,'full');
end   

a=0.1;
Sigma_d=mdiag(P,P,(exp(-a)-1)*E'*P*E);

for i=1:1:l
    for j=1:1:s
        X1A=[mu1*H1*GK1{j} X11*C1{i}+X12*A{i}+mu1*H1*GK0{j}*C1{i} X12*B{i}-H1*G{j}];
        X2A=[sqrt(mu2)*H2*GK1{j} sqrt(mu2)*H2*GK0{j}*C1{i} zeros(3*(n_y+n_x+n_u),n_u)];
        Xi{i}{j}=Sigma_d...
            +[-X11 -X12 -H1*G{j} -X21 -X22 -H2*G{j} X1A+X2A]...
            +[-X11 -X12 -H1*G{j} -X21 -X22 -H2*G{j} X1A+X2A]';
    end
end


LMI_S=[lmi(P>=0):'P'];  
for i=1:1:l
    for j=1:1:s
        LMI_S=LMI_S+[lmi(Xi{i}{j}<=0):'LMI_S'];
    end
end
 
options = sdpsettings('solver','mosek','shift','0.00');  
solution = optimize(LMI_S, [],options)
check(LMI_S)

for j=1:1:s
    GG{j}=value(G{j})
    GGKK1{j}=value(GK1{j})
    GGKK0{j}=value(GK0{j})
    K1{j}=inv(GG{j})*GGKK1{j}
    K0{j}=inv(GG{j})*GGKK0{j}
end

for i=1:1:l
    A{i}
    eig(A{i})
    for j=1:1:s
        AA{i}{j}=A{i}+B{i}*K1{j}*C1{i}+B{i}*K0{j}*C1{i}
        eig(AA{i}{j})
    end
end

num=100;x0=[-1;2];
u=0;
U=zeros(1,10);
XX1=zeros(2,10);
x=x0;
for i=1:1:num
    XX1(:,i)=x;
    %计算隶属函数
    w1_down=(x(1)+2)/10;
    w1_up=(x(1)+3)/10;
    w2_down=(-x(1)+2)/10;
    w2_up=(-x(1)+4)/10;
    w1=(1-(sin(x(1)))^2)*w1_down+((sin(x(1)))^2)*w1_up;
    w2=0.5*w2_down+0.5*w2_up;
    w3=1-w1-w2;
    w{1}=w1;w{2}=w2;w{3}=w3;
    K1_global=0;K0_global=0;
    for j=1:1:s
        for k=1:1:l
            K1_global=K1_global+K1{j}*C1{k}*w{j};
            K0_global=K0_global+K0{j}*C1{k}*w{j};
        end
    end
    if i==1
        u=K0_global*x;
    else
        u=K1_global*XX1(:,i-1)+K0_global*x;
    end
    A_global=0;B_global=0;
    for j=1:1:l
        A_global=A_global+A{j}*w{j};
        B_global=B_global+B{j}*w{j};
    end
    x=A_global*x+B_global*u;
    U(:,i)=u;
end
XX2=zeros(2,10);
x=x0;
for i=1:1:num
    XX2(:,i)=x;
    %计算隶属函数
    w1_down=(x(1)+2)/10;
    w1_up=(x(1)+3)/10;
    w2_down=(-x(1)+2)/10;
    w2_up=(-x(1)+4)/10;
    w1=(1-(sin(x(1)))^2)*w1_down+((sin(x(1)))^2)*w1_up;
    w2=0.5*w2_down+0.5*w2_up;
    w3=1-w1-w2;
    w{1}=w1;w{2}=w2;w{3}=w3;
    A_global=0;B_global=0;
    for j=1:1:l
        A_global=A_global+A{j}*w{j};
        B_global=B_global+B{j}*w{j};
    end
    x=A_global*x;
end
time=zeros(1,10);
for i=1:1:num
    time(i)=i;
end

figure()
plot(time, U, '-r' , 'linewidth',2);
xlabel('Time in samples');ylabel('u');

figure()
plot(time, XX1(1,:),'-r', time,XX2(1,:),'-.b',  'linewidth',2);
xlabel('Time in samples');ylabel('x1');
h=legend('$$System\_withController(t)$$','$$System\_withoutController(t)$$');
set(h,'interpreter','latex');

figure()
plot(time, XX1(2,:),'-r', time,XX2(2,:),'-.b',  'linewidth',2);
xlabel('Time in samples');ylabel('x2');
h=legend('$$System\_withController(t)$$','$$System\_withoutController(t)$$');
set(h,'interpreter','latex');
%我不理解，为什么会出现非线性，一个转折？

在上一个的基础上考虑扰动

clear all;yalmip('clear');

A{1}=[0.1 0.4;-0.8 0.9];B{1}=[1;1];D1{1}=[0.1;0.1];
C1{1}=[1 1];
C2{1}=[1 1];D2{1}=[0.1];
A{2}=[0.85 1;-0.3 0.9];B{2}=[1;1];D1{2}=[0.1;0.1];
C1{2}=[1 1];
C2{2}=[1 1];D2{2}=[0.1];
A{3}=[1 0;-0.3 1.3];B{3}=[1;1];D1{3}=[0.1;0.1];
C1{3}=[1 1];
C2{3}=[1 1];D2{3}=[0.1];

[~,l]    =size(A);s=2;%被控对象模糊规则数，控制器模糊规则数
[n_x,n_u]=size(B{1});
[n_y,~]  =size(C1{1});
[n_z,n_d]=size(D2{1});

mu1=0.8;mu2=0.05;

mu_1=[1;1;0];
mu_2=[1;1;0]; 
E   =[eye(n_y+n_x) zeros(n_y+n_x,n_u);zeros(n_u,n_y+n_x) zeros(n_u,n_u)];
H1  =[zeros(n_y,n_u);zeros(n_x,n_u);mu_1(1)*eye(n_u);...
      zeros(n_y,n_u);zeros(n_x,n_u);mu_1(2)*eye(n_u);...
      zeros(n_y,n_u);zeros(n_x,n_u);mu_1(3)*eye(n_u);zeros(n_d,n_u)];
H2  =[zeros(n_y,n_u);zeros(n_x,n_u);mu_2(1)*eye(n_u);...
      zeros(n_y,n_u);zeros(n_x,n_u);mu_2(2)*eye(n_u);...
      zeros(n_y,n_u);zeros(n_x,n_u);mu_2(3)*eye(n_u);zeros(n_d,n_u)];

P  =sdpvar(n_y+n_x+n_u,n_y+n_x+n_u);
X11=sdpvar(3*(n_y+n_x+n_u)+n_d,n_y,'full');
X12=sdpvar(3*(n_y+n_x+n_u)+n_d,n_x,'full');
X21=sdpvar(3*(n_y+n_x+n_u)+n_d,n_y,'full');
X22=sdpvar(3*(n_y+n_x+n_u)+n_d,n_x,'full');
GAMMA2=sdpvar(1,1,'full');
for j=1:1:s
    G{j}  =sdpvar(n_u,n_u);
    GK1{j}=sdpvar(n_u,n_y,'full');
    GK0{j}=sdpvar(n_u,n_y,'full');
end   

a=0.1;

for i=1:1:l
    C2Bar=[zeros(n_z,n_y) C2{i} zeros(n_z,n_u)];
    Sigma_d=mdiag(P,P,[(exp(-a)-1)*E'*P*E+C2Bar'*C2Bar C2Bar'*D2{i};D2{i}'*C2Bar D2{i}'*D2{i}-GAMMA2*eye(n_d)]);
    for j=1:1:s
        X1A=[mu1*H1*GK1{j} X11*C1{i}+X12*A{i}+mu1*H1*GK0{j}*C1{i} X12*B{i}-H1*G{j}];
        X2A=[sqrt(mu2)*H2*GK1{j} sqrt(mu2)*H2*GK0{j}*C1{i} zeros(3*(n_y+n_x+n_u)+n_d,n_u)];
        Xi{i}{j}=Sigma_d...
            +[-X11 -X12 -H1*G{j} -X21 -X22 -H2*G{j} X1A+X2A X12*D1{i}]...
            +[-X11 -X12 -H1*G{j} -X21 -X22 -H2*G{j} X1A+X2A X12*D1{i}]';
    end
end

LMI_S=[lmi(P>=0):'P'];  
for i=1:1:l
    for j=1:1:s
        LMI_S=LMI_S+[lmi(Xi{i}{j}<=0):'LMI_S'];
    end
end
 
options = sdpsettings('solver','mosek','shift','0.00');  
solution = optimize(LMI_S, [],options)
check(LMI_S)

gamma=sqrt(value(GAMMA2))
for j=1:1:s
    GG{j}=value(G{j})
    GGKK1{j}=value(GK1{j})
    GGKK0{j}=value(GK0{j})
    K1{j}=inv(GG{j})*GGKK1{j}
    K0{j}=inv(GG{j})*GGKK0{j}
end

for i=1:1:l
    A{i}
    eig(A{i})
    for j=1:1:s
        AA{i}{j}=A{i}+B{i}*K1{j}*C1{i}+B{i}*K0{j}*C1{i}
        eig(AA{i}{j})
    end
end

num=100;x0=[-1;2];
u=0;
U=zeros(1,10);
XX1=zeros(2,10);
x=x0;
for i=1:1:num
    XX1(:,i)=x;
    %计算隶属函数
    w1_down=(x(1)+2)/10;
    w1_up=(x(1)+3)/10;
    w2_down=(-x(1)+2)/10;
    w2_up=(-x(1)+4)/10;
    w1=(1-(sin(x(1)))^2)*w1_down+((sin(x(1)))^2)*w1_up;
    w2=0.5*w2_down+0.5*w2_up;
    w3=1-w1-w2;
    w{1}=w1;w{2}=w2;w{3}=w3;
    K1_global=0;K0_global=0;
    for j=1:1:s
        for k=1:1:l
            K1_global=K1_global+K1{j}*C1{k}*w{j};
            K0_global=K0_global+K0{j}*C1{k}*w{j};
        end
    end
    if i==1
        u=K0_global*x;
    else
        u=K1_global*XX1(:,i-1)+K0_global*x;
    end
    A_global=0;B_global=0;
    for j=1:1:l
        A_global=A_global+A{j}*w{j};
        B_global=B_global+B{j}*w{j};
    end
    x=A_global*x+B_global*u;
    U(:,i)=u;
end
XX2=zeros(2,10);
x=x0;
for i=1:1:num
    XX2(:,i)=x;
    %计算隶属函数
    w1_down=(x(1)+2)/10;
    w1_up=(x(1)+3)/10;
    w2_down=(-x(1)+2)/10;
    w2_up=(-x(1)+4)/10;
    w1=(1-(sin(x(1)))^2)*w1_down+((sin(x(1)))^2)*w1_up;
    w2=0.5*w2_down+0.5*w2_up;
    w3=1-w1-w2;
    w{1}=w1;w{2}=w2;w{3}=w3;
    A_global=0;B_global=0;
    for j=1:1:l
        A_global=A_global+A{j}*w{j};
        B_global=B_global+B{j}*w{j};
    end
    x=A_global*x;
end
time=zeros(1,10);
for i=1:1:num
    time(i)=i;
end

figure()
plot(time, U, '-r' , 'linewidth',2);
xlabel('Time in samples');ylabel('u');

figure()
plot(time, XX1(1,:),'-r', time,XX2(1,:),'-.b',  'linewidth',2);
xlabel('Time in samples');ylabel('x1');
h=legend('$$System\_withController(t)$$','$$System\_withoutController(t)$$');
set(h,'interpreter','latex');

figure()
plot(time, XX1(2,:),'-r', time,XX2(2,:),'-.b',  'linewidth',2);
xlabel('Time in samples');ylabel('x2');
h=legend('$$System\_withController(t)$$','$$System\_withoutController(t)$$');
set(h,'interpreter','latex');
%出现问题，无解，可能确实存在问题，也可能是系统确实无解
%因此，我需要在状态反馈上尝试

在状态反馈上考虑干扰，仍然无解。

clear all;yalmip('clear');
 
% A{1}=[0.1 0.4;-0.8 0.9];B{1}=[1;1];
% A{2}=[0.85 1;-0.3 0.9];B{2}=[1;1];
% A{3}=[1 0;-0.3 1.3];B{3}=[1;1];
A{1}=[0.1 0.4;-0.8 0.9];B{1}=[1;1];D1{1}=[0.1;0.1];
C2{1}=[1 1];D2{1}=[0.1];
A{2}=[0.85 1;-0.3 0.9];B{2}=[1;1];D1{2}=[0.1;0.1];
C2{2}=[1 1];D2{2}=[0.1];
A{3}=[1 0;-0.3 1.3];B{3}=[1;1];D1{3}=[0.1;0.1];
C2{3}=[1 1];D2{3}=[0.1];
 
[~,l]=size(A);s=2;%被控对象模糊规则数，控制器模糊规则数
[n_x,n_u]=size(B{1});
[n_z,n_d]=size(D2{1});
 
mu_1=[1;0]; 
E=[eye(n_x) zeros(n_x,n_u);zeros(n_u,n_x) zeros(n_u,n_u)];
H=[zeros(n_x,n_u);mu_1(1)*eye(n_u);zeros(n_x,n_u);mu_1(2)*eye(n_u);zeros(n_d,n_u)];
 
P=sdpvar(n_x+n_u,n_x+n_u);
XX1=sdpvar(2*(n_x+n_u)+n_d,n_x,'full'); 
GAMMA2=sdpvar(1,1,'full');
for j=1:1:s
    G{j}=sdpvar(n_u,n_u);
    GK{j}=sdpvar(n_u,n_x,'full');
end

for i=1:1:l
    C2Bar=[C2{i} zeros(n_z,n_u)];
    Sigma_d=mdiag(P,[-E'*P*E+C2Bar'*C2Bar C2Bar'*D2{i};D2{i}'*C2Bar D2{i}'*D2{i}-GAMMA2*eye(n_d)]);
    for j=1:1:s
        Xi{i}{j}=Sigma_d...
        +[-XX1 -H*G{j} [XX1*A{i}+H*GK{j} XX1*B{i}-H*G{j} XX1*D1{i}]]...
        +[-XX1 -H*G{j} [XX1*A{i}+H*GK{j} XX1*B{i}-H*G{j} XX1*D1{i}]]';
    end
end
 
LMI_S=[lmi(P>=0):'P'];  
for i=1:1:l
    for j=1:1:s
        LMI_S=LMI_S+[lmi(Xi{i}{j}<=0):'LMI_S'];
    end
end
 
options = sdpsettings('solver','mosek','shift','0.00');  
solution = optimize(LMI_S, [],options)
check(LMI_S)
 
gamma=sqrt(value(GAMMA2))
for j=1:1:s
    GG{j}=value(G{j})
    GGKK{j}=value(GK{j})
    K{j}=inv(GG{j})*GGKK{j}
end
 
for i=1:1:l
    A{i}
    eig(A{i})
    for j=1:1:s
        AA{i}{j}=A{i}+B{i}*K{j}
        eig(AA{i}{j})
    end
end
 
num=100;x0=[-1;2];
 
XX1=zeros(2,10);
x=x0;
for i=1:1:num
    XX1(:,i)=x;
    %计算隶属函数
    w1_down=(x(1)+2)/10;
    w1_up=(x(1)+3)/10;
    w2_down=(-x(1)+2)/10;
    w2_up=(-x(1)+4)/10;
    w1=(1-(sin(x(1)))^2)*w1_down+((sin(x(1)))^2)*w1_up;
    w2=0.5*w2_down+0.5*w2_up;
    w3=1-w1-w2;
    w{1}=w1;w{2}=w2;w{3}=w3;
    h{1}=w1/(w1+w2);h{2}=w2/(w1+w2);
    K_global=0;
    for j=1:1:s
        K_global=K_global+K{j}*w{j};
%         K_global=K_global+K{j}*h{j};%这应该才是对的，为什么最后不稳定在0处？
    end
    u=K_global*x;
    A_global=0;B_global=0;
    for j=1:1:l
        A_global=A_global+A{j}*w{j};
        B_global=B_global+B{j}*w{j};
    end
    x=A_global*x+B_global*u;
end
XX2=zeros(2,10);
x=x0;
for i=1:1:num
    XX2(:,i)=x;
    %计算隶属函数
    w1_down=(x(1)+2)/10;
    w1_up=(x(1)+3)/10;
    w2_down=(-x(1)+2)/10;
    w2_up=(-x(1)+4)/10;
    w1=(1-(sin(x(1)))^2)*w1_down+((sin(x(1)))^2)*w1_up;
    w2=0.5*w2_down+0.5*w2_up;
    w3=1-w1-w2;
    w{1}=w1;w{2}=w2;w{3}=w3;
    A_global=0;B_global=0;
    for j=1:1:l
        A_global=A_global+A{j}*w{j};
        B_global=B_global+B{j}*w{j};
    end
    x=A_global*x;
end
time=zeros(1,10);
for i=1:1:num
    time(i)=i;
end
 
figure(1)
plot(time, XX1(1,:),'-r', time,XX2(1,:),'-.b',  'linewidth',2);
xlabel('Time in samples');ylabel('x1');
h=legend('$$System\_withController(t)$$','$$System\_withoutController(t)$$');
set(h,'interpreter','latex');
 
figure(2)
plot(time, XX1(2,:),'-r', time,XX2(2,:),'-.b',  'linewidth',2);
xlabel('Time in samples');ylabel('x2');
h=legend('$$System\_withController(t)$$','$$System\_withoutController(t)$$');
set(h,'interpreter','latex');

也许是系统不行，但更可能是方法不对——我试一试。

如果是系统不行，我将γ设置成10000，系统应该是有解的（因为之前的仿真中已经证明出了结果）。因此，我的方法存在问题。

我找到了问题所在，问题在于我忘了在增广系统中考虑 z(k)，修改之后的程序如下

Jack_Swan
关注
两种方法仿真的比较
Jack_Swan的博客
12-11 850
Shur补无疑是简单好用的方法，但是投影定理有例程，更炫酷。 非奇异描述方式当然更好，但是奇异可以减少耦合项。 现在我只掌握了Shur补的一般系统仿真，奇异和投影都用不好。 综上，一共有三种： 1、Shur补一般系统（常见，已可）； 2、Shur补奇异系统（这一种我不确定是否可行，我认为可行，但是没有见过）； 3、投影定理奇异系统（王老师用，有王老师程序）。 ...
蛋蛋的问题（一）
最有趣的编程
06-05 318
蛋蛋的问题（一） 时间限制：1000 ms  |  内存限制：65535 KB 难度：2 描述 2013年的小菜鸟们终于到来了，为了你们能在南工里茁壮的成长，学长们夜以继日的给他们教c语言知识，经过一周的学习，这群小菜鸟们终于知道什么是c语言了，为了考验你们学的怎么样了，学长们给你们出了一堆题，这就是其中一个。。。。   题目很简单，就是考验你们的排序水
蛋蛋成长
哈尼熊熊的博客
10-15 374
1. 第一次秋游（2018-10-12 辰山植物园 10点到，14点回） 不愉快。早上开秋游动员会的时候被前排小朋友来了一拳，开始哭，哭的没法照相。中间非要抱抱，不抱就哭，刚到展览馆就哭睡着了，睡到走，叫醒，继续哭。 和小朋友一起参加学校活动还是很好的，可以观察在校期间，团体活动，交朋友，自理能力等。1.集体活动时死活不参加，非要在我身边。（不过后来有问老师他平时在学校里还是可以的，应该是请假...
蛋蛋的问题（二）
最有趣的编程
06-05 282
蛋蛋的问题（二） 时间限制：100 ms  |  内存限制：65535 KB 难度：2 描述 教新生的时候，也是学长们最兴奋的时候，但是只局限于讲课的那位帅帅的学长，在最后一排坐着听课的学长确实无比的苦闷啊，为了让小学弟知道，你们的二进制转换学的很烂，学长为你们创造了这道很难的题，没看错，是很难。。。 二进制的转换很简单，下面是搜搜里面的介绍：
蛋蛋读NVMe之二: 吉祥三宝
Linux阅码场
07-06 4160
蛋蛋读NVMe之一：为什么刘备需要NVMe上回书说道，NVMe有三宝：SQ,CQ和DB。接下来我们就详细的看看这吉祥三宝。Host往SQ中写入命令， SSD往CQ中写入命...
[NOI2011] 兔兔与蛋蛋游戏题解
Sinnan
03-20 4502
这些天，兔兔和蛋蛋喜欢上了一种新的棋类游戏。 这个游戏是在一个nn行mm列的棋盘上进行的。游戏开始之前，棋盘上有一个格子是空的，其它的格子中都放置了一枚棋子，棋子或者是黑色，或者是白色。 每一局游戏总是兔兔先操作，之后双方轮流操作，具体操作为： 兔兔每次操作时，选择一枚与空格相邻的白色棋子，将它移进空格。 蛋蛋每次操作时，选择一枚与空格相邻的黑色棋子，将它移进空格。 第一个不能按照规则操作的人输掉游戏。为了描述方便，下面将操作“将第x行第y列中的棋子移进空格中”记为M(x,y)M(x,y)。
蛋蛋读NVMe之一：为什么刘备需要NVMe
Linux阅码场
06-22 2815
没有前戏，直接进入。蛋蛋就是这么个人。NVMe是一种Host与SSD之间通讯的协议，它在协议栈中隶属高层。 NVMe在协议栈中处于应用层或者命令层，它是指挥官，军师，在...
P1971 [NOI2011] 兔兔与蛋蛋游戏 题解
yaosichengalpha的博客
06-29 973
P1971 [NOI2011] 兔兔与蛋蛋游戏 题解
[NOI2011]兔兔与蛋蛋游戏
aoanping0730的博客
01-04 153
Description Input 输入的第一行包含两个正整数 n、m。 接下来 n行描述初始棋盘。其中第i 行包含 m个字符，每个字符都是大写英文字母"X"、大写英文字母"O"或点号"."之一，分别表示对应的棋盘格中有黑色棋子、有白色棋子和没有棋子。其中点号"."恰好出现一次。 接下来一行包含一个整数 k（1≤k≤1000） ，表示兔兔和蛋蛋各进行了k次操作...
蛋蛋助手v1.46+增加了DBA设计
06-23
此外，从v1.45版本的“增加了导出SQL语句数据的功能”来看，蛋蛋助手已经具备了强大的数据导出能力。结合v1.46的DBA设计，我们可以推测，现在用户不仅可以方便地导出数据，还能进行更为复杂的SQL操作，如数据分析、...
JAVA小游戏源码---蛋蛋快跑
12-14
从给定的文件信息来看，这是一款名为“蛋蛋快跑”的Java小游戏的源代码片段。这款游戏通过Java语言实现，具备一定的可玩性，主要涉及到了面向对象编程、线程控制、用户界面交互以及基本的游戏逻辑处理等核心IT知识点...
深入java虚拟机（七）深入源码看java类加载器ClassLoader 1
08-04
《深入Java虚拟机（七）深入源码看java类加载器ClassLoader》 Java类加载器（ClassLoader）在Java运行环境中扮演着至关重要的角色。它负责将类的字节码加载到Java虚拟机（JVM）中，使得程序能够运行。ClassLoader是...
Matlab如何配置小波工具（Wavelet Toolbox）
weixin_63655126的博客
09-11 295
返回上一目录，直接复制wavelet文件夹，找到你下载matlab是的文件夹位置，如我的文件夹在D盘program files中。没有这个工具就要去下载了。已经有了小波工具了，运行代码也可以正常运行了，不会再报错了。点击toolbox，把文件夹粘贴在其中。找到文件夹路径后保存即可。圈圈转完点击确定就好了。最后检测是否配置成功。
Matlab,vivado/modelsim,python的数据交互
FPGA/数字IC/通信算法
09-11 756
在进行verilog模块的开发过程中，常常需要对数据处理过程进行行为级仿真，以验证代码的功能逻辑是否正确，因此需要将RTL级仿真的结果与MATLAB或python的处理结果进行对比验证。目前随着python的流行，利用python 生成仿真数据，并与verilog进行数据交互的情形也非常之多。Verilog在读取仿真数据的时候，有时候也需要进行仿真数据文件的输出，常用的函数有$fwrite(file_name)、$fdisplay(file_name),相应的脚本这里一并给出。
MATLAB中的异常处理机制：掌握错误和警告的管理
最新发布
2401_85812026的博客
09-14 613
异常处理是确保程序在遇到错误时能够优雅地恢复或终止的一种方法。在MATLAB中，异常处理主要通过try/catch语句实现。try块包含可能会引发异常的代码，而catch块则包含异常发生时执行的代码。
MATLAB图像处理
shiming8879的博客
09-12 773
同时，随着机器学习和深度学习技术的不断发展，MATLAB在图像处理领域的应用也将更加深入和广泛。MATLAB，作为美国MathWorks公司出品的商业数学软件，以其强大的矩阵运算能力和丰富的函数库，在图像处理领域得到了广泛的应用。MATLAB不仅提供了基础的图像处理功能，还通过图像处理工具箱（Image Processing Toolbox）等高级工具，为用户提供了从图像读取、显示、转换到高级分析和处理的一系列功能。MATLAB提供了多种图像分割技术，如基于阈值的分割、基于区域的分割、基于边缘的分割等。
分类预测|基于蜣螂优化随机森林的数据分类预测Matlab程序DBO-RF 多特征输入多类别输出 含基础RF程序
2401_86241083的博客
09-14 1271
分类预测|基于蜣螂优化随机森林的数据分类预测Matlab程序DBO-RF 多特征输入多类别输出 含基础RF程序
刚接触无处下手？水下航行器AUV/UUV六自由度模型/控制器设计matlab/simulink参考代码，基础的/进阶的，入门到顺利毕业/完成课题/发表论文。
xiang_123_的博客
09-13 393
水下航行器AUV/UUV六自由度数学matlab代码模型/控制器设计
蛋蛋 北京浮生记 java
12-23
《蛋蛋 北京浮生记》是一部由北京作家蛋蛋创作的小说。小说以北京市为背景，描述了主人公在这个繁华都市中的生活、情感和成长。小说描绘了北京的城市风景和人文景观，让读者感受到了这座城市的独特魅力。同时，小说...
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值