cuijizhe's DSP Using MATLAB

프랙탈 주식구조 분석코드

2011-05-23T15:27:00.001+09:00

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Fractal stock price simulator
% M. Hanchak, 18OCT10
% Reference:
% Benoit B. Mandelbrot, A Multifractal Walk down Wall Street,
% Scientific American, February 1999
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear
% User inputs %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
low_price = 15;
high_price = 55;
start_time = 0;
end_time = 10;
multiplier = 2/3;
num_dimensions = 7;
% setup %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
y = [low_price high_price];
x = [start_time end_time];
tol = (high_price-low_price) * multiplier;
hold all

% Loop through number of fractal levels %%%%%%%%%%%%%%%%
for i = 1:num_dimensions
N = 3^i+1; % keep subdividing by three
xn = linspace(x(1),x(end),N); % make new x vector
y = interp1(x,y,xn); % interpolate values

% add/subtract at key points to create fractal %%%%%%%%%%%%%%%%%%%%%%%%
for j = 2:N-1
if mod(j-1,3) ~= 0 % don't affect previous points.
% add to existing values - alternate sign and reduce magnitude
y(j) = y(j) + tol/exp(i/2)*(-1)^(j);
end
end

% plotting %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
x = xn;
plot(x,y);drawnow;pause(0.2)
axis([start_time,end_time,low_price,high_price])
end

디지탈 신호처리, 영상처리, 이미지 처리, 이미지 프로세싱, 공간기하학 정수론, 수학기초, 3D, mesh, 그래프 이론, 3차원 이미지, 카오스와 프랙탈, 정보의 표현, 코드, MATLAB, 프로그래밍, 워터마킹, 퀴즈, 이야기, 변환, 푸리에, DCT, MATLAB기초, 영상처리코딩, JPEG, 허프만, 중국어자료, 웨이블렛, 필터링, 전문용어, IT용어, 오디오,큐브,점처리, 미디언 필터, RST변환, 로테이션,스케일링, Rotation, Scaling, reverse

Matlab에서 corrcoef로 얻은 Array

2011-01-22T09:58:00.001+09:00

  correlation=corrcoef(V1, Vmodel1)

 correlation  =1.0000   0.9317   （V1와V1의 상관도는 1，V1와 Vmodel1의 상관도는 0.93317）

   0.9317    1.0000  V1와Vmodel1의 상관도는0.93317,Vmodel1 자신과의 상관도는 1)

콘벌루션: conv([1 1 1],[1 1 1]) 실행결과는 1 2 3 2 1에 관한 해석

2011-01-21T17:15:00.001+09:00

콘벌루션의 의미를 간단하게 이해하면 학교에서 배웠던 다항식의 적이다.

만약 아래의 식을 실행한다면,
conv([a b],[c d])
결과는 다음과 같다.
ac ad+bc bd

따라서 위 식의 해를 구하는 과정은
conv([1 1 1],[1 1 1])
1*1=1
1*1+1*1=2
1*1+1*1+1*1=3
1*1+1*1=2
1*1=1

이렇게 된다.

matlab에서 help conv로 찾아보면 다음과 같다.

CONV Convolution and polynomial multiplication.
C = CONV(A, B) convolves vectors A and B. The resulting vector is
length MAX([LENGTH(A)+LENGTH(B)-1,LENGTH(A),LENGTH(B)]). If A and B are
vectors of polynomial coefficients, convolving them is equivalent to
multiplying the two polynomials.

C = CONV(A, B, SHAPE) returns a subsection of the convolution with size
specified by SHAPE:
'full' - (default) returns the full convolution,
'same' - returns the central part of the convolution
that is the same size as A.
'valid' - returns only those parts of the convolution
that are computed without the zero-padded edges.
LENGTH(C)is MAX(LENGTH(A)-MAX(0,LENGTH(B)-1),0).

콘벌루션: conv([1 1 1],[1 1 1]) 실행결과는 1 2 3 2 1에 관한 해석

2011-01-21T15:51:00.001+09:00

 콘벌루션의 의미를 간단하게 이해하면 학교에서 배웠던 다항식의 적이다.

 만약 아래의 식을 실행한다면,  conv([a b],[c d])  결과는 다음과 같다. ac ad+bc bd

 따라서 위 식의 해를 구하는 과정은 conv([1 1 1],[1 1 1])  1*1=1 1*1+1*1=2 1*1+1*1+1*1=3 1*1+1*1=2 1*1=1

 이렇게 된다.

 matlab에서 help conv로 찾아보면 다음과 같다.

 CONV Convolution and polynomial multiplication.     C = CONV(A, B) convolves vectors A and B.  The resulting vector is     length MAX([LENGTH(A)+LENGTH(B)-1,LENGTH(A),LENGTH(B)]). If A and B are     vectors of polynomial coefficients, convolving them is equivalent to     multiplying the two polynomials.       C = CONV(A, B, SHAPE) returns a subsection of the convolution with size     specified by SHAPE:       'full'  - (default) returns the full convolution,       'same'  - returns the central part of the convolution                 that is the same size as A.       'valid' - returns only those parts of the convolution                  that are computed without the zero-padded edges.                  LENGTH(C)is MAX(LENGTH(A)-MAX(0,LENGTH(B)-1),0).

corrcoef은 어떻게 구한 것인가?

2011-01-21T15:44:00.001+09:00

 corrcoef는 상관계수를 구하는 것임 확율과 통계에서의 ρ(rou) 공식은 그림 참고

Matlab에서 corrcoef로 얻은 Array

2011-01-21T15:41:00.001+09:00

 correlation=corrcoef(V1,Vmodel1) correlation =     1.0000    0.9317   （V1와V1의 상관도는 1，V1와 Vmodel1의 상관도는 0.93317）     0.9317    1.0000   （V1와Vmodel1의 상관도는0.93317,Vmodel1 자신과의 상관도는 1）

PN Code, PN Sequence 만드는 함수

2010-11-22T18:51:00.001+09:00

function [P,l] = pnseq(pu,varargin)

% PNSEQ PN codes.

% [P,l] = PNSEQ(PU) generates a pseudonoise (PN) sequence P from a

% binary linear feedback shift register. l gives the length of the

% code.

% The shift register is described by the generator polynomial PU.

% PNSEQ(...,RU) generates a PN sequences given RU, the initial state of

% the shift register. RU must be a vector of the same length as the

% generator polynomial.

% Note that a {0,1}-valued binary sequence corresponds to a

% cos(pi*{0,1}) = {1,-1}-transmitted sequence.

% Joe Henning - Fall 2006

P = 0;

l = 0;

error(nargchk(1,2,nargin));

[ru,n,msg] = parseinput(pu,varargin{:});

error(msg);

l = 2^n - 1;

P = zeros(1,l);

for k = 1:l

