文件1 ds_main.m

clear all;close all;clc; %%%%%%%%%% D_S算法設(shè)置 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%num_Sensor=3; %傳感器數(shù)目 num_Object=6; %識別對象數(shù)目 四個目標(機載雷達: 火控 瞄準 導航 預(yù)警)+全集+空集 num_Period=3; %測量周期數(shù)Info=zeros(num_Sensor,num_Object,num_Period); % 一次判決所需要的信息 3*6*3 %% Info 3*6 Info(:,:,1)=[0.30 0.40 0.15 0.00 0.15 0.00;0.30 0.50 0.10 0.00 0.10 0.00;0.30 0.30 0.20 0.00 0.20 0.00;]; Info(:,:,2)=[0.40 0.20 0.20 0.00 0.20 0.00;0.50 0.20 0.20 0.00 0.10 0.00;0.50 0.30 0.10 0.00 0.10 0.00;]; Info(:,:,3)=[0.50 0.20 0.15 0.00 0.15 0.00;0.40 0.30 0.10 0.00 0.20 0.00;.40 0.20 0.10 0.00 0.30 0.00;]; Info1=zeros(num_Period,num_Object); %% 3*6%% 周期內(nèi)各傳感器間的數(shù)據(jù)融合for i=1:num_PeriodInfo1(i,:)=Info(1,:,i); %%Info表示 第i個周期的第1行直接賦值給Info1(1*6）for j=1:num_Sensor-1Info1(i,:)=DS_fusion(Info1(i,:),Info(j+1,:,i));%% DS_fusion(1*6 , 第i個周期的j+1行（行向量）)%% DS_fusion（第i行，第i+1行）end end%% 周期之間的數(shù)據(jù)融合Result=Info1(1,:);for i=1:num_Period-1 Result=DS_fusion(Result,Info1(i+1,:)); end %% ****決策判斷**************************** ec1=0.1; %% 門限值ec ec2=0.1; DS_judg(Result,ec1,ec2,num_Object)

文件2 ds_fusion.m

%% 功能：融合x,y兩行向量(經(jīng)典Dempster-Shafer組合公式) % x,y的格式形如[m1 m2 m3, ... , mk, m(全集), m(空集)] % 要求m1 m2 m3 ...之間互相無交集 % m(全集)可不為0，表示不確定度 % m(空集)肯定是0 %% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %% （前一行，后一行） function x=DS_fusion(x,y)%% DS_fusion(1*6 , 1*6) %% 判斷輸入x、y為行向量；列數(shù)相等 [nx,mx]=size(x); if 1~=nxdisp('x應(yīng)為行向量');return; end [ny,my]=size(y); if 1~=nydisp('y應(yīng)為行向量');return; end if mx~=mydisp('x,y列數(shù)應(yīng)相等');return; end%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% temp=0; for i=1:mx-1 %% nx=1 mx=6if i==mx-1x(1,i)=x(1,i)*y(1,i); %% 對全集的特殊處理elsex(1,i)=x(1,i)*y(1,i)+x(1,i)*y(1,mx-1)+y(1,i)*x(1,mx-1);endtemp=temp+x(1,i); end for i=1:mx-1x(1,i)=x(1,i)/temp; end x(1,mx)=0;

文件3 ds_judg.m

function x=DS_judg(Result,ec1,ec2,num_Object) res=zeros(1,num_Object); a1=max(Result); [c,r]=find(Result==a1); res=Result; res(c,r)=0; a2=max(res); if a1-a2>ec1 %% &&a1>Result(1,num_Object-1)&&Result(1,num_Object-1)<ec2y1=1if Result(1,num_Object-1)<ec2y2=1if a1>Result(1,num_Object-1)y3=1x=[c,r,a1];endend else x=[0,0,0]; end

總結(jié)

以上是生活随笔為你收集整理的d-s证据理论 matlab 完整代码的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。