X0=zeros(1,19);
Y0=zeros(1,19);
Z0=zeros(1,19);
i=1;
for t=0:pi/3:6*pi
X0(i)=cos(t);
Y0(i)=sin(t);
Z0(i)=t;
i=i+1;
end
plot3(X0,Y0,Z0,'k-')
hold on
阅读全部
|
gl1218hl1022
贴于 2013年4月9日 21:56
hide
bbsi
fftimage=fftshift(fft2(double(I))); % 对条纹图像 I 进行二维傅里叶变换
figure,plot(mat2gray(abs(fftimage(257,:)))); % 显示频谱图
arrinput = ginput; % 鼠标选取高低频之间的波谷位置，点一下，按回车
arrtemp = fftimage(257,:);
[FringeMax,FringeMaxIndex] = max(arrtemp(arrinput(1,1):512));
% 得到从选取点开始向右的最大值(波峰)，并记录位置
FringeMaxIndex = round(arrinput(1,1) + FringeMaxIndex - 1);
% 得到 1 级谱(基频)的波峰位置
maxvalue=max(fftimage(257,FringeMaxIndex:512));
% 取出 1 级谱的位置。257 表示第 257 行(中间行，谱的位置)
[m,n]=find(fftimage==maxvalue); % 取出 1 级谱最大值位置
fringe_num=abs(n-m); % 确定条纹个数
......................
阅读全部
|
shdlgjmm2007
贴于 2013年3月27日 14:50
hide
bbsi
%-------------------------------------------------------------------------
clear all
x=imread('test_image_9.jpg'); %读取原始图像
figure(1);subplot(1,4,1);imshow(x);title('原图像1');
y=rgb2gray(x);
subplot(1,4,2);imshow(y);title('图像1的灰度图'); % 图1灰度图
u1=imnoise(y,'salt & pepper',0.07);
subplot(1,4,3);imshow(u1);title('图像1加噪声图'); %给图1加椒盐噪声
zz=medfilt2(u1,[3 3]);z=medfilt2(zz,[5 5]);%(2次中值滤波)
subplot(1,4,4);imshow(z);title('图像1中值滤波'); %图1中值滤波（3*3窗口）
% figure(6);subplot(1,2,1);imshow(zz);title('一次滤波');
% subplot(1,2,2);imshow(z);title('二次滤波');
......................
阅读全部
|
guoxiuxiao
贴于 2012年4月9日 10:44
hide
bbsi
%给图像加扰动，扰动方式以zernike多项式的形式加入，Image为源图像，Co为要加的zernike系数，为行向量，
%输出分别为扰动后的图像，扰动PV值，扰动RMS值
function [OUT,dis_pv,dis_rms]=Image_disturb(Image,Co)
Image=imread('1.bmp');
M=size(Image,1);N=size(Image,2);
%--------所加扰动，用zernike多项式表示--------------------
Co=[1,1,1,0.02,0.02,0.02,0.02,0.02,0.02];
N_Zer=size(Co,2);%所用zernike阶数
disturb=zeros(M,N);
Zer=zernike_poly(M,N);%1-36阶
for q=1:N_Zer
disturb=disturb+Co(1,q)*Zer{q};
......................
阅读全部
|
guoxiuxiao
贴于 2012年4月9日 10:33
hide
bbsi
在MATLAB 命令中输入主函数
Syms t
f=t^4-t^2-2*t+5;
[x1,x2] = minJT(f,0,0.1)
进退法确定搜索区间函数文件minJT如下：
function [minx,maxx]= minJT(f,x0,h0,eps)
%目标函数: f ;
%初始点 : x0;
%初始步长: h0;
%精度 :esp;
%区间左端点 : minx;
%区间右端点: maxx;
......................
阅读全部
|
零度元素
贴于 2011年11月29日 17:55
hide
bbsi
%*************** 生成输入信号 *****%
f = 2000 %设定基频率为2000HZ
Fs=44100; %设定采样频率Fs=44.1kHz
N=256; %取的样本点数N
n=0:N-1; %建立矩阵
t=n/Fs; %采样时间t
w = 2*pi*f;
T = 2*pi/w; %求出正弦函数的周期
x1 = sin(w*t);
x2 =awgn(x1,5,'measured'); %加入SNR为5dB的高斯白噪声得到输入信号x2
xn =x2-x1; %间接获得白噪声xn
......................
阅读全部
|
sugeladijq
贴于 2011年10月30日 10:26
hide
bbsi
%*************** 生成输入信号 *****%
f = 2000 %设定基频率为2000HZ
Fs=44100; %设定采样频率Fs=44.1kHz
N=256; %取的样本点数N
n=0:N-1; %建立矩阵
t=n/Fs; %采样时间t
w = 2*pi*f;
T = 2*pi/w; %求出正弦函数的周期
x1 = sin(w*t);
x2 =awgn(x1,5,'measured'); %加入SNR为5dB的高斯白噪声得到输入信号x2
xn =x2-x1; %间接获得白噪声xn
......................
阅读全部
|
sugeladijq
贴于 2011年10月30日 10:26
hide
bbsi