temp = mod( fliplr(pu)*ru.', 2 );

P(k) = ru(n);

ru = [temp ru(1:n-1)];

end

%--------------------------------------------------------------------------

function [ru,n,msg] = parseinput(pu,varargin)

% Parse the input and determine optional parameters:

% Outputs:

% ru - initial values of register u

% n - number of register bits

% msg - possible error message

% Set some defaults:

ru = 0;

n = 0;

msg = '';

[mu,nu] = size(pu);

if ( (mu ~= 1 & nu ~= 1) | (mu == 1 & nu == 1) )

msg = 'pu must be a vector.';

return

end

n = length(pu);

% Assume no register values until proven otherwise

found_ru = 0;

ru = ones(1,n);

switch nargin

case 2

% is ru

if isnumeric(varargin{1})

if length(varargin{1}) == n

found_ru = 1;

ru = varargin{1};

end

if ~found_ru

fprintf('Using default ru.\n');

end

논문데이터 정리: 삽입알고리즘 절차

2010-11-21T20:10:00.001+09:00

논문에 사용될 이미지

데이터분석: 발표할 논문데이터 정리(2010년11월)

2010-11-21T20:05:00.001+09:00

아직까지 정리할 것이 너무 많다.

개발이라고 한 알고리즘의 성능이 기대이하이다...

좀 더 분발해야겠는데...

최종욱: 워터마크 기술을 미래의 저작권 보호기술로 인식

2010-11-21T19:58:00.001+09:00

지난해 8월 영국에서 열렸던 "디지털 워터마킹 국제워크숍(IWDW(International Workshop on Digital Watermarking)"에 다녀왔다. 이 국제 워크숍은 워터마크 기술을 연구하는 학자와 연구자들이 한자리에 모여 새로 개발한 기술을 발표하는 자리다. 이 워크숍은 한국에서 시작됐고 지금도 조직과 진행을 한국이 주도하고 있다. 문제는 한국에서 시작된 이 워크숍이 중국의 연구 논문으로 도배하다시피 진행되고 있으며 한국은 완전히 소수자그룹으로 전락하고 있다는 점이다.
26편의 논문 중에서 중국교수들이나 연구원들이 발표한 논문은 16편으로 다수를 점하고 있고 한국대학에서 발표한 논문은 단 1편 뿐이었다. 발표자 명단 74명 중에서 한국 사람은 단지 두 명 뿐이라는 점도 충격적이다.

워터마크 기술은 미래의 저작권 보호기술로 인식되고 있다. 중국정부가 이를 ‘미래기술’로 인식하고 연구개발을 늘리고 있고, 엄청난 연구인력들이 달려들고 있다. 더구나 이 기술은 중요한 메시지를 음악이나 영화 속에 숨겨서 보낼 수 있기 때문에 군사적인 용도로도 아주 중요하고, 북한에서 나온 수퍼 노트와 같은 정교한 위조 지폐를 근본적으로 막을 수 있는 유일한 기술이기도 하다. 한국에서 주도적으로 시작된 이 분야에 왜 한국교수들의 참여는 줄어들다 못해 전멸해가고 있는가. 한국정부의 R&D조급증과 연구과제 심사에 엄격하게 적용되는 ‘중복과제 방지’ 명분 때문이다. 최근 연구제안서를 제출하면 ‘워터마크 기술은 이미 다 개발됐다’는 논리와 ‘국민의 혈세를 가지고 같은 이름의 연구에 중복 투자할 수 없다’는 정부의 방침 때문에 대부분의 워터마크 연구 제안서가 심사에서 기각됐다.

연구원들이 중복과제라는 것을 피하기 위해 워터마크라는 이름대신 ‘데이터 은닉기술’ ‘포렌직 마크기술’ 등으로 위장하기도 했다. 워터마크 분야는 아직도 해결되지 않은 기술적 난제도 많고, 새로운 세부분야도 지속적으로 등장하고 있다. 지금도 워크숍과 학술대회를 통해 수많은 논문이 발표되고 토론이 벌어지고 있는 이유이기도 하다. 현재 전자 업계에서 개발하고 있는 3DTV는 새로운 워터마크 연구분야다. 워터마크에 은닉되어 송수신되고 있는 메시지 은닉 기술이나 화폐삽입기술도 계속적인 연구가 필요하다. 단지 워터마크라는 이름 하나로 묶어 중복과제 방지의 이유로 지속적인 연구를 중단시킨다면 R&D를 통한 우리 나라의 미래는 없다고 본다.

문제는 훌륭한 프로젝트를 발굴해내고, 어떻게 좋은 성과를 내느냐에 달려 있다. 한국처럼 프로젝트를 수행하기 위해 중복과제인지 먼저 심의하고, 수행 후에는 예산 집행이 적정했는지 여부로 프로젝트 성공을 따지는 풍조로는 좋은 프로젝트 수행이 어렵다는 것이다. 이제 한국은 더 이상 다른 집단을 잠재적 범죄자 취급을 하는 불신사회가 아니다.

연구 제목을 가지고 중복과제인지 아닌지를 판단할 수 있는 단세포적인 사회는 더욱 아니다. 연구를 집행하는 부처의 관리들이 모든 권한을 가지고 과제의 미래 가치와 성공여부를 판단할 수 있도록 해야 한다. 지금같은 방식을 고집하는 한 미래의 대한민국 R&D는 쓸 데없는 중복과제 판단에 매달려 정말로 필요한 연구개발을 놓치고, 엉뚱한 혈세의 낭비로 흐를 수 있다. 미국 유니버셜뮤직에 납품할 정도로 국내의 기술이 인정받고 있지만 이런 추세가 지속되면 저작권 보호를 위해 중국 기술과 제품을 사다 써야 할는지 모른다는 우려를 하게 된다.

[미래포럼] 대한민국의 R&D조급증

대만 Matlab신호처리 매니아 사이트

2010-11-17T12:58:00.001+09:00

http://mirlab.org/jang/

툴박스 다운로드: 오디오신호처리관련 Utility Toolbox(1)

2010-11-17T12:51:00.001+09:00

툴박스 다운로드: Speech and Audio Processing Toolbox(1)

2010-11-17T12:47:00.001+09:00

http://mirlab.org/jang/matlab/toolbox/sap.rar

이미지 워터마킹(삽입: 2-D wavelet transform)

2010-11-17T08:59:00.001+09:00

%Implementing Watermarking using discrete 2-D wavelet transform.

%Input image is watermarked with a key having Mean = 0 & Variance = 1

clc;clear all;close all;

img = imread('home.jpg'); %Get the input image

img = rgb2gray(img); %Convert to grayscale image

img = double(img);

c = 0.01; %Initialise the weight of Watermarking

figure,imshow(uint8(img)),title('Original Image');

[p q] = size(img);

%Generate the key

n = awgn(img,4,3,'linear');

N = imabsdiff(n,img);

figure,imshow(double(N)),title('Key');

[Lo_D,Hi_D,Lo_R,Hi_R] = wfilters('haar');%Obtain the fiters associated with haar

[ca,ch,cv,cd] = dwt2(img,Lo_D,Hi_D); %Compute 2D wavelet transform

%Perform the watermarking

y = [ca ch;cv cd];

Y = y + c*abs(y).* N;

p=p/2;q=q/2;

for i=1:p

for j=1:q

nca(i,j) = Y(i,j);

ncv(i,j) = Y(i+p,j);

nch(i,j) = Y(i,j+q);

ncd(i,j) = Y (i+p,j+q);

end

%Display the Watermarked image

wimg = idwt2(nca,nch,ncv,ncd,Lo_R,Hi_R);

figure,imshow(uint8(wimg)),title('Watermarked Image');

diff = imabsdiff(wimg,img);

figure,imshow(double(diff));title('Differences');

wavread를 이용하여 오디오를 원하는 길이로 자르기

2010-11-14T16:13:00.001+09:00

[original_audio, FS, NBITS, OPTS] = wavread('music\Yesterday Once More.wav',[1, 4000000]);

wavwrite(original_audio, FS, NBITS, 'Yesterday Once More.wav');

% 위에서 40000000으로 설정했는데(약1분30초), 이부분 수치를 원하는대로 바꿔서 짜르면 됨.

Matlab으로 투명한 실행창 만들기(가능할지?)

2010-11-14T12:16:00.001+09:00

안녕하세요?

Matlab으로 프로그래밍하다가 원론적인 문제에 부닥쳐서 여기에 글을 올립니다.

아래와 같은 일을 하려는데요.. 여러분의 의견을 부탁드립니다.

문제: 실행창을 투명하게 설정하는 방법

- GUI를 사용해도 되고, 어떤 방식을 사용하던 상관없이 투명한 실행창을 만들 수 있어야 됩니다.

설명을 드리면, 프로그램을 실행했을 때, 백그라운드 화면(예를 들면, Windows XP화면)도 실행창을 통해서 볼 수 있어야 됩니다. 실행은 유리(투명성 조절가능)와 같은 면에서 진행하게 되는거죠?

가능할가요?

의견부탁드립니다.

본 블로그의 임의의 곳을 선택해서 댓글 달아주시면 됩니다. 감사합니다.

E-mail: cuijizhe@gmail.com

신호 대 잡읍 비 구하기 코드: SNR(오디오용)

2010-11-12T15:03:00.002+09:00

%신호대 잡음 비 구하기

신호파워 (오지지날)^2
SNR=----------=-------------------------------
잡음파워 (잡음포함신호- 오리지날)^2

function snr=SNR_singlech(I,In)
% 2010.11.12 by cuijizhe
% I :original signal
% In:noisy signal(ie. original signal + noise signal)

snr=0;
Ps=sum(I^2); %signal power
Pn=sum((I-In).^2); %noise power
snr=10*log10(Ps/Pn);

% 끝

%사용사례

%예1:

clear all; clc; close all;
[filename,pathname]=uigetfile('*.wav','파일선택하세유：');
[X,fs]=wavread([pathname filename]);
[Y,NOISE] = noisegen(X,10);
subplot 311; plot(X);
subplot 312; plot(NOISE);
subplot 313; plot(Y);
mn=mean(NOISE)
snr=SNR_singlech(X,Y)

%예2

clear all; clc; close all;
[filename,pathname]=uigetfile('*.wav','오디오파일：');
[filename1,pathname1]=uigetfile('*.wav','잡음포함오디오파일：');
filepath_name=[pathname1 filename1];
[X,fs]=wavread([pathname filename]);
[Y,NOISE] = add_noisem(X,filepath_name,10,fs);
subplot 311; plot(X);
subplot 312; plot(NOISE);
subplot 313; plot(Y);
mn=mean(NOISE)
snr=SNR_singlech(X,Y)

Noise_CancellationUsing_LMS_algorithm

2010-11-06T11:37:00.001+09:00

Noise_CancellationUsing_LMS_algorithm

FIR filter coefficient estimation

2010-11-06T11:36:00.003+09:00

FIR filter coefficient estimation

AR_process_estimation(code)

2010-11-06T11:36:00.001+09:00

AR_process_estimation 메틀렙 코드

신호 대 잡읍 비 구하기 코드: PSNR(이미지용)

2010-11-04T17:02:00.001+09:00

%========================================%
function snr=SNR(I,In)
% 신호잡음비 구하기 함수

% I :original signal
% In:noisy signal(ie. Original signal + noise signal)
% snr=10*log10(sigma2(I2)/sigma2(I2-I1))

[row,col,nchannel]=size(I);

snr=0;
if nchannel==1 %gray image
Ps=sum(sum((I-mean(mean(I))).^2));%signal power
Pn=sum(sum((I-In).^2));%noise power
snr=10*log10(Ps/Pn);
elseif nchannel==3 %color image
for i=1:3
Ps=sum(sum((I(:,:,i)-mean(mean(I(:,:,i)))).^2));%signal power
Pn=sum(sum((I(:,:,i)-In(:,:,i)).^2));%noise power
snr=snr+10*log10(Ps/Pn);
end
snr=snr/3;
end

Matlab에서 마우스로 좌표값을 얻는 방법

2010-11-04T15:37:00.001+09:00

第一种途径：ginput()函数

ginput提供了一个十字光标使我们能更精确的选择我们所需要的位置，并返回坐标值。函数调用形式为：

[x,y] = ginput(n) % n개 필요할 때..
[x,y] = ginput
[x,y,button] = ginput(...)

对于[x,y] = ginput(n)，能使你从当前的坐标系中读取n个点，并返回这n个点的x，y坐标，均为nX1的向量。可以按回车提前结束读数。

[x,y] = ginput 可以无限的读取坐标直到按下回车键。

[x,y,button] = ginput(...) 返回x和y的坐标，以及button值（1=左键，2=中，3=右）或者按键的ASXII码值。

第二种途径：利用figure的WindowButtonDownFcn属性。当你在图上按下鼠标的时候，可通过该属性定义一个回调程序。回调程序可以是一个有效的Matlab表达式或者一个M文件。

那么为显示当前鼠标按下时的坐标值，我们可以将其定义为一个坐标获取和显示程序。

例如：

% 主函数

function test_mouse_track()
figure;
axis([-10,10,0,5]);
set(gcf,'WindowButtonDownFcn',@ButttonDownFcn);

% 回调函数

function ButttonDownFcn(src,event)
pt = get(gca,'CurrentPoint');
x = pt(1,1);
y = pt(1,2);
fprintf('x=%f,y=%f\n',x,y);

간단한 애니메이션 동영상(출처: http://pinkwink.kr/192)

2010-11-04T15:25:00.001+09:00

AniTest = figure('name','Animation Test','numbertitle','off');

x = [-3 -3.5 -4 -4 -3];

y = [1 1.5 1 -1 -1];

box = fill(x,y, 'r','EraseMode','normal');

grid on

hold on

axis([-5 5 -5 5]);

t = 0:0.01:2*pi;

plot(3*cos(t), 3*sin(t))

plot(17^0.5*cos(t), 17^0.5*sin(t))

plot((3.5^2+1.5^2)^0.5*cos(t), (3.5^2+1.5^2)^0.5*sin(t), 'c')

for t = 0:0.01:2*pi

updatedX = [10^0.5*cos(-(t+pi+atan(1/3))) (3.5^2+1.5^2)^0.5*cos(-(t+pi+atan(3/7))) 17^0.5*cos(-(t+pi+atan(1/4))) 17^0.5*cos(-(t+pi-atan(1/4))) 10^0.5*cos(-(t+pi-atan(1/3)))];

updatedY = [10^0.5*sin(-(t+pi+atan(1/3))) (3.5^2+1.5^2)^0.5*sin(-(t+pi+atan(3/7))) 17^0.5*sin(-(t+pi+atan(1/4))) 17^0.5*sin(-(t+pi-atan(1/4))) 10^0.5*sin(-(t+pi-atan(1/3)))];

set(box, 'Xdata', updatedX,'Ydata', updatedY);

drawnow;

for i=1:5000000

temp = i;

end

영화 만드는 방법(사례)

2010-11-04T15:02:00.001+09:00

% 给你个例子:

%影片生成例子: 旋转一个三维表面绘图

[X, Y, Z]=peaks(50); %创建山峰图形数据

surfl(X, Y, Z) %绘制带光照的表面

axis([-10 10 -10 10 -10 10]) %限制绘图范围

axis vis3d off %三维坐标修正,关闭坐标轴箭头

axis equal %等比例显示三维坐标

shading interp %加个影子美化

colormap(copper) %着色为铜色

for i=1:360 %旋转山峰,从1到360度,每一度捕捉一帧

view(-37.5+i,30) %从水平-37.5+i度,垂直30度的方向看山峰

m(i)=getframe; %每从一个角度看到一张图像(帧),就存储到m

end %每次看完,水平转动1度,直到1周360度;

cla %为播放影片清除坐标

movie(m) %把刚才存到m里的图像连续播放一遍,就是影片了.

Matlab 함수: movie 사용사례

2010-11-04T14:56:00.001+09:00

% 푸리에 변환 과정 보여주기

axis equal
m=moviein(16,gcf);
set(gca,'nextplot','replacechildren');
b=uicontrol('style','slider','position',[100 10 380 15],'min',1,'max',16);

for h=1:16
plot(fft(eye(h+16)));
set(b,'value',h);
m(:,h)=getframe(gcf);
end

clf
axes('position',[0 0 1 1]);
movie(m);

Arnold反置乱(解密)程序之二

2010-11-01T11:14:00.001+09:00

function I1=ArnoldITrans(I,N)

%I1=AmoldITrans(I,N)

%I is original image, I1 is result image, N is the number of transform

% by dynamic of Matlab技术论坛

% see also http://www.matlabsky.com

% contact me matlabsky@gmail.com

% 2009-12-11 14:41:30

I=double(I);

I1=I;

% determine size of cover object

Mc=size(I,1); %Height

Nc=size(I,2); %Width

for k=1:N

for x=1:Mc

for y=1:Nc

x1=2*x-y;

y1=y-x;

if x1>Mc

x1=mod(x1,Mc);

end

if y1>Nc

y1=mod(y1,Nc);

end

if x1<=0

x1=x1+Mc;

end

if y1<=0

y1=y1+Nc;

end

I1(x1,y1)=I(x,y);

end;

I=I1;

end;

Arnold置乱(加密)程序(之一)

2010-11-01T11:12:00.001+09:00

function I1=ArnoldTrans(I,N)

%I1=AmoldTrans(I,N)

%I is original image, I1 is result image, N is the number of transform

% by dynamic of Matlab技术论坛

% see also http://www.matlabsky.com

% contact me matlabsky@gmail.com

% 2009-12-11 14:41:12

I=double(I);

I1=I;

% determine size of cover object

Mc=size(I,1); %Height

Nc=size(I,2); %Width

for k=1:N

for x=1:Mc

for y=1:Nc

x1=x+y;

y1=x+2*y;

if x1>Mc

x1=mod(x1,Mc);

end

if y1>Nc

y1=mod(y1,Nc);

end

if x1==0

x1=Mc;

end

if y1==0

y1=Nc;

end

I1(x1,y1)=I(x,y);

end;

I=I1;

end;

图像置乱度计算MATLAB代码（转存）

2010-11-01T11:10:00.001+09:00

ImageScramblingDegree.m是入口主程序，运行时需要调用ImageEntropy.m程序对图像商进行计算

function S=ImageScramblingDegree(I)

%置乱度计算

%参考文献：张华熊, 吕辉, 翁向军. 基于信息熵的图像置乱程度评价方法[J]. 电路与系统学报, 2007, 12(6): 95-98.

H=ImageEntropy(I);

[M,N]=size(I);

k=1;

uHD=0; uTotal=0;

for i=1:16:M

if (i+15)>M

i=M-15;

end

for j=1:16:N

if (j+15)>N

j=N-15;

end

I16=I(i:i+15,j:j+15);

HD(k)=ImageEntropy(I16);

u=NeiborDiff(I16);

uHD=uHD+u*HD(k);

uTotal=uTotal+u;

k=k+1;

end

S=uHD/(uTotal*H);

%计算分块中互为相邻像素差的平均值,考虑了行相邻和列相邻

function u=NeiborDiff(D)

[M,N]=size(D);

D=double(D);

u=0;

for i=1:M

for j=1:N-1

u=u+abs(D(i,j)-D(i,j+1));

end

for j=1:N

for i=1:M-1

u=u+abs(D(i,j)-D(i+1,j));

end

u=u/(N*(M-1)+M*(N-1));

Matlab图像混沌置乱(转存)-采用的是混沌置乱算法

2010-11-01T11:07:00.001+09:00

首先说明，这个混沌序列置乱其实和很多置乱方式比如随机数置乱的原理相同，但是缺点就是在图像置乱之前首先得生成一个你图像长乘宽大小的混沌序列，如果图像的分辨率比较大，那么在matlab下产生混沌序列的时间比较长，但是图像置乱的效果很好

%图像置乱程序

%作者：ltx1215
%日期：2010年8月7日
%采用的是混沌置乱算法
clear;
[filename, pathname] = uigetfile('*.jpg', '打开原始图像')
filename= [pathname filename];
J=imread(filename);
info=imfinfo(filename);
[m,n,p]=size(J);

%产生混沌序列
x(1)=0.5;
for i=1:m*n-1
x(i+1)=3.7*x(i)*(1-x(i));
end
[y,num]=sort(x);%将产生的混沌序列进行排序
%如果原图为灰度图
if info.BitDepth==8

Scambled=uint8(zeros(m,n));%产生一个与原图大小相同的0矩阵
for i=1:m*n
Scambled(i)=J(num(i));
end

%如果原图为二值图像
elseif info.BitDepth==1
    S=uint8(zeros(m,n));
        for i=1:m
            for j=1:n
        if J(i,j)==1
           S(i,j)=255;
        end
            end
        end
        Scambled=uint8(zeros(m,n));%产生一个与原图大小相同的0矩阵
for i=1:m*n
    Scambled(i)=S(num(i));
end
    end

    %如果为真彩图
    if p==3
    Red=uint8(zeros(m,n));
    Green=uint8(zeros(m,n));
    Blue=uint8(zeros(m,n));
    RedNew=J(:,:,1);
    GreenNew=J(:,:,2);
    BlueNew=J(:,:,3);

    Scambled=uint8(zeros(m,n,p));%产生一个与原图大小相同的0矩阵
for i=1:m*n
    Red(i)=RedNew(num(i));
    Green(i)=GreenNew(num(i));
    Blue(i)=BlueNew(num(i));
end
Scambled(:,:,1)=Red;
Scambled(:,:,2)=Green;
Scambled(:,:,3)=Blue;
    end
imwrite(Scambled,'Scambled.jpg','quality',100);
imshow(Scambled);

………………………分割线………………………

图像反置乱程序：

%图像反置乱程序

%作者：ltx1215

%2010年8月10日

[filename, pathname] = uigetfile('*.jpg', '打开原始图像')

filename= [pathname filename];

J=imread(filename);

info=imfinfo(filename);

[m,n,p]=size(J);

%产生混沌序列

x(1)=0.5;

for i=1:m*n-1

x(i+1)=3.7*x(i)*(1-x(i));

end

[y,num]=sort(x);%将产生的混沌序列进行排序

%如果原图为灰度图

if info.BitDepth==8

IScamble=uint8(zeros(m,n));%产生一个与原图大小相同的0矩阵

for i=1:m*n

IScamble(num(i))=J(i);

end

%如果原图为二值图像

elseif info.BitDepth==1

S=uint8(zeros(m,n));

for i=1:m

for j=1:n

if J(i,j)==1

S(i,j)=255;

end

IScamble=uint8(zeros(m,n));%产生一个与原图大小相同的0矩阵

for i=1:m*n

IScamble(num(i))=S(i);

end

%如果为真彩图

if p==3

Red=uint8(zeros(m,n));

Green=uint8(zeros(m,n));

Blue=uint8(zeros(m,n));

RedNew=J(:,:,1);

GreenNew=J(:,:,2);

BlueNew=J(:,:,3);

IScamble=uint8(zeros(m,n,p));%产生一个与原图大小相同的0矩阵

for i=1:m*n

Red(num(i))=RedNew(i);

Green(num(i))=GreenNew(i);

Blue(num(i))=BlueNew(i);

end

IScamble(:,:,1)=Red;

IScamble(:,:,2)=Green;

IScamble(:,:,3)=Blue;

end

imwrite(IScamble,'IScambled.jpg','quality',100);

imshow(IScamble);

Matlab toolbox: 기하도형 그리기

2010-10-25T15:45:00.001+09:00

Test Image(원과 직선)

2010-10-22T10:49:00.001+09:00

--
Cui Jizhe(崔基哲)
Mobile:+86-189-0433-2188
MSN: korea_china@hotmail.com

ugh 변환_허프만방법으로 원, 직선을 검측하기(07.the_max1 )

2010-10-22T10:46:00.001+09:00

function [rhao,nn]=the_max1(XX,YY,ZZ)
[res]=gather_point(XX,YY,ZZ);nn=0;res
[m,n3]=size(YY);hao=ones(1,3);
for no1=1:1:n3
n=0;nst=n
if no1=(n3)/5
rhao=hao;break;
elseif nn

Hough 변환_허프만방법으로 원, 직선을 검측하기(06.pll.m)

2010-10-22T10:43:00.001+09:00

function [X,Y]=pll(X1,Y1,X2,Y2)

% 直线相交求交点

A1=Y1(1)-Y1(2);
B1=X1(2)-X1(1);
C1=Y1(2)*X1(1)-Y1(1)*X1(2);
A2=Y2(1)-Y2(2);
B2=X2(2)-X2(1);
C2=Y2(2)*X2(1)-Y2(1)*X2(2);
D=det([A1,B1;A2,B2]);
X=det([-C1 B1;-C2 B2])/D;
Y=det([A1 -C1;A2,-C2])/D;

Hough 변환_허프만방법으로 원, 직선을 검측하기(05.machZ_try.m)

2010-10-22T10:40:00.001+09:00

clear all; close all; clc;

i=imread('E:\MATLAB6p5\work\直线在其中.bmp');

[row,range]=size(i);

p=row*range;a=0;Y=zeros(1,p);Q=zeros(1,p);

for k1=1:1:row

for k2=1:1:range

if i(k1,k2)==0

a=a+1;

XXX=k2-1;YYY=row+1-k1;

Y(a)=YYY;Q(a)=YYY-XXX;

end

% for t=1:1:a

% o=[Y(t) Q(t)]

% end

X1=[0 1];Y1=[Y(1) Q(1)];

plot(X1,Y1,'r');hold on;

YY=zeros(1,a-1);QQ=zeros(1,a-1);

for t=2:1:a

X2=[0 1];Y2=[Y(t) Q(t)];

plot(X2,Y2);

[YY(t-1),QQ(t-1)]=pll(X1,Y1,X2,Y2);

% yyy=YY(t-1)

% qqq=QQ(t-1)%看交点

plot(YY(t-1),QQ(t-1),'ko');

end

% [res]=gather_point(YY,QQ,QQ);nn=0;res

% [m,n3]=size(YY);

% for no1=1:1:n3

% n=0;

% for no2=(no1+1):1:n3

% if (abs(res(no2,1)-res(no1,1))<=0.0001&abs(res(no2,2)-res(no1,2))<=0.0001&abs(res(no2,2)-res(no1,2))<=0.0001)

% n=n+1;hao=res(no1,:);

% end

% if nn>=n3/5

% rhao=hao;break;

% elseif nn

Hough 변환_허프만방법으로 원, 직선을 검측하기(04.gather_point.m)

2010-10-22T10:39:00.001+09:00

function [res]=gather_point(xx,yy,zz)%分别为交点x,y,z的集合
[m,n]=size(xx);
res=ones(n,3);
for i=1:1:n
for i1=1:1:3
if i1==1
res(i,i1)=xx(i);
elseif i1==2
res(i,i1)=yy(i);
else
res(i,i1)=zz(i);
end
end
end

Hough 변환_허프만방법으로 원, 직선을 검측하기(03.circle_try_hu.m)

2010-10-22T10:37:00.001+09:00

clear all; close all; clc;

% [cr,st]=circle_product;%st(1)->x,st(2)->y,st(3)->r

cr=imread('E:\MATLAB6p5\work\circle1.bmp');%circle1.bmp

figure;

imshow(cr);

[row,range]=size(cr);

p=row*range;a=0;Y=zeros(1,p);Q=zeros(1,p);

for k1=1:1:row

for k2=1:1:range

if cr(k1,k2)==0

a=a+1;

XXX=k2-1;YYY=row+1-k1;

Y(a)=YYY;Q(a)=XXX;

end

for k3=1:1:p

if Y(k3)==0&Q(k3)==0

break;

end

% z=ones(1,5);

% a=1:1:300;b=1:1:300;

% z1=sqrt((a-Q(1)).^2+(b-Y(1)).^2);

% z2=sqrt((a-Q(2)).^2+(b-Y(2)).^2);

% z3=sqrt((a-Q(3)).^2+(b-Y(3)).^2);

% z4=sqrt((a-Q(4)).^2+(b-Y(4)).^2);

% z5=sqrt((a-Q(5)).^2+(b-Y(5)).^2);

% % z1=sqrt((a-Q(1))*(a-Q(1))+(b-Y(1))*(b-Y(1)));

% % Z1=(a-Q(1)).^2+(b-Y(1)).^2;

% r0=(abs(z1-z2)<=0.1&abs(z1<=z3)<=0.1&abs(z1<=z4)<=0.1&abs(z1<=z5)<=0.1&abs(z2<=z3)<=0.1&abs(z2<=z4)<=0.1&abs(z2<=z5)<=0.1);

% r=r0.*z1;

% aa=r0.*a;

% bb=r0.*b;

% rrr=round(r(r0~=0))

% aaa=aa(r0~=0)

% bbb=bb(r0~=0)

k7=floor(k3/6);k5=2;

a=1:1:300;b=1:1:300;rrr=zeros(1,p);aaa=zeros(1,p);bb=zeros(1,p);k6=0;as=0;k11=1;

for k5=1:1:k7

% k6=0;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

z1=sqrt((a-Q(k5)).^2+(b-Y(k5)).^2);

z2=sqrt((a-Q(k5+1)).^2+(b-Y(k5+1)).^2);

z3=sqrt((a-Q(k5+2)).^2+(b-Y(k5+2)).^2);

z4=sqrt((a-Q(k5+3)).^2+(b-Y(k5+3)).^2);

z5=sqrt((a-Q(k5+4)).^2+(b-Y(k5+4)).^2);

z6=sqrt((a-Q(k5+5)).^2+(b-Y(k5+5)).^2);

r0=(abs(z1-z2)<=.1&abs(z1<=z3)<=.1&abs(z1<=z4)<=.1&abs(z1<=z5)<=.1&abs(z1<=z6)<=.1&abs(z2<=z3)<=1&abs(z2<=z4)<=.1&abs(z2<=z5)<=.1&abs(z2<=z6)<=.1);

rr=r0.*z1;

aa=r0.*a;

bb=r0.*b;

[m,n]=size(rr);

if rr==zeros(1,n);

as=as+1; continue;

end

[pz1,pz]=size(rr(r0~=0));

rrr=round(rr(r0~=0))

aaa=aa(r0~=0)

bbb=bb(r0~=0)

pz2=pz+k11-1;k9=1;

for k8=k11:pz2

RR(k8)=rrr(k9);

XX(k8)=aaa(k9);

YY(k8)=bbb(k9);

k9=k9+1;

end

k11=k8+1;

% if rr==zeros(1,n);

% as=as+1; continue;

% end

end

% if aa==zeros(1,300)

% as=1

% end

rangex=1:1:range;rowy=1:1:row;[RX,RY]=meshgrid(rangex,rowy);

RT=RX;figure;

mesh(RX,RY,RT);hold on;

plot3(XX,YY,RR,'ko');%%%%%%%%%%%%%%%5

[rhao,nn]=the_max1(XX,YY,RR);

i1=ones(row,range);

rhao

center=rhao;

% t2=1.5;t3=1.5;

for k1=1:1:row

for k2=1:1:range

XXX=k2-1;YYY=row+1-k1;

if abs(((XXX-center(1))*(XXX-center(1))+(YYY-center(2))*(YYY-center(2)))-center(3)*center(3))<=101

i1(k1,k2)=0;

end

imwrite(i1,'E:\MATLAB6p5\work\找出的圆.bmp','bmp');

i2=imread('E:\MATLAB6p5\work\找出的圆.bmp');

figure;

imshow(i2);%st

Hough 변환_허프만방법으로 원, 직선을 검측하기(02.circle_try.m)

2010-10-22T10:34:00.001+09:00

clear all; close all; clc;

[cr,st]=circle_product;%st(1)->x,st(2)->y,st(3)->r

[row,range]=size(cr);

p=row*range;a=0;Y=zeros(1,p);Q=zeros(1,p);

for k1=1:1:row

for k2=1:1:range

if cr(k1,k2)==0

a=a+1;

XXX=k2-1;YYY=row+1-k1;

Y(a)=YYY;Q(a)=XXX;

end

for k3=1:1:p

if Y(k3)==0&Q(k3)==0

break;

end

% z=ones(1,5);

% a=1:1:300;b=1:1:300;

% z1=sqrt((a-Q(1)).^2+(b-Y(1)).^2);

% z2=sqrt((a-Q(2)).^2+(b-Y(2)).^2);

% z3=sqrt((a-Q(3)).^2+(b-Y(3)).^2);

% z4=sqrt((a-Q(4)).^2+(b-Y(4)).^2);

% z5=sqrt((a-Q(5)).^2+(b-Y(5)).^2);

% % z1=sqrt((a-Q(1))*(a-Q(1))+(b-Y(1))*(b-Y(1)));

% % Z1=(a-Q(1)).^2+(b-Y(1)).^2;

% r0=(abs(z1-z2)<=0.1&abs(z1<=z3)<=0.1&abs(z1<=z4)<=0.1&abs(z1<=z5)<=0.1&abs(z2<=z3)<=0.1&abs(z2<=z4)<=0.1&abs(z2<=z5)<=0.1);

% r=r0.*z1;

% aa=r0.*a;

% bb=r0.*b;

% rrr=round(r(r0~=0))

% aaa=aa(r0~=0)

% bbb=bb(r0~=0)

k7=floor(k3/6);k5=2;

a=1:1:300;b=1:1:300;rrr=zeros(1,p);aaa=zeros(1,p);bb=zeros(1,p);k6=0;as=0;k11=1;

for k5=1:1:k7

% k6=0;%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

z1=sqrt((a-Q(k5)).^2+(b-Y(k5)).^2);

z2=sqrt((a-Q(k5+1)).^2+(b-Y(k5+1)).^2);

z3=sqrt((a-Q(k5+2)).^2+(b-Y(k5+2)).^2);

z4=sqrt((a-Q(k5+3)).^2+(b-Y(k5+3)).^2);

z5=sqrt((a-Q(k5+4)).^2+(b-Y(k5+4)).^2);

z6=sqrt((a-Q(k5+5)).^2+(b-Y(k5+5)).^2);

r0=(abs(z1-z2)<=.1&abs(z1<=z3)<=.1&abs(z1<=z4)<=.1&abs(z1<=z5)<=.1&abs(z1<=z6)<=.1&abs(z2<=z3)<=1&abs(z2<=z4)<=.1&abs(z2<=z5)<=.1&abs(z2<=z6)<=.1);

rr=r0.*z1;

aa=r0.*a;

bb=r0.*b;

[m,n]=size(rr);

if rr==zeros(1,n);

as=as+1; continue;

end

[pz1,pz]=size(rr(r0~=0));

rrr=round(rr(r0~=0))

aaa=aa(r0~=0)

bbb=bb(r0~=0)

pz2=pz+k11-1;k9=1;

for k8=k11:pz2

RR(k8)=rrr(k9);

XX(k8)=aaa(k9);

YY(k8)=bbb(k9);

k9=k9+1;

end

k11=k8+1;

% if rr==zeros(1,n);

% as=as+1; continue;

% end

end

% if aa==zeros(1,300)

% as=1

% end

rangex=1:1:range;rowy=1:1:row;[RX,RY]=meshgrid(rangex,rowy);

RT=RX;figure;

mesh(RX,RY,RT);hold on;

plot3(XX,YY,RR,'ko');%%%%%%%%%%%%%%%5

[rhao,nn]=the_max1(XX,YY,RR);

i1=ones(row,range);

rhao

center=rhao;

% t2=1.5;t3=1.5;

for k1=1:1:row

for k2=1:1:range

XXX=k2-1;YYY=row+1-k1;

if abs(((XXX-center(1))*(XXX-center(1))+(YYY-center(2))*(YYY-center(2)))-center(3)*center(3))<=101

i1(k1,k2)=0;

end

imwrite(i1,'E:\MATLAB6p5\work\找出的圆.bmp','bmp');

i2=imread('E:\MATLAB6p5\work\找出的圆.bmp');

figure;

imshow(i2);st

Hough 변환_허프만방법으로 원, 직선을 검측하기(01.circle_product.m)

2010-10-22T10:32:00.001+09:00

% clear all; close all; clc;

function [ncir,st]=circle_product

ncir=ones(300,300);

[m,n]=size(ncir);

center=[0 0];%x=(b-a)*r + a

center(1)=round(((0.6-0.5)*rand+0.5)*m);%y

center(2)=round(((0.6-0.5)*rand+0.5)*n);%x

ol=ones(1,5);ol(1)=abs(center(1)-1);ol(2)=abs(center(1)-m);ol(3)=abs(center(2)-1);ol(4)=abs(center(2)-n);

for i=1:1:4

if ol(i)<=ol(i+1)

ol(5)=ol(i);

ol(i)=ol(i+1);

ol(i+1)=ol(5);

end

ol(5)=ol(i+1);

end

r=round(ol(5)*(3/4));

for k1=1:1:m

for k2=1:1:n

xxx=k2-1;yyy=m+1-k1;

if abs(((xxx-center(2))*(xxx-center(2))+(yyy-center(1))*(yyy-center(1)))-r*r)<=101

ncir(k1,k2)=0;

end

imwrite(ncir,'E:\MATLAB6p5\work\circle.bmp','bmp');

i2=imread('E:\MATLAB6p5\work\circle.bmp');

figure;

imshow(i2);

st=[center(2),center(1),r]

MATLAB으로 원을 그리고 N등분하기

2010-10-22T10:25:00.001+09:00

clc;clear R=1;t=0:pi/20:2*pi; x=R*cos(t);y=R*sin(t); plot(x,y),axis equal n=48;a=2*pi/n; for k=0:n-1 hold on plot([-cos(pi+k*a),cos(pi+k*a)],[-sin(pi+k*a),sin(pi+k*a)]) end

높이가 다른 동심원 그리기(3차원보기)

2010-10-22T10:20:00.001+09:00

r=input('반지름 r입력：'); t=-pi:pi/50:pi; hold on for i=1:r     z=i.*t./t;     x=i.*sin(t);     y=i.*cos(t);     pause(1);     plot3(x,y,z); end hold off  view(3)

동심원 그리기(소스코드)

2010-10-22T10:14:00.001+09:00

clear;clc

r=input('원의 반지름 r을 입력하세요(예, 5): ');
n=input('그릴 동심원의 개수를 입력하세요.(예, 30): ');
a=input('동심원 원의 중심의 x축 좌표를 입력하세요(예, -2):');
b=input('동심원 원의 중심의 y축좌표를 입력하세요(예, 3):');
t=linspace(-pi,pi);
x=sin(t)'*linspace(0,r,n+1)+a;
y=cos(t)'*linspace(0,r,n+1)+b;
plot(x,y)
axis equal

walsh_demo관련 결과데이터

2010-10-11T17:38:00.001+09:00

변환과정 보기

--
연변대학교 경제관리학원 정보관리&정보시스템학과
延边大学经济管理学院信息管理与信息系统

walsh_demo(하다마드, 왈쉬, 오리지날이미지, 변환이미지, 역변환이미지)

2010-10-11T17:34:00.001+09:00

% 장파작성, 최공 정리

RGB=imread('(Gray_256)lena.jpg');

I=RGB;

figure(1);imshow(I);title('오리지날');

F=double(I);

J=hadamard(256);

figure(2);imshow(J);title('하다마드');

H=walsh_xulie(256);

W=H*F*H;

figure(3);imshow(H);title('왈쉬코드');

figure(4);imshow(W);title('변환이미지');

N=H*W*H;

figure(5);imshow(N,[]);title('역변환이미지');

--
연변대학교 경제관리학원 정보관리&정보시스템학과
延边大学经济管理学院信息管理与信息系统

하다마드(Hadamard)에서 왈쉬코드(Walsh Code)만들기(Walsh_xulie)

2010-10-11T17:29:00.001+09:00

%주: 장파에게 내준 숙제내용... 최공 부분 수정

% 해당 내용을 Function형 M파일로 저장하여 실행

function X=walsh_xulie(N)

N=ceil(log2(N));

H=hadamard(2^N); %function provided by Matlab to generate Hadamard martrix

Hp=zeros(2^N);

Ep=zeros(2^N);

Wg=zeros(2^N);

Eg=zeros(2^N);

E=eye(2^N);

%matrix with Paley order

for n=1:2^N

p=0;

for i=1:N

p=bitset(p,i,bitget(n-1,N-i+1));

end

Hp(p+1,:)=H(n,:);

Ep(p+1,:)=E(n,:);

end

Hp; % Matrix Hp is the matrix with Paley order

Ep; % the transform matrix for transforming H matrix to Hp matrix

% Matrix with Walsh order

for n=1:2^N

g=0;

% if (bitget(n-1,N)==0)

g=bitset(g,N,bitget(n-1,N));

%else

%g=bitset(g,N,1);

%end

for i=2:N

if (bitget(n-1,N-i+1)==0)

g=bitset(g,N-i+1,bitget(g,N-i+2));

else

g=bitset(g,N-i+1,1-bitget(g,N-i+2));

end

Wg(g+1,:)=Hp(n,:);

Eg(g+1,:)=Ep(n,:);

end

Wg; % Matrix Wg is the matrix with Walsh order

Eg; % the transform matrix for transforming Hp matrix to Wg matrix

%matrix with inversed walsh order

for n=1:2^N

m=0;

for i=1:N

m=bitset(m,i,bitget(n-1,N-i+1));

end

Hm(m+1,:)=Wg(n,:);

Em(m+1,:)=Eg(n,:);

end

Hm; % Matrix Hm is the matrix with inversed Walsh order

Em; % the transform matrix for transforming Wg matrix to Hm matrix

X=Wg;

논문작업: 2D Barcode(dot code A) 테스트 sample만들기

2010-10-08T16:28:00.001+09:00

파라메터

DPI120 dpi 가로54 dots 세로 54 dots

--
연변대학교 경제관리학원 정보관리&정보시스템학과
延边大学经济管理学院信息管理与信息系统

바코드 중국전문가 블로그

2010-10-07T19:52:00.001+09:00

http://blog.csdn.net/changemyself

대만 국립교통대학공업공정관리교수 사이트

2010-10-07T19:35:00.001+09:00

http://sparc.nfu.edu.tw/~chh/nctuiem.htm

바코드 소개 문서 다운로드하기(특허기준, 설명서)

2010-10-07T19:15:00.001+09:00

http://www.adams1.com/stack.html

2D 바코드 만드는 프로그램(barcode-generator)

2010-10-07T16:57:00.001+09:00

http://www.glinttech.com/downloads/barcode-generator/Setup.exe

--
연변대학교 경제관리학원
정보관리&정보시스템학과 최기철

延边大学经济管理学院
信息管理与信息系统崔基哲

1차원, 2차원 바코드 만들기 사이트

2010-10-07T16:48:00.001+09:00

어떤 코드를 만들지는 이미지를 참고하기 바람.

사이트: http://www.hdf.com.pl/exec/barcode.php?lang=en

1D and 2D Barcode(1차원, 2차원 바코드) Generator

2010-10-07T16:16:00.001+09:00

1D Barcode Symbologies

2D Barcode Symbologies

http://www.racoindustries.com/barcodegenerator/

匿名连接技术：洋葱路由技术(Onion Routing)

2010-10-06T18:08:00.000+09:00

洋葱路由技术(Onion Routing)[21]是美国海军研究实验室的研究者们采用多次混淆的办法，提出的一种新的匿名通信技术，洋葱路由是为了阻止在公用网络上进行窃听和流量分析，以提供双向、实时的匿名连接，可以在公开的计算机网络中隐蔽网络的结构，对在互连网上进行的跟踪、窃听和流量分析有很强的抵抗作用。

通信双方用洋葱包代替通常的TCP/IP 数据包，利用代理技术实现与目标系统间的连接。通信过程可以简单地描述如下：主机A 与主机B 的路由路径以次是路由器W、X、Y、Z。当主机A 要与主机B 进行通信时，主机A向代理路由器W 发送一个数据包，代理路由器对该数据进行封装，封装是按分层进行：首先是对代理路由器Z 的地址和要发送的数据进行封装，把刚封装好的包与路由器Y 的地址再进行二次封装，然后再把所得的包与路由器X 的地址进行封装，这样封装好的数据包类似于洋葱的结构，故称为洋葱包。

当代理路由器W 把这个洋葱包传到路由器X 时，路由器X 用自己的密钥对所得的数据包进行解密，得出下一站是路由器Y，于是路由器X 就剥掉包中指明Y 的地址及包头信息，并对洋葱包进行填充，使洋葱包的大小不变，然后把处理后的洋葱包传到下一个路由器Y，路由器Y 收到了洋葱包后，按照同样的原理把洋葱包又给下一个洋葱路由器，最终传到终点路由器，即代理路由器Z。Z 用自己的密钥解密洋葱包，得到主机B 的地址后，直接把数据包发送给主机B。

窃听者如果监听链路中的某个节点，如路由器 Y，只能得到洋葱包是从路由器X 发来的，以及下一站是路由器Z 的信息，不可能获得到其它任何信息。在节点转发时可以加入填充技术，剥掉最外层的洋葱连同产生的随机数填充组成与原来同样大小的洋葱数据包，这样在中间任何一个路由器上看到的信息都不相同，但包的大小相同，从而加大了攻击者流量分析的难度。通过对洋葱路由技术的进一步研究表明，洋葱路由技术很好的实现匿名通信。美国政府已于2001年7月24日给Onion Routing 授予了专利[8]。

--
연변대학교 경제관리학원 정보관리&정보시스템학과 최기철
延边大学经济管理学院信息管理与信息系统崔基哲

论文阅读：匿名数字指纹的发展方向

2010-10-06T17:34:00.001+09:00

匿名数字指纹技术是一门涉及密码学、图形学、视觉科学、中文信息处理、图像处理、通信、信息安全等交叉的边缘学科。近五六年来，匿名数字指纹的研究得到了国外研究界的重视，并已获得部分成果；而在国内，据公开的资料看，匿名数字指纹的研究尚未得到充分重视。整体看来，数字指纹的研究仍然处于起步阶段，其应用还有很大的局限性。

以下几个是匿名数字指纹技术中值得研究的

方向：

1、结合数字水印嵌入技术的发展，进一步提出更合理的数字指纹的"嵌入假设"，并通过对用户攻击手段的研究，提出更好的数字指纹编码方案；并且设计出既有较好抗合谋攻击能力，又有较高实现效率的非对称指纹方案。

2、对特殊的密码学协议进行研究，研究有效的匿名指纹体制的构造方法。

3、将特殊的数字签名体制用于匿名数字指纹体制的构造，或者结合匿名通信技术构造匿名指纹体制。

4、将数字指纹与数字水印技术、加密技术、数字签名技术和访问控制技术相结合，将其应用于实际的不同形式内容保护。

--
연변대학교 경제관리학원 정보관리&정보시스템학과
延边大学经济管理学院信息管理与信息系统

论文阅读：匿名数字指纹的研究现状和问题

2010-10-06T17:30:00.001+09:00

自 B.Pfitzmann 于1997 年被提出匿名数字指纹的概念以后，对其的研究得到了各方面的重视。在文献[4]中，B. Pfitzmann 给出了实现匿名数字指纹体制的模型，并讨论了几种变形；同时给出了在某些假设下实现匿名指纹的一种框架，但该框架过于依赖零知识证明而不具备实用价值。

J. Domingo-Ferrer 于1998 年提出了一种发行商一旦发现重新分发的拷贝，无需被怀疑用户及登记中心的参与，就能确定非法分发者的匿名指纹体制[11]，该体制的安全性基于计算离散对数的困难性及安全多方计算的可行性。

1999 年Dominigo-Ferrer 在文献中[12]利用Oblivious Transfer 协议来设计非对称匿名指纹体制，其核心思想就是发行商不能得到用户参与发行的拷贝，基于Oblivious Transfer 协议的安全性，不可能得到用户所获得的拷贝信息，发行商发现非法拷贝时，首先应提取出非法拷贝中的指纹信息，获得相应证据后，在登记中心的帮助下来确认用户身份。

文献[13]和[14]中，采用类似于文献[11]中的匿名实现策略，但是对于部分协议的实现效率做了改进，并且实现都是基于多方安全计算，不切实可行。

不过在文献[10]中，A. Sadeghi 指出了文献[13]和[14]以及类似的匿名协议实际上并不能提供真正的匿名性。因为登记中心和发行商如果进行联合，可以将他们就一个用户的有关记录进行对比和联系从而确定出用户的身份，因此成为是准匿名的。

首个具有严格意义上匿名性的指纹方案设计仍是由B.Pfitzmann 在文献[15]提出的。该文中，作者基于电子货币的不可跟踪性，将用户的分发过程视为重复花费设计了一种匿名指纹体制。该体制使用两类比特承诺方案，一类是基于离散对数的位承诺方案，用于对嵌入内容进行承诺并用来构造证明嵌入内容合理性的零知识证明协议；另一类基于二次剩余的位承诺方案用于生成带指纹拷贝的步骤。加入零知识证明协议使得后者带指纹拷贝的生成效率很低，因而在实际中不实用。另外在该体制中，用户需要参加审判协议才能最终确定其是否有罪。

文献[16]在该方案的基础上进行了改进，使用户不需要参与审判协议，但带指纹拷贝的生成方法仍然沿袭了文献[15]中的方法。

2001 年，M.Kuribayashi 和H.Tanaka 利用具有同态性质的Okamoto-Uchiyama 加密体制提高了文献[15,16]中指纹嵌入协议的效率[17]。但该文中对匿名性的讨论并不完整，并未真正讨论具体的登记协议。

最近，中科院冯登国、孙中伟等人提出基于加同态性质的Bresson 公钥密码体制的匿名数字指纹方案[18]，这种方案是一种避免使用安全多方计算的思路。该文没有设立独立运行的注册机构RA，已经假设在PKI 体系结构中注册管理的职能由CA 来完成。其基本假设就是：协议中引入的参与实体越多，效率就越低。经验证明，这种假设是成立的，但该体制仍沿袭了M.Kuribayashi 和H.Tanaka 的方法，其主要贡献在于提出了效率较高的具体的非对称指纹嵌入机制。

J.Camenisch 在文献[19]中利用群签名构造了一种匿名指纹体制，用户的身份是在群签名的签名群体中，因此可以被有效隐藏。发行商

可以从非法拷贝中提取出用户的秘密信息，并行使撤销管理人(revocation manager)的职能。如果使用的群签名体制中撤销管理人能够不需成员管理人(member-ship manager)的帮助即可实现用户身份确认，则发行商也可以不需登记中心的帮助直接实现对用户身份的跟踪。

从以上发表的文章看，数字指纹的研究得到了国外研究界的重视，并已获得部分成果。但总体来说匿名数字指纹的研究还处于理论摸索状态，实际应用的例子不多，有必要加大研究力度，使基于数字指纹的版权保护技术早日走向应用。然而，随着版权保护意识的不断深入人心，数字指纹这一新型版权保护手段将会得到越来越多的关注和重视。

--
연변대학교 경제관리학원 정보관리&정보시스템학과
延边大学经济管理学院信息管理与信息系统

论文阅读：匿名数字指纹的基本概念

2010-10-06T17:21:00.001+09:00

匿名信件、匿名投票等匿名行为在现实生活中早已存在。不同的学者对匿名性(anonymity) 、无关联性(unlinkability)和假名(pseudonymity)等概念有着不同的理解，因而这三个概念也经常被混淆[8]。Marit Koehntopp 和Andreas Pfitzman 于2000 年指出标准化定义这些概念的必要性[9]，并给出以下定义：

匿名性(anonymity)：一个行为所对应的实体是匿名的，是指对应该行为的实体在特定的、具有一定相同特性的实体集中的不确定性，这个集合称为匿名集，匿名性的强度与该集合的大小、各集合元素的可能概率相关。

关联性(linkability)：是指系统中的对象之间在某一行为发生前与发生后的联系。匿名性与关联性是紧密耦合的两个概念，匿名性往往以关联的形式表述。在一个匿名通信系统中的无关联性可分为消息对发送者的无关联性、消息对接收者的无关联性、消息对发送者与接收者联系的无关联性和消息之间的无关联性，对应的匿名性可表述为发送者匿名、接收者匿名、关系匿名和无关联性（unlinkability）。这里的无关联性是指不能判定两次不同的消息发送（接收）是否来自同一通信实体。

假名（Pseudonymity）：使用假身份来实现匿名的一种办法，在这种实现中，同一个实体的不同的行为有可能被发现是关联的，即可联接的。综合上所述的匿名的几个基本定义，匿名数字指纹应满足以下基本要求[4]：

1、正确性（Correctness）：只要协议的各个部分都能够成功执行，最后用户将得到其要购买的拷贝。

2、抗陷害性（Exculpability）：无辜的用户不会受到陷害。

3、可跟踪性（Traceability）：当发现非法拷贝时，发行商能够凭借其中的指纹信息对用户进行跟踪。

4、匿名性 (Anonymity)：如果没有拿到用户非法分发的拷贝，发行商(即使与TTP 联合)也不能确定用户的身份。

5、无关联性(Unlinkability)：发行商(即使与TTP 联合)也不能将同一个用户的不同购买进行联系。

匿名数字指纹中匿名的含义有两种：用户在购买带有自己指纹的拷贝时可以对自身身份进行保密；但如果发行商和TTP 进行联合则可以确定出用户的身份，即这种意义上的匿名不能抵抗发行商和TTP 等实体的合谋，因此称为准匿名的[10](semi-anonymous)，这种匿名的含义较弱。如果即使发行商和TTP 或其他实体进行联合，也不能确认无辜用户的身份，并且也不能对同一用户的不同购买进行联系，那么这种匿名可以认为是较强意义的匿名。

我们把匿名数字指纹协议中可能的角色描述如下：

1、数据发布者(Data supplier 或Merchant，简称M )：提供数字产品的商家。

2、合法用户(User 或Buyer，简称B )：购买数据发布者产品的正版用户。

3、盗版者(Traitor，简称T )：向非法用户提供个人密钥(指纹)的一个合法用户。

4、共谋者(Traitors)：一组通过比较各自个人密钥进行侵权活动的盗版者。

5、非法使用者(Pirate user)：使用traitors 提供的个人密钥获得发布者产品的盗版用户。

6、仲裁者(Judge，简称J )：仲裁merchant 对traitor(s)的起诉。

--
연변대학교 경제관리학원 정보관리&정보시스템학과
延边大学经济管理学院信息管理与信息系统

论文阅读：匿名指纹模式的必要性

2010-10-06T17:14:00.001+09:00

安全和匿名业务的需求使得匿名技术的研究逐步深入，由于匿名指纹模式既能保护用户的隐私，又能使买卖双方的权益得到保障，因此，它成为目前数字多媒体版权保护问题的一个研究热点[3, 5]。

--
연변대학교 경제관리학원 정보관리&정보시스템학과
延边大学经济管理学院信息管理与信息系统

论文阅读：数字水印和数字指纹的区别

2010-10-06T17:08:00.001+09:00

随着多媒体技术和计算机互联网络的飞速发展，人们获取数字信息已变得越来越便利和快捷，但随之出现的对数字多媒体作品的版权保护问题也日显突出。早期，人们通过加密、访问受限等方法保护它们的版权，可一旦恶意的用户获得了这些数据，就无法阻止其进行非法复制。数字水印技术是近几年发展起来的新型数字版权保护技术。数字水印技术通过在所销售的拷贝中嵌入与创作者有关的特定信息来到达保护版权的目的，可以用来证明一个数字产品的原始版权所有者是谁，但是不能够追踪到非法拷贝与传播者；数字指纹则是通过在所销售的拷贝中嵌入与购买者有关的特定信息（称为数字指纹），销售商可以在收缴到盗版拷贝时对非法拷贝者进行追踪[3]。

比较而言，数字水印是向数字产品中嵌入版权拥有者的一些信息，当发生争议时能够有效确认出版权归属，对相同的作品嵌入的水印信息是相同的。而数字指纹是在原产品中嵌入与用户有关的信息，产品提供者(也称发行商)能够根据该信息对非法用户进行跟踪。嵌入的内容对不同购买者是不同的。数字指纹技术作为内容保护的重要手段之一，与数字水印、数字签名、数据加密等技术都是相关非常紧密的。总的来说，数字指纹与数字水印技术区别如下：

1、关注点不同：数字水印技术研究主要关注一个算法在数字产品(主要针对数字图像、视频和音频)中如何嵌入一串信息(水印)并能够完整地检测恢复。而数字指纹主要关注的是发现盗版拷贝后，通过数字指纹技术，能提供不可否认的证据。

2、嵌入信息的目的不同：数字水印技术总嵌入的信息可以是发布商的名称，也可以是任意的有特定意义的BIT 串。而数字指纹嵌入的指纹信息，对每个发行的拷贝来说，都是不一样的，而嵌入的主要目的，是为了跟踪盗版拷贝。

3、嵌入信息的要求不同：数字水印嵌入的要求是不可见性、鲁棒性，还要求能抵抗一定的主动攻击。而数字指纹技术对嵌入的信息要具有合谋容忍性[6, 7]，能够抵抗一定的合谋攻击。

[1] Wagner N R. Fingerprinting. In: Proceedings of the 1983 Symposium on Security and Privacy. 1983, 18-22

[2] Douglas R S. Cryptography: Theory and Practice, Third Edition. 1st ed. Chapman & Hall/CRC, 2005: 616

[3] 吕述望, 王彦, 刘振华. 数字指纹综述. 中国科学院研究生院学报, 2004,21(03): 289-298

[4] Pfitzmann B, Waidner M. Anonymous Fingerprinting. In: Advances in Cryptology-EUROCRYPT '97, vol. 1233 of LNCS. Berlin/Heidelberg:Spinger-Verlag, 1997, 88-102

[5] 颜浩 , 陈克非. 数字指纹协议的研究与发展. 计算机科学, 2002, 29(12):146-149

[6] Wu M, W T, Wang Z J, et al. Collusion-resistant fingerprinting for multimedia.Signal Processing Magazine, IEEE, 2004, 21(2): 15-27

[7] Boneh D, Shaw J. Collusion-secure fingerprinting for digital data. IEEE Transactions on Information Theory, 1998, 44(5): 1897-1905

[8] 王继林, 伍前红, 陈德人等. 匿名技术的研究进展. 通信学报, 2005, 26(02):112-118

[9] Andreas P, Marit K. Anonymity, Unobservability, and Pseudonymity-A Proposal for Terminology. In: International workshop on Designing privacy enhancing technologies: design issues in anonymity and unobservability, vol. 2009 of LNCS. Berlin/Heidelberg: Springer-Verlag, 2001, 1-9

--
연변대학교 경제관리학원 정보관리&정보시스템학과 최기철
延边大学经济管理学院信息管理与信息系统崔基哲

论文阅读：数字指纹方案可以分为三种基本类型

2010-10-06T16:59:00.001+09:00

数字指纹方案可以分为三种基本类型，它们分别是对称指纹模式、非对称指纹模式和匿名指纹模式[3]。

1、对称指纹模式：在这种模式中，发行商通常将带有指纹的拷贝发放给用户。由于发行商和用户都知道该拷贝以及该拷贝中的指纹，当发现被非法分发的带有某用户指纹的拷贝时，将无法确定谁应该对它负责。因为拷贝可能是该用户分发的，也可能是发行商本人分发的以对该用户进行陷害。

2、非对称指纹模式：针对这一问题，Pfitzmann 和Schunter 引入了非对称(asymmetric)指纹的概念[4]，类似于非对称的加密体制能够实现不可否认性，非对称指纹体制最主要的特点是实现非法用户的不可否认性。这其中又包括以下三个方面：

首先协议如果能正常执行，用户就可以得到含有其指纹的合法拷贝；

其次对发行商而言，在一定的合谋尺寸下[7]，发行商能够从非法拷贝中跟踪出至少一个非法分发者，同时能够提供证明用户有罪的不可伪造的证据；

最后对用户而言，无论合谋人数的多少，无辜用户不能受到陷害。

3、匿名指纹模式：无论是对称还是非对称指纹协议中，用户均需在购买过程中提交自己的身份信息，这破坏了购买过程的隐秘性。正是在这种背景下，Pfitzmann 和Waidner 在文献[4] 中提出了匿名数字指纹(Anonymous Fingerprinting)的概念。这种指纹机制中，用户在购买拷贝的过程中不会泄漏自己的身份信息，但如果用户进行非法分发活动，凭借非法拷贝中的信息，发行商可以识别非法者的身份。

匿名指纹协议的实现通常是引入一个可信第三方TTP[5]（Trusted Third Part），负责为用户的真实身份进行登记，同时为用户发放购买过程中需要的一些验证信息(通常是假名及其相应的证书)。

匿名指纹技术实际上和其他匿名技术是密切相关的，必须有相应的匿名技术作为前提，才能保证用户的私人信息不被泄露。本文将在后面介绍相关的匿名技术在匿名数字指纹中的应用。

--
연변대학교 경제관리학원 정보관리&정보시스템학과
延边大学经济管理学院信息管理与信息系统

教你电脑自动开机并播放音乐提醒你起床,并且晚上自动关机！

2010-09-14T15:44:00.000+09:00

使用电脑的自动开机,哇,好爽,并且音乐可以随便你定哦，一天都有好的心
　　.自动定时开机
　　此功能需要在BIOS中设定,并且你的主板的高级电源管理设置必须支持这个功能,不过现在大多数的主板全支持,放心使用吧.
　　步骤一:
　　启动系统按DEL,进入BIOS,然后选中"POWERMANAGEMENTFEATURES"
　　(电源管理设置),回车进入.在电源窗口中找到"RESTOREONAC/POWER
　　LOSS"项,按回车,然后通过方向键设置其值为"POWERON".
　　步骤二:
　　选中"RESUMEONRTCALARM"并按回车,接着通过方向键设置其值为"ENABLE"
　　步骤三:
　　选中"RTCALARMDATE"按回车,然后将其值设定为"EVERYDAY",表示每天都进行一样的操作,按照同样的方法,分别将"RTCALARMHOURS""RTCALARMMINUTE""RTMALARMSECOND"设置为"07""15""00",这表示开机时间为早上的"7:15:00",全部设置好后,按F10,保存并重新启动.
　　(由于不同的主板的BIOS不一样,有些主板的设置可能与上面的有一点点不同，但只要记住在电源管理中找到"WAKEUP"字样的就行了，因为这是设置定时开机的选项)
　　2自动登陆系统
　　电脑开机了，但是我们有些人设置了开机密码,不能自动登陆,怎么办?
　　我们可以通过修改注册表来实现自动登陆.
　　步骤一:
　　在运行中输入"REGEDIT",打开HKEY_LOCAL_MACHINE\SOFTWARE\MICROSOFT\WINDOWSNT\CURRENTVERSION\WINLOGON.
　　步骤二:
　　在右方的窗口中找到"DEFAULTUSERNAME"子键,将数值数据改为用户登陆帐户名,如果不存在该键值,就可以在右边的窗口中右键新建字符串,接着将"新值#1"重命名为"DefaultUserName"后进行修改.
　　步骤三:
　　在右侧窗口中找到DefaultPassword,如果不存在,可以新建,然后将其键值设定为你的密码
　　步骤四:
　　按照以上的方法,我们在右边的窗口中新建一个名为"AutoAdminLogon"的键值,将其值设为1,然后退出重新启动就可以实现自动登陆
　　二:计划任务帮助我自动打开音乐
　　1让音乐吵醒你的耳朵
　　如果希望早上开机后可以听到美丽的音乐,可以按照下面的步骤进行,
　　步骤一:
　　运行WINDOWSMEDIAPLAYER,挑选自己喜欢的音乐添加到播放列表里,然后在菜单中选择文件→播放列表另存为将当前播放的列表进行保存并命名,例如取名为"早上好"
　　步骤二:
　　依次点击开始,所有程序,附件,系统工具,任务计划,这时将会打开任务计划窗口,在窗口中双击添加任务计划
　　步骤三:
　　在弹出的窗口中单击下一步,在出现的窗口中点浏览,然后在随后打开的"选择程序以进行计划"的窗口中找到之前创建的"早上好.wpl"播放列表，点击打开按钮添加,并将执行这个任务设定为"计算机启动时"
　　步骤四:
　　点下一步,在出现的窗口中输入自动执行任务的登陆帐户的用户名及密码,最后点完成!
　　OK拉!这样你就可以每天早上自动听到美妙的音乐了!
　　3.自动关机
　　如果你工作很投入，晚上经常忘记休息,或者忘记关机,那么这可以帮助你
　　步骤一:
　　在"任务计划"窗口中双击添加任务计划，下一步,接着点浏览,在选择程序以进行计划的窗口中的C:\WINDOWS\SYSTEM32"目录下找到SHUTDOWN.EXE程序并点打开按钮.
　　步骤二:
　　将执行这个任务的时间设置为"每天",点下一步,将起始时间设置为"23:00"并点下一步,按照要求输入用户名及密码,在单击完成时,打开此任务计划的高级属性，并点完成
　　步骤三:
　　在弹出的SHUTDOWN窗口中点"任务",然后在"行"拦中输入"C:\WINDOWS\system32\shutdown.exe-s-t60"(其中-S表示关闭计算机,-T60表示60秒后自动关闭计算机,关机时间可根据自己的需要设置)点击应用,并重新输入用户密码，最后点确定!
　　(如果想取消自动关机,在运行中输入"Shutdown-a"命令回车就可以了
　　不同的电脑呢,有不同的设置方法,但都是同样在电源管理选项POWERMANAGEMENTFEATURES里面
　　如果选项里设置是用日期的,先0那就是天天都自动开机了
　　如果怕你的低音炮浪费电呢,就自己做一块小一点的功放板,电源接在主机的电源输出接口上面,这样就可以在电脑开机的同时也给功放供电,就不会再浪费电啦,如果不会做呢,花几块去买一块小功放回来就OK了.别忘了给功放板接上一个扬声器哦!

상명대학교 저작권보호학과(Copyright Protection Sangmyung Univ)

2010-09-12T12:48:00.001+09:00

소개

지식사회의 급성장으로 지적재산권 및 저작권 보호에 대한 관심과 필요성이 증대되고 있어 문화기술(CT)과 정보기술(IT)을 융합하여 저작권보호와 관련된 법, 정책, 경영, 정보기술에 대한 총체적인 지식을 갖춘 국제적인 저작권 보호 인재를 양성하고 저작권 보호 관련 기술 표준화를 선도하기 위하여 본 학과가 설립되었다. 글로벌 인재 양성을 위하여 본 학과는 학생들에게 영어와 한자능력을 기본적으로 요구하고 있으며, 실사구시를 추구하기 위하여 다양한 프로그램을 운영하고 있다.

※ 연혁
2009년 3월 1일

＊저작권보호학과 30명 신설 인가
＊김종원교수 부임

※ 교육목적 및 목표교육목적
오늘날 문화산업에는 저작권 및 저작권 보호에 관한 법과 다양한 정책, 기술 등이 요구되고 있으며, 그에 따라 저작권 및 저작권 보호에 대한 전문적인 지식과 문화산업에서 필요한 기획개발 능력을 겸비한 인재를 요구한다. 본 전공에서는 저작권 및 저작권 보호에 대한 폭넓은 지식을 바탕으로 이를 창의적으로 해석할 수 있는 능력과 아울러 문화산업에서 활용할 수 있는 법, 정책, 기술에 대한 국제적 수준의 전문적 능력을 갖추는 것을 교육목적으로 한다.

교육목표
저작권 관련 기술 표준화를 선도하며 저작권보호와 관련된 융합형 실용주의, 글로벌 수준의 전문 인재양성

※ 졸업 요건

저작권보호학과는 글로벌 실용인재 양성을 목표로 하고 있기 때문에 내부적으로 학기별 토플점수와 한자능력검정을 요구하고 있음.
토플은 iBT로 1학년 1학기 45점, 2학기 61점, 2학년 1학기 79점, 2학년 2학기 100점을 요구하고 있으며, 한자능력 검정은 1학년 때 2급, 2학년 때 1급을 요구하고 있음.

※ 졸업 후 진로

저작권보호학과를 이수한 이후의 진로는 법학전문대학원을 진학하여 법조계로 진출하는 경우와 저작권보호와 관련된 다양한 산업체에진출이 가능함. 저작권신탁단체나 온라인서비스제공자는 저작권보호 전문 인력을 필요로 하고 있으며, 국내 150여개의 온라인서비스사업자, 20여개에 달하는 저작권신탁단체, 법무법인 등에 취업이 가능함. 특히, “저작권관리업법”이 통과되면 “저작권관리사”제도가 신설되고 본 학과에서 저작권관리사를 취득하여 전문 인력으로서의 산업계 진출이 가능함.

포렌식 워터마킹(Forensic Watermarking)

2010-09-12T12:21:00.001+09:00

상영되고 있는 영화의 비디오 카메라 녹화를 통한 콘텐츠 유출을 방지하는 기술. 눈에 보이지는 않지만 녹화된 콘텐츠의 분석을 통해 어느 스크린에서 촬영되었는지, 필요하다면 극장 내 어느 위치에서 촬영되었는지 등까지도 추적이 가능하기 때문에 콘텐츠 유출의 경로까지 파악할 수 있다.

콘텐츠에 대한 접근을 사전에 통제하는 DRM을 대체하기 위해 새로 개발된 것이 워터마킹/포렌식마크 기반의 저작권 기술이다. 워터마크는 콘텐츠에 사람이 인지할 수 없는 저작권 정보를 삽입하고 검출기를 통해 삽입 정보를 식별하는 기술인데, 이러한 워터마크 기술에 구매자 정보나 유통경로 및 사용자 정보를 삽입하여 유포자와 배포경로를 추적할 수 있는 기술이 포렌식마크이다.

현재 UMG(Universal Music Group)는 워터마크를 삽입하여 음악 트랙을 제공한 판매자 정보를 관리하고 있는데, 2009년 8월 국내 전문업체 마크애니는 UMG에 향후 3년간 매년 70만 달러에 달하는 워터마크 기술 공급계약을 체결한 바 있다. 또한 IPTV, 지상파 방송국 등에서 포렌식마크 기술을 이용한 불법 콘텐츠 추적과 모니터링 시스템 구축을 준비하는 등 워터마크/포렌식마크 기술수요가 점차 증대되고 있는 상황이다.

종래 저작권 기술은 DRM 등 단일한 기술을 이용하여 콘텐츠를 개별적으로 보호하는 것이 일반적이었다면, 향후에는 디지털 워터마킹/포렌식마크 기술 등 다양한 기술들이 결합된 시스템으로 진화될 것으로 전망된다. 또한 IPTV, VOD, 스마트폰 등 새로운 콘텐츠 유통모델의 개발과 시장 확대로 인해 이에 필요한 기술을 선점하기 위한 경쟁은 더욱 심화 될 것으로 예상된다. DRM기술은 MS, Adobe등 글로벌 IT업체가 독점적 지위에서 시장을 선점하고 있어 우리나라에서는 많은 자금이 기술 로열티로 해외에 유출되고 있는 실정이다.

현재 우리나라의 저작권 보호·유통 기술수준은 미국의 83.4% 수준으로 3.8년의 기술격차를 보이고 있다. 기술 로열티에 의한 국부유출을 막고 저작권기술을 하나의 산업으로 육성하기 위해서 국가차원의 노력이 필요하다. 우리나라 정보통신기술은 세계최고 수준으로 기술개발에 필요한 기본인프라는 이미 구축되어 있어 국책 R&D사업으로 추진된 'DRM 상호 연동기술'과 같이 국가 차원에서 저작권기술을 적극적으로 개발·지원한다면 저작권기술 관련 산업이 국가경쟁력을 키울 수 있는 새로운 산업으로 발전될 것이다.

포렌식 워터마킹(Forensic Watermarking)

2010-09-12T12:19:00.001+09:00

상영되고 있는 영화의 비디오 카메라 녹화를 통한 콘텐츠 유출을 방지하는 기술. 눈에 보이지는 않지만 녹화된 콘텐츠의 분석을
통해 어느 스크린에서 촬영되었는지, 필요하다면 극장 내 어느 위치에서 촬영되었는지 등까지도 추적이 가능하기 때문에 콘텐츠
유출의 경로까지 파악할 수 있다.

콘텐츠에 대한 접근을 사전에 통제하는 DRM을 대체하기 위해 새로 개발된 것이 워터마킹/포렌식마크 기반의 저작권 기술이다.
워터마크는 콘텐츠에 사람이 인지할 수 없는 저작권 정보를 삽입하고 검출기를 통해 삽입 정보를 식별하는 기술인데, 이러한
워터마크 기술에 구매자 정보나 유통경로 및 사용자 정보를 삽입하여 유포자와 배포경로를 추적할 수 있는 기술이 포렌식마크이다.
현재 UMG(Universal Music Group)는 워터마크를 삽입하여 음악 트랙을 제공한 판매자 정보를 관리하고 있는데,
2009년 8월 국내 전문업체 마크애니는 UMG에 향후 3년간 매년 70만 달러에 달하는 워터마크 기술 공급계약을 체결한 바
있다. 또한 IPTV, 지상파 방송국 등에서 포렌식마크 기술을 이용한 불법 콘텐츠 추적과 모니터링 시스템 구축을 준비하는 등
워터마크/포렌식마크 기술수요가 점차 증대되고 있는 상황이다.

종래 저작권 기술은 DRM 등 단일한 기술을 이용하여 콘텐츠를 개별적으로 보호하는 것이 일반적이었다면, 향후에는 디지털
워터마킹/포렌식마크 기술 등 다양한 기술들이 결합된 시스템으로 진화될 것으로 전망된다. 또한 IPTV, VOD, 스마트폰 등
새로운 콘텐츠 유통모델의 개발과 시장 확대로 인해 이에 필요한 기술을 선점하기 위한 경쟁은 더욱 심화 될 것으로 예상된다.
DRM기술은 MS, Adobe등 글로벌 IT업체가 독점적 지위에서 시장을 선점하고 있어 우리나라에서는 많은 자금이 기술
로열티로 해외에 유출되고 있는 실정이다.

현재 우리나라의 저작권 보호·유통 기술수준은 미국의 83.4% 수준으로 3.8년의 기술격차를 보이고 있다. 기술 로열티에 의한
국부유출을 막고 저작권기술을 하나의 산업으로 육성하기 위해서 국가차원의 노력이 필요하다. 우리나라 정보통신기술은 세계최고
수준으로 기술개발에 필요한 기본인프라는 이미 구축되어 있어 국책 R&D사업으로 추진된 'DRM 상호 연동기술'과 같이 국가
차원에서 저작권기술을 적극적으로 개발·지원한다면 저작권기술 관련 산업이 국가경쟁력을 키울 수 있는 새로운 산업으로 발전될
것이다.

오디오 워터마킹 알고리즘: MP3, WAV에 워터마크 삽입

2010-08-23T16:11:00.001+09:00

% matlab을 이용하여 MP3 파일을 읽는 부분은 잠시 여기서 고민하지 않는다.

% 본 프로그램은 MP3를 WAV로 전환한 후 워터마크를 삽입한다.

% 이하 소스 코드임

% 시작

clc;clear;close all;

% 워터마크 삽입 프로그램

% 오디오 신호 3.wav를 읽어서 변수 A에 저장한다.

[A,fs,nbits] = wavread('3.wav');

% 오리지날 오디오 신호를 보여준다.

subplot(211); plot(A);

axis([0 350000 -2 2]);

title('오리지날 오디오신호의 파형');

% 변수 L에 오디오 신호의 길이를 저장한다.

L = size(A);

% 이미지 Lena.bmp를 변수 M에 저장하여 워터마크 신호로 사용한다.

M = imread('Lena.bmp');

% 이미지 Lena.bmp를 이진화하며 변수 BW에 저장한다.

BW = im2bw(M);

% 워터마크 매트릭스의 크기를 계산한다.

[M1,M2] = size(BW);

% M12은 중간 변수이다. 번마다 M1*M2를 계산하는 작업량을 줄인다.

M12 = M1*M2;

% 차원을 낮춘다(2차원->1차원), 얻은 1차원 워터마크

%신호를 C에 저장한다.

C = reshape(BW,1,M12);

n = M12;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% 얻은 워터마크 신호에 대하여 스프레드 스펙트럼 처리를 진행함

% 효과는 그닥잖음, 나중에 수정해야 됨.

% 스프레드 스펙트럼 처리계수는 2로 한다.

n = M12*2;

M = zeros(n,1);

% 비밀키 열을 생성M

for k = 1 : n

if mod(k,4) == 0

M(k) = 1;

else

M(k) = 0;

end

% 워터마크 열을 비밀키 열과 비트XOR를 함

l = ceil(k/2);

S(k) = bitxor(C(l),M(k));

end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% 워터마크 하나를 삽입하기 위한 오디오 데이터를 N이라 함

N = 10;

length = n*10;

% 오리지날 오디오 신호를 2개 부분으로 분해한다. 즉 Ae 와 Ar

i = 1 : length;

j = [1];

% 매트릭스 A의 l부터 length행까지로 매트릭스 Ae를 만든다.

Ae = A(i,j);

% Ae(i,j)

i = length+1 : L;

% 매트릭스 A의 length부터 L 제1렬로 매트릭스 Ar를 만든다.

Ar = A(i,j);

% 셀B를 만듬，오디오 데이터 Ae(m)는 B의 하나의 원소임

k = 1;

% 행수가 M1 x M2이고 열수가 l인 셀을 만듬

B = cell(n,1);

th = n*N;

% k가 Ae의 길이보다 작을 때, 매트릭스 Ae의 매 10행마다

% 오디오 데이터 단위로 셀B에 저장한다.

while ( k < th )

i = k : k+9;

m = (k+9)/10;

B{m,1} = Ae(i,j);

k = k+10;

end

% 셀 D를 만듬. 그리고 셀B의 이산코사인변환(DCT)한 원소들을 D에 저장함.

D = cell(n,1);

for i = 1 : n

D{i,1} = dct(B{i,1});

end

% 셀E를 만듬, E에는 워터마크 신호가 삽입된 것을 저장(中频系数)

E = cell(n,1);

E = D;

for i = 1 : n

% % 워터마크 정보 C를 오디오 데이터에 삽입

% E{i,1}(3) = (D{i,1}(3))*(1+2*C(i));

% 스프리드 스펙트럼 처리를 거친 워터마크 정보S를 오디어 데이터에 삽입

E{i,1}(3) = (D{i,1}(3))*(1+2*S(i));

end

% 셀 F를 만듬

F = cell(n,1);

% 셀E의 IDCT된 원소들을 셀F에 저장함

for i = 1 : n

F{i,1} = idct(E{i,1});

end

% 쪼개서 처리했던 부분들을 통합하여 매트릭스 G에 저장한다.

G = F{1,1};

% i는 3으로부터 M1*M2까지 차례로 읽어드림

for i = 2 : n

% 셀G와 셀F의 제i조 원소를 묶어서 G에 적어넣는다.

G=[G; F{i,1}];

end

% 셀G와 행렬Ar를 묶어서 G에 적어넣는다.

G =[G; Ar];

wavwrite(G,fs,nbits,'test.wav');

subplot(212); plot(G);

title('워터마킹돈 오디오신호파형');

axis([0 350000 -2 2]);

runprogram.m(Barcode.jpg)

2010-08-16T15:38:00.001+09:00

I=imread('barcode.jpg');%read image from drive;
I=im2bw(I);
imshow(I);
hold on;
[x,y]=ginput(2);        %get data from click mouse
u1=x(1);
u2=y(1);
v1=x(2);
v2=y(2);
line(x,y);              %paint the line between 2 point at click mouse
a = [u1:1:v1];          %find out function line between 2 point
if (u1~=v1) & (u2~=v2)
    g = round(((v2-u2)/(v1-u1))*a + u2 -(u1/(v1-u1)));
elseif (u1==v1)
    a = u1;
end
i=1;                    %find the value each pixel that the line go through
h=v1-u1+1;
b=[1:1:h];
for a=u1:v1
    if (I(g,a)==0)
        b(i)=1;
        i=i+1;
    elseif (I(g,a)==1)
        b(i)=0;
        i=i+1;
    end
end
c=b
i=1;                %find number pixels of each line in barcode picture;
s=[1:1:60];
for k=1:60
    j=0;
    if c(i)==0
        while (c(i)==0)&(i<=h)
            j=j+1;
            i=i+1;
            s(k)=j;
        end
    elseif c(i)==1
        while (c(i)==1)&(i<=h)
            j=j+1;
            i=i+1;
            s(k)=j;
        end
    end
end
mau=s(2);       %the first line is the sample for barcode,in another line is ratio with this
q=s./mau;
p=round(q);
doc1=[1:1:6];   %decode
k=1;
for i=5:4:25
    if (p(i)==3)& (p(i+1)==2) &(p(i+2)==1) &(p(i+3)==1)
        doc1(k)=0
    elseif (p(i)==2)& (p(i+1)==2) &(p(i+2)==2) &(p(i+3)==1)
        doc1(k)=1
    elseif (p(i)==2)& (p(i+1)==1) &(p(i+2)==2) &(p(i+3)==2)
        doc1(k)=2
    elseif (p(i)==1)& (p(i+1)==4) &(p(i+2)==1) &(p(i+3)==1)
        doc1(k)=3
    elseif (p(i)==1)& (p(i+1)==1) &(p(i+2)==3) &(p(i+3)==2)
        doc1(k)=4
    elseif (p(i)==1)& (p(i+1)==2) &(p(i+2)==3) &(p(i+3)==1)
        doc1(k)=5
    elseif (p(i)==1)& (p(i+1)==1) &(p(i+2)==1) &(p(i+3)==4)
        doc1(k)=6
    elseif (p(i)==1)& (p(i+1)==3) &(p(i+2)==1) &(p(i+3)==2)
        doc1(k)=7
    elseif (p(i)==1)& (p(i+1)==2) &(p(i+2)==1) &(p(i+3)==3)
        doc1(k)=8
    elseif (p(i)==3)& (p(i+1)==1) &(p(i+2)==1) &(p(i+3)==2)
        doc1(k)=9
    end
    k=k+1;
end
doc2=[1:1:6];
k=1;
for i=34:4:54
    if (p(i)==3)& (p(i+1)==2) &(p(i+2)==1) &(p(i+3)==1)
        doc2(k)=0
    elseif (p(i)==2)& (p(i+1)==2) &(p(i+2)==2) &(p(i+3)==1)
        doc2(k)=1
    elseif (p(i)==2)& (p(i+1)==1) &(p(i+2)==2) &(p(i+3)==2)
        doc2(k)=2
    elseif (p(i)==1)& (p(i+1)==4) &(p(i+2)==1) &(p(i+3)==1)
        doc2(k)=3
    elseif (p(i)==1)& (p(i+1)==1) &(p(i+2)==3) &(p(i+3)==2)
        doc2(k)=4
    elseif (p(i)==1)& (p(i+1)==2) &(p(i+2)==3) &(p(i+3)==1)
        doc2(k)=5
    elseif (p(i)==1)& (p(i+1)==1) &(p(i+2)==1) &(p(i+3)==4)
        doc2(k)=6
    elseif (p(i)==1)& (p(i+1)==3) &(p(i+2)==1) &(p(i+3)==2)
        doc2(k)=7
    elseif (p(i)==1)& (p(i+1)==2) &(p(i+2)==1) &(p(i+3)==3)
        doc2(k)=8
    elseif (p(i)==3)& (p(i+1)==1) &(p(i+2)==1) &(p(i+3)==2)
        doc2(k)=9
    end
    k=k+1;
end

walsh transform(1D) 왈쉬 변환(소스코드)

2010-08-13T12:40:00.001+09:00

%The function implement the 1D Walash Transform which can be used in

%signal processing,pattern recognition and Genetic algorithms. The Formula of

%1D Walsh Transform is defined as :

%% N-1 q-1

%% ---- --

%% 1 \ | | b[i](m)*b[q-1-i](u)

%% W(u) = --- / f(m)| | (-1) ,u = 0,...,N-1

%% N ---- | |

%% m=0 i=0

%% where for instance ,A[i] is the ith indices of A.

%The definition :

%%% f(m) : One dimentional Image(Sequence) ,m = 0,...,N-1;

% q

%%% N = 2 : where N is the size of Image

%%% W(u) : Walsh Transform;

%% b[k](u): kth bit(from LSB) in the binary representation of u;

%% For instance if u = 6 where in binary it becomes 110 then

%% b[0](6) = 0,b[1](6) = 1 and b[2](6) = 1,

%** Example :

%% W = walsh1d([1 2 1 1]);

%% Author : Ahmad poursaberi

%% e-mail : a.poursaberi@ece.ut.ac.ir

%% Control and Intelligent Processing Center of Excellence

%% Faulty of Engineering, Electrical&Computer Department,

%% University of Tehran,Iran,August 2004

function W = walsh1d(I);

warning off

siz = size(I);

siz = siz(2);

q = log2(siz);

if sum(ismember(char(cellstr(num2str(q))),'.'))~=0

disp(' Warning!... ');

disp('The size of Vector must be in the shape of 2^N ..');

return

else

for u = 1:siz

binu = dec2bin(u-1,q);

binsize = size(binu);

binsize = binsize(2);

Wtemp = 0;

%%Inner Loop

for m = 1:siz

binm = dec2bin(m-1,q);

binsize = size(binm);

binsize = binsize(2);

temp = 1;

for i = 1:q

temp = temp * (-1)^(binm(q+1-i)*binu(i));

end

Wtemp = I(m)*temp + Wtemp;

end

W(u) = inv(siz)*Wtemp;

end

%%End of Loop

end

핑거프린팅: 평균공모공격 방식 구조도

2010-08-13T12:16:00.001+09:00

핑거프린팅: 평균공모공격 방식 구조도

핑거프린팅 시스템구조도

2010-08-13T12:08:00.001+09:00

핑거프린팅 시스템구조도

RSA 암호화알고리즘(Java내부함수를 MATLAB에서 사용)03fo03

2010-08-13T11:12:00.001+09:00

function varargout=rsakeys(varargin)

% RSAKEYS - rapidly generate RSA private and public keys up to 2048 bits

% [prikey,pubkey]=rsakeys

% [prikey,pubkey]=rsakeys(strength)

% [prikey,pubkey,details]=rsakeys(...)

% prikey = struct with fields "modulus" and "exponent" where the

% exponent has significant length, plus a field "intent" which

% is set to 'private' as a reminder

% pubkey = struct with fields "modulus" and "exponent" where the

% exponent is always 65537, and a field "intent" which is

% set to the string 'public' as a reminder

% strength = an RSA strength of 512, 1024, or 2048 bits (default = 2048)

% details = optional output string containing computational details about

% the key pair

% Notes: (1) This function uses Java, which contains the calculation

% routines

% (2) An RSA "key" consists of two quatities: a modulus and an

% exponent. RSA keys always come in pairs, with one key denoted

% "public" and the other as "private." A message that is

% encrypted using EITHER key must be decrypted using the OTHER

% key. You keep the private key secret and broadcast the public

% key. Anyone can then encrypt any message to you using the

% public key, but no one but you can decrypt it (not even the

% sender) since only you have the secret key. Also, you can

% encrypt a message using your private key, and anyone can

% decrypt it using the public key. The benefit is that they then

% know for sure that you are the person who generated the

% message.

% (3) RSA modulus and exponent selection have certain nuances

% beyond the basic numerical calculations. In particular,

% certain values for various components of the key generation

% can result in encryption weaknesses. For this reason, you

% may choose to avoid generating RSA key pairs yourself and

% instead rely upon "official" key generation code.

% (4) As of 2005, "RSA Laboratories currently recommends key

% sizes of 1024 bits for corporate use and 2048 bits for

% extremely valuable keys...." For more information, see

% http://www.rsasecurity.com/rsalabs/node.asp?id=2218

% (5) As of 2005, the US National Institue of Standards and

% Technology claims that 2048 bit keys for RSA will remain

% resistant to cracking until about 2030.

% (6) Tested but no warranty; use at your own risk.

% (7) Michael Kleder, Nov 2005

% EXAMPLE

% [pr,pu,d]=rsakeys(1024)

strength=2048;

if nargin > 0

strength=varargin{1};

if ~any(ismember(strength,[512 1024 2048]))

error('Supported RSA strengths are 512, 1024, and 2048 bits')

end

kpg=java.security.KeyPairGenerator.getInstance('RSA');

kpg.initialize(strength,java.security.SecureRandom.getInstance('SHA1PRNG'));

kp=kpg.generateKeyPair;

pubkey.modulus=char(kp.getPublic.getModulus.toString);

pubkey.exponent=char(kp.getPublic.getPublicExponent.toString);

pubkey.intent='public';

prikey.modulus=pubkey.modulus;

prikey.exponent=char(kp.getPrivate.getPrivateExponent.toString);

prikey.intent='private';

varargout{1}=prikey;

varargout{2}=pubkey;

if nargout > 2

varargout{3}=char(kp.getPrivate.toString);

end

return

RSA 암호화알고리즘(Java내부함수를 MATLAB에서 사용)02fo03

2010-08-13T11:11:00.001+09:00

function outp = rsaenc(inp,key)

% RSAENC - rapidly encrypt a short message using RSA

% ct=rsaenc(pt,key)

% pt = plaintext message to encrypt (a vector of type char or uint8)

% up to the following size:

% 245 bytes if using a 2048-bit key

% 117 bytes if using a 1024-bit key

% 53 bytes if using a 512-bit key

% key = private or public key created by the RSAKEYS function

% ct = cyphertext (uint8 vector)

% Notes: (1) This function uses Java, which contains the calculation

% routines

% (2) The usual use of RSA is for encryption of keys in a "key

% exchange," so the plaintext size limits provided by the Java

% routines (which are adequate for that purpose) have not been

% expanded here. (For example, a 128-bit AES encryption key

% requires only 16 bytes.)

% (3) An RSA "key" consists of two quatities: a modulus and an

% exponent. RSA keys always come in pairs, with one key denoted

% "public" and the other as "private." A message that is

% encrypted using EITHER key must be decrypted using the OTHER

% key. You keep the private key secret and broadcast the public

% key. Anyone can then encrypt any message to you using the

% public key, but no one but you can decrypt it (not even the

% sender) since only you have the secret key. Also, you can

% encrypt a message using your private key, and anyone can

% decrypt it using the public key. The benefit is that they then

% know for sure that you are the person who generated the

% message.

% (4) RSA modulus and exponent selection have certain nuances

% beyond the basic numerical calculations. In particular,

% certain values for various components of the key generation

% can result in encryption weaknesses. For this reason, you

% may choose to avoid generating RSA key pairs yourself and

% instead rely upon "official" key generation code such as the

% Java routines used in this function.

% (5) As of 2005, "RSA Laboratories currently recommends key

% sizes of 1024 bits for corporate use and 2048 bits for

% extremely valuable keys...." For more information, see

% http://www.rsasecurity.com/rsalabs/node.asp?id=2218

% (6) As of 2005, the US National Institue of Standards and

% Technology claims that 2048 bit keys for RSA will remain

% resistant to cracking until about 2030.

% (7) Tested but no warranty; use at your own risk.

% (8) Michael Kleder, Nov 2005

% EXAMPLE:

% [pri,pub]=rsakeys(1024);

% char(rsadec(rsaenc('This is a secret.',pub),pri))

% char(rsadec(rsaenc('This is a secret.',pri),pub))

if ischar(inp) | islogical(inp)

inp=uint8(inp(:));

else

inp=typecast(inp(:),'uint8');

end

modulus = java.math.BigInteger(key.modulus);

priexp = java.math.BigInteger(key.exponent);

kfac = java.security.KeyFactory.getInstance('RSA');

% For consistency, always use Java private key encryptor. We are not doing

% key intent verification here (we just want the math):

prikey=kfac.generatePrivate(java.security.spec.RSAPrivateKeySpec(modulus,priexp));

c=javax.crypto.Cipher.getInstance('RSA');

c.init(1,prikey);

ct=typecast(c.doFinal(inp),'uint8');

outp=typecast(ct,'uint8');

outp=outp(:)';

return

RSA 암호화알고리즘(Java내부함수를 MATLAB에서 사용)01fo03

2010-08-13T11:10:00.001+09:00

function outp = rsadec(inp,key)

% RSADEC - rapidly decrypt a short message using RSA

% pt=rsadec(ct,key)

% ct = cyphertext (uint8 vector) created by RSAENC

% key = private or public key created by the RSAKEYS function

% pt = plaintext message (a vector of type uint8) which you may

% wish to convert to text using "char(pt)"

% Notes: (1) This function uses Java, which contains the calculation

% routines

% (2) The usual use of RSA is for encryption of keys in a "key

% exchange," so the plaintext size limits provided by the Java

% routines (which are adequate for that purpose) have not been

% expanded here. (For example, a 128-bit AES encryption key

% requires only 16 bytes.) See the RSAENC function.

% (3) An RSA "key" consists of two quatities: a modulus and an

% exponent. RSA keys always come in pairs, with one key denoted

% "public" and the other as "private." A message that is

% encrypted using EITHER key must be decrypted using the OTHER

% key. You keep the private key secret and broadcast the public

% key. Anyone can then encrypt any message to you using the

% public key, but no one but you can decrypt it (not even the

% sender) since only you have the secret key. Also, you can

% encrypt a message using your private key, and anyone can

% decrypt it using the public key. The benefit is that they then

% know for sure that you are the person who generated the

% message.

% (4) RSA modulus and exponent selection have certain nuances

% beyond the basic numerical calculations. In particular,

% certain values for various components of the key generation

% can result in encryption weaknesses. For this reason, you

% may choose to avoid generating RSA key pairs yourself and

% instead rely upon "official" key generation code such as the

% Java routines used in this function.

% (5) As of 2005, "RSA Laboratories currently recommends key

% sizes of 1024 bits for corporate use and 2048 bits for

% extremely valuable keys...." For more information, see

% http://www.rsasecurity.com/rsalabs/node.asp?id=2218

% (6) As of 2005, the US National Institue of Standards and

% Technology claims that 2048 bit keys for RSA will remain

% resistant to cracking until about 2030.

% (7) Tested but no warranty; use at your own risk.

% (8) Michael Kleder, Nov 2005

% EXAMPLE:

% [pri,pub]=rsakeys(1024);

% char(rsadec(rsaenc('This is a secret.',pub),pri))

% char(rsadec(rsaenc('This is a secret.',pri),pub))

if ischar(inp) | islogical(inp)

inp=uint8(inp(:));

else

inp=typecast(inp(:),'uint8');

end

modulus = java.math.BigInteger(key.modulus);

pubexp = java.math.BigInteger(key.exponent);

kfac = java.security.KeyFactory.getInstance('RSA');

% For consistency, always use Java public key decryptor. We are not doing

% key intent verification here (we just want the math):

pubkey=kfac.generatePublic(java.security.spec.RSAPublicKeySpec(modulus,pubexp));

c=javax.crypto.Cipher.getInstance('RSA');

c.init(2,pubkey);

ct=typecast(c.doFinal(inp),'uint8');

outp=typecast(ct,'uint8');

outp=outp(:)';

return

RSA 암호화 알고리즘(MATLAB 소스코드)

2010-08-13T10:39:00.001+09:00

function RSA2=decod(miwen)

miwen=input('암호문을 입력하세요: ');

Q=input('Q=(p-1)(q-1)=');

e=input('e=');

n=input('n=');

k=0;

while mod(k*Q+1,e)~=0

k=k+1;

end

d1=(k*Q+1)/e; %계산 e의 모듈역

disp(d1)

l1=length(miwen);

b=dec2bin(d1);

l2=length(b);

for j=1:1:l1

c=0;

d=1;

for i=1:1:l2

c=2*c;

d=mod(d*d,n);

if b(i)=='1'

c=c+1;

d=mod(d*miwen(j),n);

end

if miwen(j)>=97&&miwen(j)<=122

minwen(j)=char(mod(d,26)+97); %암호문은 소문자

elseif miwen(j)>=65&&miwen(j)<=90

minwen(j)=char(mod(d,26)+65); %암호문은 대문자

else minwen(j)=d; %암호문은 숫자

end

disp(minwen)

%%%%%%%%%%%%%%%%%%%

function RSA1=encod(minwen)

minwen=input('평문을 입력하세요: ');

m=1;

for i=300:1:2000

k=0;

for j=2:1:(i-1)

if rem(i,j)==0

k=1;

break;

end

if k==0

f(m)=i;

m=m+1;

end

n1=10;

h=ceil(n1.*rand(1,2));%임의 두개의 숫자를 생성

p=f(h(1,1));

q=f(h(1,1)+h(1,2)); %두개의 서로 다른 소수를 생성

n=p*q; %공개키계산n

Q=(p-1)*(q-1);

for j=1:1:100000

k=0;

E=ceil(500*rand);

for i=2:1:(Q-1)

if rem(Q,i)==0&&rem(E,i)==0

k=1;

break;

end

if k==0

e=E;

break;

end

end %계산 공개키e

l1=length(minwen);

b=dec2bin(e);

l2=length(b);

for j=1:1:l1

c=0;

d=1;

for i=1:1:l2

c=2*c;

d=mod(d*d,n);

if b(i)=='1'

c=c+1;

d=mod(d*minwen(j),n);

end

if minwen(j)>=97&&minwen(j)<=122

miwen(j)=char(mod(d,26)+97); %평문은 소문자

elseif minwen(j)>=65&&minwen(j)<=90

miwen(j)=char(mod(d,26)+65); %평문은 대문자

else miwen(j)=d; %평문은 숫자

end

disp(Q) %Q=(p-1)(q-1)

disp(e)

disp(n) %비밀키(프라이빅 키)n

disp(miwen)

% 두개의 함수임 RSA1과 RSA2로 따로따로 저장하여 실행

%--

%연변대학교 경제관리학원 정보관리&정보시스템학과

%延边大学经济管理学院信息管理与信息系统

沃尔什编码（Walsh Code）

2010-08-12T18:32:00.001+09:00

沃尔什编码（Walsh Code）

沃尔什编码（Walsh Code）是一组具有优质自相关器的性能和劣质的交叉相关性能的传播代码。Walsh 编码是 CDMA
系统的支柱，它用于 CDMA 个人信道的研发。就 IS-95
而言，这里有64个代码可以使用。代码0作为引导信道，代码32作为同步信道，代码1到7作为控制信道，其他代码则作为通信信道。如果它们不被需要，代码2到7也可用来作通信信道。对
CDMA2000，存在着大量的 Walsh 编码，以其不同的长度来调节不同的数据比率和传布不同的无线结构因数。

Walsh Code(월쉬 부호, Walsh 코드, 왈쉬 코드, 왈시 코드)

2010-08-12T18:29:00.001+09:00

Walsh Code 월쉬 부호, Walsh 코드, 왈쉬 코드, 왈시 코드

1. 개요

ㅇ Walsh 부호는, 주로 CDMA 시스템에서 모든 사용자들이 동일한 주파수 대역을 사용하고 있기 때문에, 순방향 채널 상에서 상호간섭을 피하고 각각의 사용자 채널을 분리구분하기 위해 사용되는 직교성을 갖는 확산코드 세트를 말한다.
* 직교성 : 코드 간에 상호 상관성이 없음
2. 역사 및 성질
  ㅇ 1923년 왈쉬(J. L. Walsh)는 정규화 된 구간 (0,1)에서 +1과 -1만의 값을 가지는 완벽한 직교(Orthogonal) 함수계를 제시
ㅇ 성질
- Walsh 함수는 서로 다른 코드를 곱하면(Exclusive-OR), 0 (또는 -1) 과 1 이 섞여서 나오고 이를 모두 평균하면 0 이 되도록 되어 있고, 같은 코드를 곱하면 모두  1 이 나옴.
- 코드 간에 비교하면, 같은 비트와 다른 비트의 수가 같게됨
3. CDMA 방식 이동전화에의 응용
  ㅇ 용도 : 스펙트럼확산 및 채널구분
ㅇ IS-95A 방식에서 사용하는 Walsh 함수는 64 비트로 구성되어 64개 종류가 있다.
- 64x64 크기의 Hardmard 행렬에 의한 왈쉬코드를 사용하여 최대 64개의 물리적인 통신채널을 구분토록함
- 이를 W0, W1, W2, W3·· 로 구분한다.
. W0 : 파일럿 채널에 할당
. W1~7 : 페이징 채널에 할당 (W1이 기본채널에 할당, 나머지는 선택)
. W32 : 동기화 채널에 할당

4. 생성 방법
ㅇ Hardamard Matrix를 이용한 방법
ㅇ Walsh 함수의 대칭성을 이용한 방법

Orthogonal Signal, Orthogonal Code(직교 신호, 직교 코드)

2010-08-12T18:26:00.001+09:00

Orthogonal Signal, Orthogonal Code, Biorthogonal Code
직교 신호, 직교 코드, 직교 부호, 배직교 코드, 배직교 부호

1. 개요

ㅇ 직교 신호
- 그들의 내적(Inner Product)이 다음과 같이 항상 0 인 신호를 말함.
∞ *
< x₁(t), x₂(t) > = ∫ x₁(t) x₂(t) dt = 0
-∞
. (例) 직교 파형 신호의 예 : sin x 및 cos x
- 일정 시간 동안 직교 신호가 되는 기준
. 적어도 두 파형이 직교하려면 한 신호주기 동안 상호 상관 관계가 없을 것
T
∫ x₁(t) x₂(t) dt = 0
0
- 직교 신호 집합
. 상호 수직인 벡터 신호들로 형성되어지는 신호 집합.
. 시간상으로 전혀 겹침이 없는 신호 집합

ㅇ 직교 코드
- 서로 다른 코드 사이에 상호 상관도가 0 인 코드를 말함.
- 예) Walsh Code, OVSF(Orthogonal Variable Spreading Factor) 등

2. 응용 例
ㅇ 광대역 CDMA 방식은, 주파수대역확산이 PN 코드가 아니고 직교코드에 의해 이루어짐
- 동기식 CDMA 방식(cdma2000)
. 왈쉬코드(Walsh Code)로 64개 물리적 통신 채널을 구분
. 64개 이원 시퀸스(binary sequence)의 집합으로 구성
. 이 시퀸스는 Hamard 행렬을 사용하여 생성됨
- 비동기식 CDMA 방식(WCDMA)
. OVSF 기법에 의한 ...

Gold Code(골드 코드)

2010-08-12T18:23:00.001+09:00

1. 개요
Gold 코드는, 상관 특성이 좋고 코드의 종류가 다양한 코드로, 기지국 간에 동기를 맞추지 않은 비동기 방식인 W-CDMA 방식에서 기지국 식별용으로 주로 사용되는 코드를 말한다.

2. 참고사항
ㅇ Robert Gold 박사에 의해 개발됨.
ㅇ PN Code에 비해 매우 많은 수의 코드 종류를 만들 수 있다.

Scrambling Code(스크램블링 코드, 스크램블 코드)

2010-08-12T18:22:00.001+09:00

1. 개요
ㅇ 스크램블링 코드는, 비동기식 CDMA 방식에서 서로 다른 기지국이나 이동단말 등을 구분하려는 용도로 사용하는 코드에 대해
부르는 용어로써, 그 용도상 차이점으로는, 하향링크에서는, 기지국 셀/섹터 등을 구분하고, 상향링크에서는, 한 셀내에서
서로다른 이동국(단말)을 구분한다.
ㅇ 한편, 명칭상 차이로써, 동기식 CDMA 방식에서는 이를 PN 코드라고하여 m-sequence를 사용, 비동기식 CDMA
방식에서는 스크램블링 코드라고하며, Gold sequence(Gold 코드)를 사용
ㅇ 스크램블링은 확산 과정의 맨끝단에서 사용됨, 스크램블링 그 자체로는 Chip Rate에 변화 없음

PN Code, PN Sequence(PN 코드, 유사잡음코드, 유사잡음부호, 칩 코드)

2010-08-12T18:19:00.001+09:00

PN Code, PN Sequence Pseudo Noise Code, Pseudo random Noise Code, Chipping Code
PN 코드, 유사잡음코드, 유사잡음부호, 칩 코드

1. 개요
ㅇ 유사잡음코드(PN Code, Pseudo Random Noise Sequence)는, 랜덤 시퀀스와 유사한 잡음 특성을 보이면서도 재생이 가능한 코드 즉, 그 속에 일정한 규칙을 갖는 코드 시퀀스(수열)를 말함. Pseudo Randomness (의사 난수열)
2. PN 코드(유사잡음)로써 갖추어야 할 조건/특성
예리한 자기상관 특성 : 송수신 PN 코드가 동일하고 시간적 관계가 정확히 일치함
낮은 상호상관 특성: 사용자들에게 할당된 임의의 두 PN 코드사이에 상호상관특성이 작아야함
통계적 균형성 : 한 주기 코드에 "0" 과 "1"이 균형적으로 됨.
편이와 가산성 : 특정 PN 코드를 시간 지연 시켜 모듈러-2 연산하여 생긴 시퀸스는 본래 코드를 단지 시간 지연 시킨 코드와 동일한 코드가 됨.
런 특성 : 한 주기 안에 있는 심볼이 연속적으로 이어져 나오는 시퀸스를 의미함.

ㅇ 기타사항
발생의 용이성 : PN 발생기에서 긴 시퀸스를 쉽게 발생시킬 수 있어야 함
초기 동기화 용이성 : 이동통신 단말기 입장에서는 초기 동기를 재빨리 취할 수 있어야함
3. PN 수열 생성법
대부분의 PN 수열은 거의 항상 순서논리회로(피드백 쉬프트 레지스터 포함)를 이용
例) 최장 수열 생성기 (Maximum-length shift-register)
지금까지 가장 잘 알려진 PN 수열 발생기로, 간략히 m-sequence이라고도함
- 길이 : L = 2^m - 1 비트, 수열은 주기적이며 그 주기(길이)는 L
- 선형궤환(Linear Feedback)을 갖고있는 m 단계 쉬프트 레지스터에 의해 생성 즉, LFSR (Linear Feedback Shift Register) 카운터를 이용함
- 2^m -1 개의 상태를 갖음
- LFSR에 의한 수열의 상태 순서는 통상적인 이진 상태 계수 순서와는 다름
- 바로 이 점이 ‘의사 난수’계수에 의한 오류의 검출 및 수정(CRC 코드 등)에 유리하여, 많은 응용에 쓰임
- 한편, CDMA 방식의 PN 수열 발생기에도 쓰임
4. 응용의 例
ㅇ IS-95 등 CDMA 방식의 이동전화에서는 다음 2종류의 PN 코드를 사용
- Short PN Code : 기지국을 구분하기 위해 주로 사용 코드의 길이는 2의 15승 (32,768 chip)
. 순방향 채널 : 대역확산 및 기지국 구분 용도
. 역방향 채널 : 송신시점 결정 및 Long PN Code의 Mask로 사용
※ 현재 IS-95 방식에서는 기지국에 사용하는 PN 코드는 모두 동일한 코드를 사용하지만 서로를 구별하기 위해 기지국 마다 일정 시간 간격으로 Time Shift (Time Offset, PN offset) 되어 발생
- Long PN Code : 기지국이 이동국을 구분하고, 음성신호를 확산하기 위해 사용코드의 길이는 2의 42승
. 순방향 링크에서 암호화를 위해 사용
. 이는 단말기에 부여된 고유 코드를 응용하여 만들어짐

ㅇ 802.11 방식의 무선 LAN에서는 PN Code로 Barker Code (10110111000) 를 사용한다. 한편, 이를 칩 코드(Chipping Code)라고 부른다.

파일럿 채널, 파일럿 신호, 파일럿 부반송파, 파일롯 채널

2010-08-12T18:13:00.001+09:00

Pilot Channel, Pilot Signal, Pilot Subcarrier 파일럿 채널, 파일럿 신호, 파일럿 부반송파, 파일롯 채널

1. 개요
파일럿 채널이란 무선통신 기지국에서 이동국(수신단)으로하여금 반송파 위상 동기화, 기지국 정보 획득 등에 도움이 되도록 실제 데이터를 가지지 않으나, 항상 송출되는 높은 출력의 채널을 말한다.
2. 사용 例
ㅇ CDMA 방식 (IS-95 경우) 95 채널 구조

- 신호형태 : Pilot 채널은 변조되지 않은 스펙트럼 확산된 신호임
- 사용목적 : 이동국들이 초기 시스템 동기 및 기지국 신호의 위상을 추적하기 위함
- 특징
  . Pilot 신호는 이동수신기에서 타 채널 신호의 복조를 위한 반송파 위상 동기로 사용되므로, 다소 높은 레벨을 유지
  . Pilot 신호는 각 기지국 마다 서로를 구별하기 위해, 서로 다른 코드 위상 Offset (Pilot Offset)을 사용함
  . Pilot 기준신호에 맞추어 Coherent Detection을 가능하게 해줌
  ㅇ WiBro의 경우
- WiBro 수신기로 하여금 채널상태정보(CSI,Channel State Information)를 추정하고 추적하기 쉽게하기위해 삽입
- 파일럿용 부반송파 갯수 : 96 개
  ㅇ 802.11a 또는 HiperLAN/2
- 전송할 데이터에 사전에 정의된 4개의 부반송파 첨가

확산 이득, 확산처리이득, 프로세스 이득, 확산율

2010-08-12T18:06:00.001+09:00

Processing Gain, Process Gain, Spreading Ratio 확산 이득, 확산처리이득, 프로세스 이득, 확산율

1. 개요
확산이득은, 대역확산방식에서 원래 데이터 신호의 대역이 확산코드(Spreading Code)에 의해서 얼마나 넓게 스펙트럼이 확산될 수 있는 지를 나타내는 파라미터를 말함. 또는, 대역확산 시스템에서 점유 대역폭이 확산됨에 따라 얻어진 신호대 잡음비에서 얻어지는 이득을 말하기도 한다.

2. 확산이득 (Gp)

ㅇ 계산식 표현
Bss
  Gp = 10 log ----- (Bss : 확산 신호의 대역폭, Bs : 데이터 신호의 대역폭)
Bs
  - 확산 전송 대역폭 대 신호대역폭 의 比
  - 확산속도(Spreading Rate) 대 사용자 데이터 속도(초당 사용자 데이터 심볼수)

3. 확산율(Spreading Ratio) 또는 프로세스 이득(Processing Gain)
  ㅇ Processing Gain(확산이득) 및 Spreading Ratio(확산율)은 통상적으로 같은 의미로 사용되지만, 확산율은 대역확산 방식에서 하나의 정보 비트(심볼)를 전송하기 위해 사용되는 칩(Chip)의 개수를 말함.
- 높은 확산율은 복원 가능성이 더 좋지만, 대신에 더 넓은 대역폭이 요구됨
  ㅇ 확산이득의 경우는 이상적이 아닌 실제 시스템을 고려하기 때문에, 확산율 보다는 약간 더 낮음. ☞ 칩 코드, C/I
  ㅇ 例 : 802.11 규격에서 최소 확산이득 요구조건을 11로 정한바 있음

Spreading Code, Spreading Sequence 확산 코드, 확산 부호

2010-08-12T18:03:00.001+09:00

Spreading Code, Spreading Sequence 확산 코드, 확산 부호

1. 개요
확산대역 통신방식에서, 주로 대역확산용으로 정보신호에 곱해지는 코드를 확산코드라 한다.
2. 기술적 특징 및 참고사항
  ㅇ 확산 코드가 곱해진 신호는 원래의 정보 신호보다 대역폭이 훨씬 커지기 때문에, 코드를 곱하는 과정은 신호의 스펙트럼을 확산시키는 과정이며, 이러한 이유로 이 과정을 대역확산 변조라고도 한다.
ㅇ 결과적으로 최종 신호는 대역확산된 신호가 된다. 전송되는 신호를 대역확산시키고 각 사용자들을 부가적인 확산 코드를 이용하여
구분할 수 있기 때문에, CDMA 통신방식은 다중접속 (Multiple Access) 능력을 갖게 된다.
ㅇ 송신기에서는 이 확산된 신호를 반송파에 실려서 공중으로 전파시킨다.
ㅇ 수신기는 이 확산 코드를 미리 알고 있어 이 코드를 이용하여 수신된 신호를 복호하며, 원래의 정보를 복원할 수 있다.
- 이것이 가능한 이유는 원하는 사용자의 코드와 다른 사용자의 코드 간에 상호 상관(Cross-Correlation)이 작기 때문이다.
3. 주요 용도별 확산코드
  ㅇ 주요 용도들
- ① 대역 확산, ② 기지국/ 채널/ 사용자의 구분, ③ 데이터 암호화 등
  ㅇ 구분
- 주파수대역 확산용(Spread Spectrum) 및 다중 접속용(Multiple Access)
. 정보 신호에 확산코드가 곱해져 대역을 확산(Spread Spectrum)시키게 됨
. 例) IS-95 : 하향링크에 직교코드의 일종인 Walsh code 사용하여 대역확산
- 기지국 구분용 (일명 PN 코드라고 부름)
. 동기식 CDMA 방식 : PN 코드, m-sequence 사용
. 비동기식 CDMA 방식 : Scrambling 코드, Gold-seqence 사용

4. 확산 코드 구현 종별
  ㅇ PN Code : 통상 H/W에 의해 생성
- Long PN Code : 기지국이 단말기를 구분할때 사용됨
   - Short PN Code : 단말기가 기지국을 구분할때 사용됨
  ㅇ Walsh Code : 통상 정해진 코드 표(Table)를 통해 생성
- 순방향 채널에서 이동단말기가 각 채널을 구분할때 사용됨

CDMA관련: FH, FH-SS(주파수 호핑, 주파수 도약 스펙트럼)

2010-08-12T17:59:00.001+09:00

FH, FH-SS, FHSS Frequency Hopping, Frequency Hopping Spread Spectrum
주파수 호핑, 주파수 도약, 주파수도약변조, 주파수 도약 스펙트럼

1. 개요
주파수호핑 방식은 통신방식에 있어서 직접확산변조(Direct Sequence) 방식과 함께 대표적인 확산스펙트럼 변조방식을 말한다.

2. DS (Direct Sequence) 방식 및 FH (Frequency Hopping) 방식
ㅇ DS 방식은 디지털 부호계열(의사잡음계열)에 의해 반송파를 변조하는 것임에 비해, FH 방식은 부호계열에 의해 정해진 패턴으로 반송주파수를 불연속적으로 편이시키는 것이다. 이 방식은 100만개의 별도 주파수를 선택할 수 있어 각자 전송하는 정보와 부호에 근거하여 선택된다.

3. FH의 송수신 방식
ㅇ FH 방식은 우선 주어진 대역폭을 많은 수의 호핑 채널(Hopping Channel)로 나누고, 송신 측에서는 1차 변조된 신호(중간주파수)를 RF 주파수대로 변환할 때 미리 정해진 순서에 따라 서로 다른 호핑 채널에 할당시킨다. 이 때 신호가 할당된 채널이 빠르게 변경되기 때문에 다중 채널 간섭 등의 잡음의 영향을 줄일 수 있다.
ㅇ 한편, 수신단에서는 여러 호핑 채널에 분산된 수신신호들을 송신 측과 같은 순서로 연결함으로써 원 신호를 복원하게 된다.

4. 동작원리 및 특징
ㅇ 호핑 코드에 따라 순간적으로 한 주파수로 튜닝(Hopping)되어 그 주파수에서 아주 짧은 시간동안 전송하고 다시 다른 주파수로 튜닝되어 전과 동일한 과정을 반복하는 기법
ㅇ 송신부와 수신부는 같은 시간대에 같은 주파수에 위치토록 동기화가 반드시 필요함. 즉, 수신측은 신호를 수신하기 위하여 송신시 사용한 호핑 코드와 동일한 코드를 이용하여 특정시간에 특정주파수로 튜닝하여야 함
ㅇ 전파방해나 잡음간섭을 최대한 방지하고, 같은 주파수를 사용하더라도 호핑 코드만 다르면 여러 스프레드 스펙트럼 시스템을 동일 장소에 사용 가능

ㅇ 부품이 상대적으로 저렴하고 높은 출력을 요구하지 않음

CDMA관련: DSSS, DS-SS(직접확산, 직접 시퀀스 대역 확산)

2010-08-12T17:55:00.001+09:00

DSSS, DS-SS Direct Sequence, Direct Sequence Spread Spectrum
직접 확산, 직접확산변조, 직접 확산 스펙트럼, 직접 시퀸스 대역 확산

1. 개요
직접확산 변조란 원래의 신호에 주파수가 높은 디지털 신호(확산코드)를 곱(XOR)하여 확산(Spreading)시키는
대역확산(Spread Spectrum) 변조방식 을 말한다.

2. 동작방식 및 특징
ㅇ 전송될 2진 데이터 신호를 칩 코드(Chipping Code, Chip)라고 하는 다른 2진 코드(비트 패턴)로
변조토록하여 사용주파수 전역으로 확산시켜 전송한다.
- 데이터의 각 비트(Bit)를 칩(Chip)이라는 형태의 여러 비트들로 변조하여 사용주파수 전역으로 확산시켜
전송하고, 수신시 칩이 다시 원래의 비트 단위로 변환되어 데이터 복원
- 칩이 크면 클수록 원래의 데이터로 복원될 가능성이 커진다 (단, 더넓은 대역폭이 필요함)
ㅇ 기존 RF 신호는 신호가 특정 단일 주파수에 집중되어 신호세기가 크나, 확산 대역신호는 넓은 주파수대역에 퍼져 있어
신호세기가 매우 약하므로 기존 수신기로는 수신하기 어렵다.
- 저에너지의 광대역 잡음으로 인식하여 DSSS 신호를 무시함
ㅇ 칩 코드를 이용해 신호를 전대역폭을 이용해 확산 전송하므로 대역폭 이용효율은 떨어지나 향상된 잡음비와 칩 코드가 다르면
같은 주파수대역에서 여러 스프레드 스펙트럼 시스템이 사용 가능
ㅇ 변조방식은 BPSK, QPSK, M-ary PSK 등이 많이 사용된다.

3. 관련참조용어 : CDMA, 확산율(Spreading Ratio)

CDMA관련: 대역확산기술

2010-08-12T17:52:00.001+09:00

1. 개요
대역확산기술이란 여러 사용자가 같은 시간에 같은 주파수에서 동시에 정보를 다중으 로 전송하되 서로 직교성(Orthogonality)을 가진 코드를 사용하여 정보를 전달하는 기술을 말한다.

2. 특징
ㅇ 확산은 전송대역을 정보전송에 필요한 대역폭 보다 훨씬 넓게 사용한다.
ㅇ 송수신자간에 미리 약속된 확산코드에 의한 대역확산을 하므로 비화 특성이 우수
ㅇ 대역확산은 주파수대역을 넓게 사용하여 마치 주파수 다이버시티 효과를 얻을 수 있어서, 페이딩에 강하다.

3. 주요 대역확산 방식 분류
ㅇ 직접시퀀스 (Direct Sequence, DSSS)
- 고정된 주파수 대역을 사용하여 전송, 높은 속도 가능
- 例) IS-95 (이동전화 CDMA), 802.11b 등

ㅇ 주파수도약 (Frequency Hopping, FHSS)
- 여러 주파수 대역 사이를 호핑하며 전송, 다소 저속이나 구현비용이 낮음
- 例) 블루투스, HomeRF 등

한국의 핑거프린팅 기술연구동향: 핑거프린팅 논문/특허 목록

2010-08-12T15:08:00.001+09:00

BIBD기반의멀티미디어핑거프린팅코드의공모코드들에대한공모자추적.pdf 946K
JPEG압축및공모공격에강인한디지털이미지핑거프린팅기술.pdf 354K
LDPC와BIBD를이용한공모된멀티미디어핑거프린트의검출.pdf 467K
Yong_lee의익명핑거프린팅프로토콜의안전성취약점및개선방안.pdf 225K
가우시안의차를이용하여검색속도를향상한최소오디오핑거프린팅.pdf 1154K
공모공격에강인한디지털핑거프린팅기술.pdf 681K
공모공격에안전한불확정전송프로토콜기반의디지털핑거프린팅기법.pdf 465K
기하학적공격에강인한핑거프린팅기술개발.pdf 644K
디지털콘텐츠의지적재산권보호를위한익명핑거프린팅의연구동향.pdf 631K
디지털핑거프린팅에대한공모공격기술.pdf 619K
벨소리핑거프린팅_특허.pdf 365K
부정자추적후의재분배를고려한익명핑거프린팅.pdf 295K
비밀분산법을이용한익명성보장핑거프린팅기법.pdf 456K
삼차원메쉬모델에적용한공모방지핑거프린팅기법.pdf 814K
쉬프팅기법을이용한디지털이미지핑거프린팅기술.pdf 289K
시간축웨이블릿변환에의한비디오핑거프린팅.pdf 299K
시간축웨이블릿변환을이용한블라인드비디오핑거프린팅.pdf 602K
신경회로망에의한공모된멀티미디어핑거프린트의검출.pdf 527K
신뢰기관비참여의핑거프린팅기법에관한연구.pdf 497K
신뢰기관을이용한비대칭핑거프린팅기법.pdf 198K
신뢰기관을이용한익명적비대칭핑거프린팅기법.pdf 434K
압축 도메인 특징을 이용한 강인한 오디오 핑거프린팅.pdf 434K
오디오핑거프린팅기술과응용.pdf 611K
이미지콘텐츠출력물의저작권보호를위한디지털핑거프린팅기술에관한연구.pdf 282K
익명적비대칭핑거프린팅기법의보안취약성분석및개선방안.pdf ) 765K
재분배자의자동식별기능을갖는효율적인익명성을제공하는핑거프린팅.pdf 581K
콘텐츠불법사용자추적을위한디지털비디오핑거프린팅.pdf 524K
핑거프린팅기법을이용한불법콘텐츠추적시스템.pdf 356K
핑거프린팅을이용한디지털원문보호시스템설계.pdf 571K

1020040069959.pdf
1020050086557.pdf
1020060076670.pdf
1020060129447.pdf
1020067005805.pdf
1020070041754.pdf
1020070042047.pdf
1020077025412.pdf

--
연변대학교 경제관리학원 정보관리&정보시스템학과
延边大学经济管理学院信息管理与信息系统

CDMA의 원리: IS-95의 특성(참고: 노종선 교수 문서)

2010-08-12T13:38:00.001+09:00

IS-95의 특성

CDMA 방식 이동전화시스템의 표준으로 현재 미국 및 우리나라에서 IS-95라는 것이 사용되고 있는데 그 특징을 보면 다음과 같다. 순방항채널에서 사용되는 주파수대역은 869.4 MHz-893.9 MHz이고 역방향채널에서는 824.04 MHz-848.97 MHz이고 한 채널이 점유하는 주파수 대역폭은 1.23 MHz이다. 그리고 (음성) 데이터의 전송률은 1200, 2400, 4800, 9600 bps가 사용되고 있는데 이는 음성의 특성에 따라 가변으로 음성을 부호화하여 전송한다. 그리고 전송되는 데이터가 채널에서 생기는 비트오류를 정정하기 위해서 필요한 채널부호로서 길쌈부호기(convolutional encoder) 및 비터비 복호기(Viterbi decoder)가 사용되고 있다. 순방향채널에서는 부호율이 1/2이고 구속장이 9인 길쌈부호가 사용되고 있고 역방향채널에서는 채널이 열악한 상태에 있어 부호율이 1/3이고 구속장이 9인 성능이 우수한 길쌈부호가 채택되어 사용되고 있다. 그리고 변조방식으로는 변조효율이 우수한 4-ary 변조방식이 사용되는데 순방향채널에서는 QPSK 변조방식이 표준으로 채택되어 사용되고 있으며 역방향채널에서는 envelope fluctuation이 3dB로 비교적 우수한 OQPSK 변조방식이 사용되고 있다. 따라서 역방향채널의 송신기인 이동국에서는 전력증폭기로서 가격이 저렴하고 전력효율이 우수한 비선형 증폭기를 사용할 수 있다.

1. 순방향채널
우선 순방향채널은 파일롯채널(Pilot channel), 동기채널(Sync channel), 페이징채널(Paging channel), 그리고 트래픽채널(Traffic channel)로 구성되어 있다. 우선 파일롯채널은 기지국에서 항상 송신되는 비변조된 CDMA 채널신호로서 이동국에서 수신신호의 복조를 위한 초기동기를 얻기 위해 사용되는 채널이다. 기지국에서 송신되는 여러 가지 채널 중에 가장 큰 출력으로 전송되는 채널이라 할 수 있다. 기지국마다 서로 다른 위상의 파일롯 PN 수열을 사용함으로서 기지국의 구분이 가능하도록 한다.

동기채널은 기지국에서 1200 bps의 정보전송율로 항상 송신되는 채널로서 동기채널을 통하여 전송되는 정보는 해당 기지국을 구별하는 파일롯 PN 수열의 오프셋, 긴부호 생성기의 레지스터의 상태, 현재의 시스템 타임 등이며 이를 통하여 긴부호 동기와 시스템 타이밍을 획득할 수 있다. 동기채널은 80 ms구간에 3개의 프레임이 있다.

페이징채널은 기지국에서 항상 송신되며 기지국마다 최대 7개까지 송신될 수 있는 채널로서 이동국은 이 채널을 통하여 시스템 정보, 전체 또는 특정 이동국에 대한 명령이나 페이징정보, 이동국이 기지국에 액세스를 요구하였을 때 응답메시지 등의 정보를 포함하고 있다. 이동국이 페이징채널을 선택하는 방법은 이동국의 ESN과 미리 결정된 알고리즘에 의한다. 이 채널의 프레임은 20 ms로 기지국은 슬롯모드로 동작하는 이동국을 서비스하기 위해 200 ms 길이의 페이징채널 슬롯을 나누어 운영하는 선택기능도 있다.

트래픽채널은 호가 접속된 상태에서 기지국이 특정 이동국에 데이터 혹은 1200, 2400, 4800, 9600 bps로 부호화된 음성정보를 전송하는 채널로서 프레임은 20 ms로 구성되어 있다. 전송되는 정보의 전송율에 따라 전송되는 신호 전력의 세기를 조절하여 전송하게 된다. 즉 9600 bps 음성데이터의 경우는 4800 bps 음성데이터의 경우보다 2배의 출력으로 전송하게되고 1200 bps 음성데이터보다는 8배의 출력전력으로 전송하게 된다. 따라서 낮은 정보전송율의 경우에는 간섭 및 전력소모를 줄일 수 있게 된다. 앞서 설명한 바와 같이 이러한 방식에 의해 간섭파 잡음을 감소시켜 시스템의 용량을 증대시킬 수 있게 된다.

2. 역방향채널

역방향채널은 액세스채널(Access channel)과 트래픽채널(Traffic channel)로 구성되어 있다. 액세스채널은 이동국이 기지국에 접속을 시도할 때 사용하는 채널로서 정보전송율이 4800 bps이다. 사용자가 기지국에 접속을 요구할 때 이동국은 액세스채널을 통하여 기지국에 호접속을 시도한 후 기지국으로부터 페이징채널로 전송되는 모든 메시지를 수신하여 해석한다. 각 액세스채널은 슬롯단위로 액세스 시도를 하는데 슬롯의 선택은 랜덤과정에 의해서 결정되기 때문에 다수의 이동국이 호시도 중에 동일한 채널과 슬롯으로 전송하는 이동국들이 발생할 수 있다.

트래픽채널은 순방향의 경우와 마찬가지로 데이터 또는 1200, 2400, 4800, 9600 bps의 부호화된 음성데이터를 전송하는 채널로서 프레임의 크기는 20 ms이다. 순방향채널과는 달리 역방향 트래픽채널에서는 다른 정보전송율의 경우 출력신호의 전력을 변화시키는 것이 아니라 출력전력은 일정하게 유지하면서 duty cycle을 조절하여 음성데이터를 전송한다. 즉 9600 bps 음성데이터의 경우 100 %, 4800 bps 음성데이터의 경우 50 %, 2400 bps 음성데이터의 경우 25 %, 그리고 1200 bps 음성데이터의 경우 12.5 %의 duty cycle로 전송하여 신호의 간섭파 잡음을 감소시키게 된다.

결 론
확산 스펙트럼 통신방식 중의 하나인 CDMA방식의 이동전화시스템은 그의 구현에 있어서 상당한 어려움이 있음에도 불구하고 대단히 우수한 기술임에는 틀림이 없다. 특히 세계 최초로 우리나라에서 대규모의 상용화에 성공한 시스템으로서 많은 장점을 가지고 있다. 특히 CDMA방식의 이동전화시스템은 기존의 아날로그 이동전화통신방식에 비해 상당히 용량을 증가시킬 수 있고 또한 주파수계획의 필요성이 없으며, 음성 활성도를 쉽게 활용할 수 있고, 도시지역에서 생기는 다중경로파에 의한 페이딩을 쉽게 극복할 수 있다는 장점이 있다. 더욱이 CDMA 방식은 향후 IMT-2000시스템 및 미래의 많은 다른 무선통신시스템에 활용될 수 있는 핵심기술요소로서 이에 관련된 기술개발 및 기반기술 구축은 향후 정보통신 산업발전에 중요한 역할을 할 것으로 사료된다.

참 고 문 헌
[1] Joint Technical Committee of Committee T1 R1P1.4 and TIA 46.3.3/TR45.4.4 on Wireless Access, "RF Channel Characterization and System Deployment Modeling," Paper No. JTC(AIR)/94.09.23-065R6, Sep. 23 (1994).
[2] TIA/EIA, INTERIM STANDARD, Mobile Station -Base Station Compatibility Standard for Dual -Mode Wideband Spread Spectrum Cellular System. IS-95, Jul. (1993).
[3] S. W. Golomb, Shift Register Sequences. (San Francisco, CA, Holden-Day, 1967); revised edition (Laguna Hills, CA, Aegean Park Press, 1982).
[4] R. A. Scholtz, IEEE Trans Commun. COM-30, 822-854, (1982).
[5]R. A. Scholtz, IEEE Trans Commun. COM-25, 748-755 (1977).
[6]A. J. Viterbi, IEEE Journ. Selected Areas in Comm. 8, 641-649 (1990).
[7]K. S. Gilhousen, et al., IEEE Trans. Veh. Tech. 40, 303-312 (1991).

노종선 교수는 미국 남가주대에서 공학박사학위(디지털통신 및 정보이론분야)를 취득한 후 미국 Hughes Network Systems에서 책임연구원으로 근무하고 건국대학교 부교수를 거쳐 현재 서울대학교 전기공학부에서 조교수로 재직 중이다. 또한 서울대학교 부설 뉴미디어통신공동연구소 차세대무선통신연구센터 센터장을 맡고 있다.( jsno@snu.ac.kr)

--
연변대학교 경제관리학원 정보관리&정보시스템학과
延边大学经济管理学院信息管理与信息系统

CDMA 방식의 특성(전력제어, 소프트핸드오프, 시스템용량의 증가)

2010-08-12T13:35:00.001+09:00

CDMA 방식의 이동전화시스템은 다른 방식에 비해 여러 가지 다른 특성을 가지고 있다. 즉, 전력제어, 소프트핸드오프 (soft handoff), 시스템용량의 증가 등이 있다.

1. 전력제어

DS-CDMA 방식의 이동통신시스템에서는 여러 사용자가 같은 채널에서 시간과 주파수대역을 공유하며 신호를 전송하기 때문에 다른 사용자의 신호가 나에게는 간섭파 잡음이 된다. 따라서 어느 특정한 사용자의 수신신호가 다른 사용자의 수신신호보다 큰 경우에는 간섭파 잡음이 커지기 때문에 모든 사용자의 수신신호의 전력이 수신측에서 볼 때 모두 같아야 된다는 조건이 필요하다. 특히 CDMA 방식 이동전화시스템의 역방향채널의 경우 이동국이 기지국 근처에 있을 때와 셀의 경계에 있을 경우에 기지국에서 두 이동국으로부터 수신되는 신호의 전력의 비는 통산 100[dB] 이상 되는 경우가 흔히 존재한다. 이러한 문제를 DS-CDMA 다원접속방식에서 near-far 문제라 한다. 따라서 이러한 문제를 해결하기 위해서 이동국에서는 기지국 근처에 있을 경우와 셀의 경계에 있을 경우에 출력신호의 전력을 조절해야 한다. 이를 전력제어(power control)라 한다. 실제 사용되는 이동국에서는 셀 경계에 있을 경우에는 기지국 근처에 있을 경우보다 90[dB] 이상 큰 출력전력으로 신호를 전송한다. 따라서 실제 이동전화를 사용하는 경우에 기지국으로부터의 수신신호 레벨이 작은 경우에는 이동국이 기지국으로부터 멀리 있다는 것으로 내가 사용하는 이동국에서는 매우 큰 출력전력으로 신호가 전송되기 때문에 전자파에 의해 건강에 미치는 영향은 훨씬 크다고 할 수 있을 것이다. 그리고 이동국에서의 90[dB]이상의 전력제어를 실현한다는 것은 매우 어려운 문제이며 만일 이동국의 전력제어회로의 오동작이 생기면 그 셀 또는 섹터 내에 있는 모든 이동전화는 모두 동작불능 상태가 될 수 있으므로 역방향채널에서의 전력제어는 매우 중요한 문제라 할 수 있다. 또한 이동국은 계속해서 움직이고 있기 때문에 기지국에서 특정한 이동국으로부터 수신되는 신호의 전력은 변하게 되므로 이동국에서의 전력제어는 빠르게 이루어져야 한다. 현재 CDMA 방식의 이동전화 시스템에서 사용되는 전력제어 방식은 개방루프 전력제어(open loop power control) 방식과 폐쇄루프 전력제어(closed loop power control) 방식이 동시에 사용되고 있다.

개방루프 전력제어방식은 이동국에서 기지국으로부터 수신되는 신호의 전력 세기를 측정하여 이동국 출력신호의 전력을 조정하는 방식으로서 전력제어가 빠르게 이루어지는 장점이 있다. 이 방식은 순방향채널과 역방향채널이 어느 정도 동일하다는 가정 하에 이루어지는 전력제어라 할 수 있다. 반면에 폐쇄루프 전력제어방식은 기지국에서 이동국으로부터 수신되는 신호의 세기를 측정하고 이동국에 전력제어 메시지를 전송하면 이동국에서는 그 메시지 따라 출력 신호의 전력을 조정하는 방식으로 전력제어는 느리게 이루어지지만 보다 세밀한 전력제어를 할 수 있다는 장점이 있다.

2. 소프트핸드오프

이동통신시스템에서 핸드오프(handoff)란 이동국이 현재 셀영역을 벗어나 새로운 셀로 진입하여 새로운 기지국을 이용하여 통화를 계속하거나 동일 셀 내에서 무선전파환경의 악화로 새로운 채널을 할당받아 통화의 연속성을 유지하는 절차를 말한다. 핸드오프는 크게 하드핸드오프(hard handoff)와 소프트핸드오프(soft handoff) 그리고 소프터핸드오프(softer handoff)로 나뉠 수 있다.

소프트핸드오프 및 소프터핸드오프는 CDMA 방식의 이동전화시스템에서 수신신호의 복조의 어려움 때문에 처음으로 도입된 핸드오프 방식으로 여러 가지 장점과 단점을 가지고 있다. 소프트핸드오프는 한 기지국과 통화가 끊어지기 전에 다른 기지국과 통화가 연결되는 방식(make-before-break operation)으로서 셀의 경계지역에서 핸드오프 상태에 있는 단말기는 두 개 이상의 기지국과 통화를 하는 것으로 이를 하드핸드오프를 하는 경우와 비교하면 기지국간의 디버시티에 의해 통화성능이 향상될 수 있는 반면에 핸드오버를 위해 채널카드와 같은 추가적인 장비와 순방향에서의 송신전력 증가와 같은 단점을 가진다. CDMA 방식의 이동전화시스템에서는 위의 세가지 핸드오프 방식을 포함하여 다음과 같은 핸드오프 방식이 존재하는데 그 특징 등을 살펴보면 다음과 같다.

●소프터핸드오프 : 이동국이 통화의 단절없이 동일한 기지국내에서 기존의 섹터에서 새로운 섹터로 통화를 인계하여 통화의 연속성을 유지하는 방식으로 동일한 셀 내에서만 이루어진다.
●소프트핸드오프 : 이동국이 통화의 단절 없이 새로운 기지국과 통신을 시작하는 방식으로 동일한 주파수가 할당된 CDMA 채널 사이에서만 이루어진다.
●CDMA-CDMA 하드핸드오프 : 다른 주파수 또는 다른 프레임 옵셋을 갖는 기지국과 서로 다른 기지국군 사이를 움직일 때 발생한다.
●CDMA-아날로그 핸드오프 : 이동국이 순방향 통화채널에서 아날로그 음성채널로 전환할 때 발생한다.

그리고 소프트핸드오프 상태에 있는 이동국이 두 개의 기지국 또는 섹터와 동시에 통화를 하는 경우 이를 2-way 핸드오프 그리고 세 개의 기지국 또는 섹터와 동시에 통화를 하는 경우를 3-way 핸드오프라 한다.

3. CDMA 방식의 시스템 용량
CDMA 방식이 미국에서 표준으로 채택될 수 있었던 이유 중의 하나는 시스템의 용량이 매우 크다는 것이다. 따라서 CDMA 방식의 이동전화시스템에서 시스템의 용량의 계산은 중요한 문제라 할 수 있다. 일반적으로 서로 다른 방식의 이동전화시스템에 대해서 각각의 시스템의 용량을 비교한다는 것은 여러 가지의 요소들이 고려되어야 하기 때문에 대단히 어려운 문제이다. 근본적으로 CDMA 방식 이동전화시스템에서 시스템 용량의 증대는 다른 사용자의 간섭파 잡음을 어떻게 감소시키는가에 있다고 할 수 있다. CDMA 방식에서 시스템의 용량을 증대시키는 여러 가지 요인들에 대해서 살펴보도록 하겠다.

우선 사용자가 통화중에 음성신호를 전송하는 경우 음성신호의 묵음 구간이 자주 존재하고 또한 상대방의 대화를 청취하는 경우에는 음성신호의 묵음구간이 매우 길어지는 특성이 있다. 이러한 묵음구간에 CDMA 방식에서는 신호를 전송하지 않으므로 전체적인 사용자들의 간섭신호는 감소하게 된다. 통계적으로 보면 통화중에 약 40 % 정도의 시간구간 동안에만 음성신호를 전송하게 되어 간섭파 잡음이 약 40 % 정도로 줄게 되므로 시스템용량의 증대효과가 2.5배에 이르게 된다. 이를 voice activity detection에 의한 시스템의 용량 증대라 한다. 또한 현재 사용되고 있는 각각의 셀은 보통 세개의 섹터로 분할되어 사용되고 있는데 이러한 섹터의 분할에 의해 채널 상에 존재하는 간섭파 잡음은 1/3로 감소하게 된다. 따라서 시스템의 용량은 3배로 증가하게 된다. 또한 CDMA 방식에서의 주파수재사용(frequency reuse) 효율은 1인데 이는 모든 셀 및 섹터에서 같은 주파수를 사용할 수 있다는 것으로 이에 의한 시스템의 용량 증대 효과는 매우 크다고 할 수 있다. 이러한 여러 가지 요인으로 인하여 CDMA 방식의 시스템용량은 주어진 주파수 스펙트럼에서 다른 방식들에 비해 매우 크다고 할 수 있다.

그리고 TDMA 및 FDMA 방식에서는 현재 사용 가능한 채널이 없으면 사용자의 새로운 통화연결이 불가능하지만 CDMA 방식에 있어서는 시스템 용량에 융통성이 있어 간섭파 잡음을 약간 증가시키면서 항상 사용자의 새로운 통화연결이 가능하게 될 수 있다는 특징이 있다. 이를 소프트용량(soft capacity)라 한다.

4. 확산부호

DS-CDMA 방식의 통신시스템에서 확산부호로 사용될 수 있는 것으로 PN 수열이 있다. 이는 확산대역통신방식의 핵심적인 요소기술로서 확산대역통신방식의 발전과 더불어 발전되어 왔다. 그런데 최근 들어서 PN 수열은 그 응용분야가 점점 넓어지면서 그 중요성이 증대되고 있다. 즉 PN 수열은 확산대역통신방식의 직접수열방식 및 CDMA 통신방식에 사용될 뿐 아니라 데이터나 영상신호 및 음성신호의 스크램블링을 위해서 사용되고 있다. 그리고 PN 수열은 컴퓨터에서 랜덤 데이터를 발생시키는데 사용되며 암호화과정에서 키스트림(key stream)으로 사용되고 있다. 이러한 PN 수열은 우선 다음과 같은 경우가 대부분 활용되어 왔다.

●0 또는 1의 값을 갖는 이진 PN 수열
●주기는 N=2n -1, 여기서 n은 정수, 즉
N=7, 15, 31, 63, 127, 255, 511, 1023, 2047, 4095, ...

●주기적인 PN 수열 즉,
s(t)=s(t +N), 0≤t≤N-1
그리고 PN 수열들이 가지고 있어야 할 바람직한 성질들을 보면 아래와 같다.

●낮은 위상어긋나기(out-of-phase) 자기상관 특성
●낮은 상호상관 특성
●큰 선형구간(linear span)
●심볼의 균형
●많은 서로 다른 수열군들의 존재

이러한 좋은 특성을 갖는 수열 중에 m-수열라는 것이 있는데 이들은 IS-95 CDMA 방식의 이동전화시스템에서는 주기가 32768인 Pilot PN 수열라는 이름으로 사용되고 있고 또한 주기가 242-1인 수열이 긴부호(long code)라는 이름으로 매우 중요한 역할을 하며 사용되고 있다.

또한 CDMA에서 신호를 확산시키기 위하여 사용되는 것으로 Walsh 부호 또는 하다마드 행렬(Hadamard matrix)이라는 것이 있다. 하다마드 행렬은 각 요소들이 이진인 +1, -1로 구성되는 행렬로서 IS-95의 순방향채널에서는 신호의 대역 확산 및 사용자 구분을 위해서 사용되고 있다. 하다마드 행렬은 여러 가지 방식으로 생성될 수 있는데 IS-95에서 사용되는 하다마드 행렬은 실베스터 방식으로 생성되는 64 ×64 크기의 하다마드 행렬이 사용되고 있다. 그러나 역방향채널에서는 하다마드 행렬이 신호의 변조를 위해서 사용되고 있다.

5. 경로디버시티

CDMA 방식의 이동전화시스템은 확산부호를 이용하여 전송신호를 확산시킴으로서 매우 넓은 주파수스펙트럼을 갖는 신호를 전송하게 된다. 따라서 CDMA 방식은 다중경로 페이딩을 극복할 수 있는데 이는 도심지와 같이 다중경로 페이딩이 심한 지역에서 매우 우수한 특성으로 간주될 수 있다. IS-95에서 사용되는 확산부호의 칩주기(chip period)는 약 0.8138 μsec로서 시간지연이 이보다 더 크게 되어 수신되는 다중경로파는 모두 RAKE 수신기에 의해 분리되어 복조됨으로서 페이딩 현상을 극복할 수 있다. RAKE 수신기는 여러 개의 시간지연을 갖는 복조기로 구성되어 시간지연이 칩주기보다 큰 다중경로파는 모두 복조할 수 있게 된다. IS-95의 경우 기지국에서는 수신기가 4개인 RAKE 수신기가 사용되고 있고 이동국에서는 수신기가 3개인 RAKE 수신기가 사용되고 있다. 각 복조기에서 복조된 신호는 결합하여 보다 고품질의 데이터를 얻을 수 있게 된다.

연변대학교 경제관리학원 정보관리&정보시스템학과
延边大学经济管理学院信息管理与信息系统

CDMA 다원접속 방식의 원리 및 종류

2010-08-12T13:29:00.001+09:00

확산대역통신방식은 송신기에서 정보신호의 스펙트럼을 넓게 확산(spreading)하여 전송하고 수신기에서는 수신된 신호의 주파수 스펙트럼을 역확산(despreading)하여 복원하는 방식이다.

확산대역통신방식은 정보신호의 주파수 스펙트럼을 확산하는 방식에 따라 직접수열(direct sequence) 방식과 주파수도약(frequency hopping) 방식이 있는데 이중 이동전화에서 사용되는 방식은 직접수열 방식의 부호분할다원접속(DS-CDMA) 방식이다. DS-CDMA 방식을 간단히 CDMA라고도 하는데 이 방식은 여러 사용자가 하나의 무선채널을 사용하는 경우 시간과 주파수를 공유하면서 사용자간의 간섭을 최소화하기 위하여 각 사용자에게 확산부호인 상호상관특성이 우수한 PN(pseudo-noise) 수열을 할당하고 각 사용자는 할당된 PN 수열을 이용하여 송신할 신호를 확산하여 전송하고 수신 측에서는 송신 측에서 사용한 것과 같은 PN 수열을 발생시켜서 동기를 맞추고 이를 이용하여 수신된 신호를 역확산하여 원하는 신호를 복원한다. 그림 2에서 CDMA 다원접속방식의 블럭도를 도시하였다.

CDMA 방식은 또한 다음과 같이 두 가지로 분류될 수 있다. 여러 사용자가 신호를 전송할 때 반송파 및 PN 수열의 동기를 맞추어 여러 사용자의 신호들을 전송하게 되는 경우를 동기식(synchronous) CDMA라 하고, 그렇지 않은 경우를 비동기식(asynchronous) CDMA라 한다. 이동전화시스템의 순방향채널(forward channel; 기지국에서 이동국)은 동기식 CDMA 방식이 사용되고 역방향채널(reverse channel; 기지국에서 이동국)은 비동기 방식이 사용되고 있다.

이러한 CDMA 다원접속방식의 이동전화시스템은 다음과 같은 장점을 있다.
·다른 방식에 비하여 채널의 용량을 크게 증가시킬 수 있 다.
·음성신호를 전송할 경우 음성 활성도(voice activity)를 이 용하여 채널의 사용효율을 높일 수 있다.
·시스템의 용량이 다른 방식에 비해 융통성이 있다. 즉 사 용자가 늘어나면 채널의 잡음이 증가하여 상태가 나빠지 게 되지만 시스템의 용량은 증가시킬 수 있다.
·전송신호가 PN 수열을 사용하여 확산되므로 사용자의 신호에 대한 개인비밀이 보장된다.
·다중경로 페이딩에 강한 특성을 갖는다.
·협대역 잡음신호에 강인한 특성을 갖는다.
·사용자간의 채널할당이 필요하지 않다.
·주파수설계가 필요하지 않다.

--
연변대학교 경제관리학원 정보관리&정보시스템학과
延边大学经济管理学院信息管理与信息系统

CDMA의 원리

2010-08-12T13:26:00.001+09:00

이동전화는 1980년 초반부터 아날로그통신방식을 이용하여 그 서비스가 시작되었는데 할당된 사용 가능한 주파수 스펙트럼은 한정되어 있고 가입자수는 계속적으로 증가하여 1990년 초에 이동전화의 시스템 용량 (system capacity) 증대가 매우 중요한 문제로 대두되었다. 이를 해결하기 위하여 1990년대 중반부터 전세계적으로 아날로그방식의 이동전화에서 디지털방식의 이동전화시스템으로의 진화가 이루어지게 되었다. 디지털방식의 이동전화시스템은 근본적으로 송신되는 메시지가 디지털 데이타라는 점이다. 이러한 디지털 데이터를 전송하기 위하여 디지털통신방식이 사용되어야 하는데 이 방식에는 여러 가지 요소기술들이 있다. 예를 들면 아날로그 음성을 디지털 데이터로 바꾸어 압축하는 음성부호화 (speech coding), 채널에서 전송되는 데이터에 발생되는 비트오류를 정정하기 위한 오류정정부호, 신호의 전송 및 수신을 위한 변조기 및 복조기, 그리고 가장 중요한 요소기술인 채널에서의 다원접속(multiple access)기술 등이 있다.

다원접속기술은 하나의 무선통신채널에 여러 사용자의 신호를 사용자 서로간에 간섭을 일으키지 않고 전송할 수 있도록 하는 기술로서 여러 가지 디지털통신시스템의 요소기술 중에서 가장 중요한 것이며 또한 무선통신시스템의 구조를 가장 크게 바꾸어 주는 것이다. 따라서 디지털이동통신시스템의 방식을 나눌 때 다원접속방식에 따라 주파수분할다원접속 (FDMA: frequency division multiple access), 시분할다원접속 (TDMA: time division multiple access), 부호분할다원접속 (CDMA: code division multiple access) 등으로 나뉘게 된다.

CDMA 방식은 디지털이동전화시스템이 도입될 당시에 이동전화시스템의 가장 중요한 문제인 시스템의 용량을 가장 크게 증가시킬 수 있다는 장점이 있어 미국 샌디에고에 있는 Qualcomm사에서 제안하고 미국의 CTIA (Cellular Telecommunication Industry Association)에서 여러 이동전화 표준 중의 하나로 채택된 IS-95라는 시스템이다.

CDMA 통신방식은 군통신 시스템에서 주로 사용되어 왔던 확산대역통신방식(spread spectrum communication system)의 하나로 무선통신이론 중에서 가장 최신의 이론으로서 매우 구현이 어려운 방식이라 할 수 있으며 우리 나라에서 세계 최초로 대규모의 상용화 성공한 시스템이라 할 수 있다. CDMA 방식 단말기의 전체적인 블럭도는 그림 1과 같다.

이 글에서는 우선 CDMA방식의 원리 및 여러 가지 종류의 CDMA 방식 그리고 확산대역통신시스템에 대하여 논하였고 제3장에서는 CDMA 방식의 특성으로서 전력제어, 소프트핸드오프, 시스템용량, PN 수열(sequence) 및 경로디버시티 등에 관하여 기술하였으며 제4장에서는 현재 사용 중인 CDMA 이동전화시스템의 표준인 IS-95의 순방향채널의 종류 및 특성 그리고 역방향채널의 종류 및 특성에 관하여 기술하였고 제 V장에서 결론으로 끝을 맺었다.

참고: http://jncis.egloos.com/2332863

--
연변대학교 경제관리학원 정보관리&정보시스템학과
延边大学经济管理学院信息管理与信息系统

3D BARCODE(3차원 바코드)도 존재함

2010-08-11T18:48:00.001+09:00

1차원, 2차원바코드에 이어 3차원 바코드도 사용되고 있습니다. 국내에서는 아직 적용된 사례가 거의 없으나, 개발국가인 미국에서는 이미 사용된지 몇년 되었습니다.

1, 2차원 바코드가 주로 면에 바코드를 인쇄하여 사용되는 것이라면, 3차원 바코드는 면에 바코드를 인쇄하기 어려운 철강제품이나 금형주조물, 프라스틱 금형물 등에 요철을 주어 1,2차원의 해독방법인 바의 두께에 바코드의 깊이까지 해독에 사용됩니다.

--
연변대학교 경제관리학원
정보관리&정보시스템학과 최기철

延边大学经济管理学院
信息管理与信息系统崔基哲

数字指纹的研究发展方向

2010-08-11T18:44:00.001+09:00

数字指纹技术的研究领域中，在取得成绩的同时，我们仍要意识到目前还有很多问题需要解决，本文的工作相对整个数字指纹技术领域而言，仅是冰山一角，为丰富完善该领域，大量的工作等待我们继续研究探讨，仍需要进一步的努力。文章结束之时，我们依据数字指纹的系统模型，特提出一下几个值得研究的方向：

(1)数字指纹技术理论基础的研究。数字指纹技术属于信息隐藏技术中的一种。然而，信息隐藏技术发展到今天，还没有找到自己的理论依据，没有形成理论体系，很多信息理论的经典理论在处理信息隐藏系统时不再适用。随着数字指纹技术的不断发展，它对理论指导的期待已经越来越迫切。

(2)编码方案的研究。数字指纹系统的一个特点就是运算量大，为此在保证一定的抗合谋性能前提下，设计合理的编码方案，在码字长度和抗合谋性能之间寻求合理的折中方案需要人们的不断地努力。

(3)嵌入、提取和跟踪技术的同步完善。数字指纹系统需同时满足稳健性攻击和合谋攻击。然而，目前的大多数算法中，仅仅考虑系统的抗合谋性能，而忽略了系统的稳健性，有的算法仅仅就编码设计给出性能上分析，而却没有具体的实验验证，因此在指纹系统设计中，应充分考虑三者的联系。

(4)另外实验验证算法过程中，大多采用静止图像为研究对象，而视频、音频等又有其自身的特点，比如冗余性、感知性等，因此在以后的算法设计中应将媒体产品自身的特点考虑进去。

(5)系统整体设计研究。数字指纹技术在实用化过程中所需的指纹检测体系的建立、技术衡量标准、法律的保护等问题有待解决。

(6)结合数字技术、密码学、数字签名技术和访问控制技术，将数字指纹技术应用实际内容的保护，同时在应用中使其得以发展成熟。

我们相信，随着研究的不断深入，数字指纹技术将有良好的应用前景。

参考：合肥工业大学，硕士毕业论文，抗合谋数字指纹技术研究，姓名:李敏 20090301

指导教师：朗文辉

--
연변대학교 경제관리학원 정보관리&정보시스템학과
延边大学经济管理学院信息管理与信息系统

数字指纹的跟踪：三种不同的标准实施跟踪

2010-08-11T18:28:00.000+09:00

根据实际应用的需要，可采用以下三种不同的标准实施跟踪。

Catch one：在这种设计场合，目标是要设计指纹以最大化捕捉至少一个合谋者的机会，同时试图最小化错误地指控无罪用户的可能性。针对这个期望的目标，评价准则集由误报概率和漏报概率组成。从探测器的角度看，当探测器不能识别任何合谋者(漏报)或探测器错误表明无罪用户是一个合谋者(误报)时，探测方案就失败了。当在法庭上提供证据时，这个标准是有重大意义的。

Catch many：在这种设计场合，目标是要捕捉尽可能多的合谋者，尽管可能的成本是指控更多的无罪用户。针对这个期望的目标，评价准则集由期望的被成功捕捉的合谋者部分和期望的受到错误怀疑的无罪用户部分组成。

Catch all：在这种设计场合，设计指纹以最大化捕捉所有合谋者的概率，同时保持受错误指控的无罪用户的数量是可接受的。当信息收件人的信用受到极大关注，以至于涉及信息泄露的所有用户都需要识别时，就会出现这种情况。这个评价准则集由对捕捉所有合谋者的概率和一个实用率的度量组成，后者描述了每合谋者无罪的人受指控的期望值。

数字指纹合谋攻击：非线性合谋攻击

2010-08-11T18:24:00.001+09:00

一般情况下，除了平均攻击，合谋用户还经常使用非线性攻击方法，只要采用合适的非线性攻击方法，合谋者就可以改变多媒体拷贝中的任意部分的输出值，使输出值在允许的最大值和最小值之间，而且不影响拷贝的质量。常见的非线性攻击方法就是采用获取多媒体版本中数字指纹的最大值、最小值或者中间值。

在最小值攻击、最大值攻击、中值攻击下，合谋者产生的新拷贝的每个信号分别是所有合谋者当前位置信号的最小值、最大值和中值，并用此值取代正版版本中的数字指纹以获得最后的盗版拷贝。

除此之外，对视频序列进行旋转、比例压缩、剪切(RST)变换，其效果分析类似于前述的非线性攻击。

--
연변대학교 경제관리학원 정보관리&정보시스템학과
延边大学经济管理学院信息管理与信息系统

数字指纹合谋攻击：线性合谋攻击

2010-08-11T18:20:00.001+09:00

线性合谋是现在最简单易行的合谋攻击方法，例如，当K个用户试图进行合谋攻击时，我们只需要将他们的K个拷贝进行线性迭加然后除以K即可得到盗版拷贝。由于在这种攻击方式下，每名攻击者平均地承担了被抓获的风险，故而被经常用来产生合谋版本。

1)平均攻击：新拷贝的每个信号是合谋者当前位置信号的平均

2)进行加权平均后，再加上一个小功率的噪声。这个噪声可以进一步扰乱篡改拷贝中残留的数字指纹信息。

3) 采用多重输入单输出的线性LSI滤波，对添加的相互独立、且具有相同概率分布的数字指纹进行平均攻击。

值得注意的是，这种加权平均的攻击方法会削弱数字指纹的效果，而且人数越多，每个合谋者被追踪到的可能性就越小。通过这种平均加权攻击得到的多媒体拷贝往往比正版的版本在视觉、听觉效果上更出色。

--
연변대학교 경제관리학원 정보관리&정보시스템학과
延边大学经济管理学院信息管理与信息系统

数字指纹的攻击类型：合谋攻击

2010-08-11T18:12:00.001+09:00

合谋攻击
由于同一数字作品的每份拷贝都不相同，带来的问题是，多个用户可以联合他们的拷贝，比较拷贝中的不同之处，从而可以找出拷贝中的部分指纹标记的位置，对这些位置上的标记进行修改，选择其中的一种甚至创造新的标记，从而制造出一份新的作品拷贝，将其分发，希望能够逃避跟踪。这一过程就称为用户合谋，是盗版用户针对数字指纹系统的主要攻击方式。

合谋攻击仅针对同一作品不同拷贝中不一致的内容进行修改，是专门针对指纹标记的攻击方式。对于这类攻击，数字指纹技术主要通过在数字指纹编码过程中使指纹具有抗合谋能力来对指纹加以保护。

用户合谋时，可以选择不同的合谋方式，但总的来说，合谋用户会采用某一被认为最优的策略，使他们能够成功制造出某个不属于任何用户，或者属于无罪用户的指纹，从而免于责罚。通常，合谋用户的最优策略是使他们被抓获的可能性最小，又或是成功陷害某个无罪用户的概率最大。

如果盗版用户可以伪造多个合谋指纹时，他们将选择最符合上述策略的指纹。用户合谋产生新的作品拷贝，这一行为等同于制造新的用户指纹，数学上，可以将合谋过程理解为由用户指纹通过合谋运算得到合谋指纹。具体的合谋运算可以不仅仅是简单意义上的代数运算，可以是逻辑数学中的逻辑运算，统计数学中的统计运算，也可以是综合这些运算的多次组合运算。

数字签名的攻击类型：单用户攻击

2010-08-11T18:06:00.001+09:00

单用户攻击是指单个用户通过对其拥有的数字产品进行各种操作，以期去掉数字产品中的指纹或者将指纹转换成其它用户的指纹，即通过对作品拷贝进行各种处理来削弱或者删除作品拷贝中指纹，这时数字指纹系统和数字水印系统所遭遇的攻击基本相同。
一般来说攻击方法可以分为三类：去除攻击(Removal Attack)、同步攻击(Svnchronization Attack)、协议攻击(Protocol
Attack)。
(1)去除攻击。

也称为稳健性攻击，其目的是从数字产品中去除掉指纹(水印)，并且不影响数字产品的正常使用。这类攻击大体可分为两类：信号处理攻击和分析攻击。信号处理攻击包括一些基本的信号处理操作，如去噪、有损压缩、量化、重新调制、滤波、颜色校正等；分析攻击也称直接攻击，包括在指纹(水印)的嵌入和检测提取阶段采用特殊方法来擦除或减弱产品中的指纹(水印)。
(2)同步攻击。

也称为表达攻击，它不需要利用算法来去除数字产品中指纹(水印)，而是试图通过对掩蔽载体做各种修改以使得检测器失效，破坏了载体与指纹(水印)的同步性。这类攻击主要包括几何失真攻击(Geometrical Distortion Attack)、马赛克攻击(Mosaic Attack)、抖动攻击(Jitter Attack)、Oracle攻击。其中，几何失真攻击(几何攻击)由于简单易行，是较为常见的一类攻击，主要包括旋转(Rotation)、缩放(Scaling)、平移(Translation)、剪切(cropping)、图像反转等。其中，旋转、缩放、平移简称为RST。

(3)协议攻击。

解释攻击又称为死锁攻击，它通过伪造假的原始数字产品或掩蔽载体来制造混乱，从而任何人可以声明其对产品的所有权，导
致版权争议。拷贝攻击是从掩蔽载体中估计出指纹或水印，并将估计的内容加入其它数字产品，破坏版权产品的合法性。

对于这类攻击，数字指纹技术主要通过在数字指纹嵌入过程中提高嵌入指纹的稳健性来对指纹加以保护。
--
연변대학교 경제관리학원 정보관리&정보시스템학과
延边大学经济管理学院信息管理与信息系统

핑거프린팅의 분류방식

2010-08-11T17:57:00.001+09:00

1．基于客体的分类

客体的自然属性是一个最基本的标准，这是因为它能提供一种定制的方法为客体嵌入指纹。基于客体的分类有两种：数字指纹(Digital Fingerprinting)和物理指纹(Physical Fingerprinting)

2．基于检测灵敏度的分类

基于对侵害的检测灵敏度可把指纹分为：完美指纹(Perfect fingerprinting)、统计指纹(Statistical Fingerprinting)[1 6]和门限指纹(Threshold Fingerprinting)。如果对客体的任何修改使指纹不可识别的同时，也导致了客体不可用，则称这种指纹为完美指纹。

3．基于嵌入指纹方法的分类

基本的指纹处理方法，如识别、删除、添加、修改等已经被作为另一种分类标准。如果指纹方案由识别和记录那些已成为客体一部分的指纹组成，则称之识别类型(Recognition Type)，如人的指纹和虹膜。删除类型(Deletion Type)指纹是指删除客体的一些合法部分作为客体的指纹。若在客体中加入一些新的成分来嵌入指纹，那么它就属于添加类型(Addition Type)。

添加的部分可以是敏感的，也可以是无意义的。若修改客体的某部分来嵌入指纹，则为修改类型(Modification Type)，例如变换颜色组成的地图等。

4．基于指纹属性的分类

如果生成的指纹是有限的离散取值，那么就称该指纹为离散指纹(Discrete Fingerprinting)，例如数字文件的哈希值。如果生成的指纹是无限的连续取值，那么就称该指纹为连续指纹(Continuous Fingerprinting)，大部分物理指纹属于这种类型。

--
연변대학교 경제관리학원 정보관리&정보시스템학과

延边大学经济管理学院信息管理与信息系统

数字指纹在版权保护中的应用

2010-08-11T17:34:00.001+09:00

1)对原创者的鉴别，将原创者的指纹(如创作者的信息及创作时间等)嵌入到所创作的作品中，这种指纹要求具有强鲁棒性。即使是放在Internet上，也可鉴别出作品的来源。

2)用作数据真实性鉴别，如记者用数码相机拍摄的图像，在使用这些图像前，新闻机构需要核实这些图像在拍摄后有没有进行修改或编辑。为了以后能够验证，需要在拍摄时嵌入指纹(由摄像机拍摄时自动将指纹嵌入所拍摄的图像中)。使用这些图像时，如果图像上存在指纹，表明该图像拍摄后未经修改。这种指纹要求具有低的鲁棒性，即脆弱性指纹。

3)未授权的拷贝，所有者(销售者)关心的是购买者是否会免费分发所购买的作品，这样就会降低销售者的版权收入。为了解决这一问题，销售者在销售自己的作品时，每份作品中都应嵌入自己的指纹(作品的编号和销售者的信息等)。当怀疑有未授权的拷贝时，从拷贝的作品中提取出指纹，通过对指纹的检测就可判断出是否有未授权的拷贝。这种指纹要求具有好的鲁棒性，并且能够快速提取。

4)所有权证明，在数字作品中嵌入销售者的指纹。如果有人在未付版费的情况下编辑和发行销售者的作品，销售者就可以从盗版作品中提取指纹，根据指纹来证明自己拥有该作品的版权。这种指纹除要求鲁棒性和可逆性外，还可能需要多重指纹方案。

5)侵权者的证明，除所有权证明之外，销售者还希望知道谁在非法分发自己的作品。为此，在销售时应嵌入用户的指纹(用户信息和交易时间等)。当销售者一旦发现盗版作品，就可从盗版作品中提取指纹，通过对所提取的指纹进行分析，就可证明谁是非法再分发者或合谋者，从而追究其责任。这种指纹要求具有抗合谋能力。

--
연변대학교 경제관리학원 정보관리&정보시스템학과
延边大学经济管理学院信息管理与信息系统

中国对数字指纹的研究动向

2010-08-11T17:27:00.001+09:00

1999年以来，Min wu小组在多媒体抗合谋数字指纹[8]方面一直有论文出现，且数字指纹从理论分析到具体仿真方面都取得了良好的成绩。最近，Celik等针对静止图像提出了一种先发制人的抗合谋指纹方案[9]。该技术在发行前对宿主信号的几何特征施加随机、不可感知的扭曲。Mao[10]等将这一思想推广应用于视频内容，他们采用视频帧随机重采技术，为每个用户发送一个有轻微差异的视频拷贝。国内研究力量主要集中于大学和研究所、研究院，其他如微软研究院等机构也开展了这方面的研究工作[11]

[8]W．Trappe， M． Wu， and K．J．R． Liu． Collusion・resistant fingerprinting for multimedia[C]．in Proc．IEEE Int．Con￡Acoustics，Speech，Signal Processing，2002，PP．3309-3312

[9]M U．Celik，G Sharma,A MTekalp．Collusion-Resilient Fingerprinting Using Random Pre—Warping[J】．IEEE Sig．Proc．Lett．，2004，11(1 O)：83 1・35．

[10]Y．Mao，M K Mihcak．Collusion—Resistant Intentional De-Synchronization for Digital Video Fingerprinting【J】．In Proceedings of IEEE International Conference on Image Processing(ICIP)，2005，237-240．

[11]吕述望，王彦，刘振华．数字指纹综述【J】．中国科学院研究生院学报．2004，21(3)：289．298．

广泛采用的数字水印追踪算法(TA算法)的复杂度

2010-08-11T17:23:00.001+09:00

另外，目前广泛采用的追踪算法(TA算法)的运行时间复杂度为0(，z)，即运行时间随着用户数的增加而增加。这样，整个数字指纹系统的执行效率将随着用户数的增加而大大的降低。
--
연변대학교 경제관리학원 정보관리&정보시스템학과
延边大学经济管理学院信息管理与信息系统

Pfitmann和Schunter: 一种匿名的指纹方案

2010-08-11T17:18:00.001+09:00

Pfitmann和Schunter还给出一种匿名的指纹方案，它类似于盲签字，该方案使用一个可信的称作注册中心的第三方来识别被怀疑有非法行为的买方，也就是卖方没有注册中心的帮助就无法识别买方。

B Pfitzmann，M Waidner．Anonymous fingerprinting[R]．In：Walter Fumy．Advances in Cryptology-EUROCRYPT 97．I，ectIlre Notes in Computer Science 1 233．Springer-Verlag，1997．88-102．

--
연변대학교 경제관리학원 정보관리&정보시스템학과
延边大学经济管理学院信息管理与信息系统

cuijizhe's DSP Using MATLAB

프랙탈 주식구조 분석코드

Matlab에서 corrcoef로 얻은 Array

콘벌루션: conv([1 1 1],[1 1 1]) 실행결과는 1 2 3 2 1에 관한 해석

콘벌루션: conv([1 1 1],[1 1 1]) 실행결과는 1 2 3 2 1에 관한 해석

corrcoef은 어떻게 구한 것인가?

Matlab에서 corrcoef로 얻은 Array

PN Code, PN Sequence 만드는 함수

논문데이터 정리: 삽입알고리즘 절차

데이터분석: 발표할 논문데이터 정리(2010년11월)

최종욱: 워터마크 기술을 미래의 저작권 보호기술로 인식

대만 Matlab신호처리 매니아 사이트

툴박스 다운로드: 오디오신호처리관련 Utility Toolbox(1)

툴박스 다운로드: Speech and Audio Processing Toolbox(1)

이미지 워터마킹(삽입: 2-D wavelet transform)

wavread를 이용하여 오디오를 원하는 길이로 자르기

Matlab으로 투명한 실행창 만들기(가능할지?)

신호 대 잡읍 비 구하기 코드: SNR(오디오용)

Noise_CancellationUsing_LMS_algorithm

FIR filter coefficient estimation

AR_process_estimation(code)

신호 대 잡읍 비 구하기 코드: PSNR(이미지용)

Matlab에서 마우스로 좌표값을 얻는 방법

간단한 애니메이션 동영상(출처: http://pinkwink.kr/192)

영화 만드는 방법(사례)

Matlab 함수: movie 사용사례

Arnold反置乱(解密)程序之二

Arnold置乱(加密)程序(之一)

图像置乱度计算MATLAB代码（转存）

Matlab图像混沌置乱(转存)-采用的是混沌置乱算法

Matlab toolbox: 기하도형 그리기

Test Image(원과 직선)

ugh 변환_허프만방법으로 원, 직선을 검측하기(07.the_max1 )

Hough 변환_허프만방법으로 원, 직선을 검측하기(06.pll.m)

Hough 변환_허프만방법으로 원, 직선을 검측하기(05.machZ_try.m)

Hough 변환_허프만방법으로 원, 직선을 검측하기(04.gather_point.m)

Hough 변환_허프만방법으로 원, 직선을 검측하기(03.circle_try_hu.m)

Hough 변환_허프만방법으로 원, 직선을 검측하기(02.circle_try.m)

Hough 변환_허프만방법으로 원, 직선을 검측하기(01.circle_product.m)

MATLAB으로 원을 그리고 N등분하기

높이가 다른 동심원 그리기(3차원보기)

동심원 그리기(소스코드)

walsh_demo관련 결과데이터

walsh_demo(하다마드, 왈쉬, 오리지날이미지, 변환이미지, 역변환이미지)

하다마드(Hadamard)에서 왈쉬코드(Walsh Code)만들기(Walsh_xulie)

논문작업: 2D Barcode(dot code A) 테스트 sample만들기

바코드 중국전문가 블로그

대만 국립교통대학공업공정관리교수 사이트

바코드 소개 문서 다운로드하기(특허기준, 설명서)

2D 바코드 만드는 프로그램(barcode-generator)

1차원, 2차원 바코드 만들기 사이트

1D and 2D Barcode(1차원, 2차원 바코드) Generator

1D Barcode Symbologies

2D Barcode Symbologies

匿名连接技术：洋葱路由技术(Onion Routing)

论文阅读：匿名数字指纹的发展方向

论文阅读：匿名数字指纹的研究现状和问题

论文阅读：匿名数字指纹的基本概念

论文阅读：匿名指纹模式的必要性

论文阅读：数字水印和数字指纹的区别

论文阅读：数字指纹方案可以分为三种基本类型

教你电脑自动开机并播放音乐提醒你起床,并且晚上自动关机！

상명대학교 저작권보호학과(Copyright Protection Sangmyung Univ)

포렌식 워터마킹(Forensic Watermarking)

포렌식 워터마킹(Forensic Watermarking)

오디오 워터마킹 알고리즘: MP3, WAV에 워터마크 삽입

runprogram.m(Barcode.jpg)

walsh transform(1D) 왈쉬 변환(소스코드)

핑거프린팅: 평균공모공격 방식 구조도

핑거프린팅 시스템구조도

RSA 암호화알고리즘(Java내부함수를 MATLAB에서 사용)03fo03

RSA 암호화알고리즘(Java내부함수를 MATLAB에서 사용)02fo03

RSA 암호화알고리즘(Java내부함수를 MATLAB에서 사용)01fo03

RSA 암호화 알고리즘(MATLAB 소스코드)

沃尔什编码 （Walsh Code）

Walsh Code(월쉬 부호, Walsh 코드, 왈쉬 코드, 왈시 코드)

Orthogonal Signal, Orthogonal Code(직교 신호, 직교 코드)

Gold Code(골드 코드)

Scrambling Code(스크램블링 코드, 스크램블 코드)

PN Code, PN Sequence(PN 코드, 유사잡음코드, 유사잡음부호, 칩 코드)

沃尔什编码（Walsh Code）