风雨钝剑的小窝

Tag: 个人笔记，CentOS/RHEL的远程桌面服务搭建，原始发布于：RHEL/CentOS上用x11vnc搭建vnc远程桌面服务

个人操作笔记，如有转载，请注明出处。
——by realasking

服务器端配置

1.安装

 bash 
#tar -xzvf x11vnc.tar.gz
#cd x11vnc
#./configure --prefix=/usr
#make
#make install

2.配置

 bash 
vi /etc/xinetd.d/x11vnc

输入以下内容：

 bash 
service x11vnc
{
port = 5900
type = UNLISTED
socket_type = stream
protocol = tcp
wait = no
user = root
server = /usr/bin/x11vnc
server_args = -inetd -o /var/log/x11vnc.log -display :0 -forever -auth /var/gdm/:0.Xauth -localhost -many -bg
disable = no
}

保存，然后执行：

 bash 
#service xinetd restart

编辑/boot/grub/grub.conf，在kernel行后加上5

客户端配置

3.客户端软件

若用putty

Connection-SSH-Tunnels,source port 5900,Destination 127.0.0.1:5900,Local,Auto,

然后点击Add，并保存

若用cygwin-X

假设cygwin安装到（包括必要的x组件）：

E:\cygwin

新建一个快捷方式：CygWin-X，内容是：

E:\cygwin\bin\run.exe /usr/bin/bash.exe -l -c /usr/bin/startxwin.exe -wgl

启动ssh转发

启动putty; 或者启动CygWin-X，并执行ssh username@host -L 5900:localhost:5900

启动TightVNC Viewer

Server栏填写：localhost
Connection profile选择:High-speed network
然后点击connect

参考链接

Tag：个人笔记，原始出处：intel编译器使用笔记

也是以前在百度发布过的老文，如有转载，请注明出处。

——by realasking

ifort相关的库的调用

-lifcore -lsvml -limf

fftw2xc的调用

-lfftw2xc_intel -lmkl_intel_lp64 -lmkl_intel_thread -lmkl_blacs_openmpi_lp64 -lmkl_core -lmkl_def -liomp5 -lpthread -lm

需要同时把编译器和mkl的INC目录均加入到CFLAGS里，把两者的库都加入到LDFLAGS里，需要先在interface里编译接口。

mkl的调用

同上

优化选项

最简单的优化选项

-O2，-O3

针对指令集的优化选项

-xSSE2，-xSSE3，-xSSSE3，-xSSE4.1，-xSSE4.2，-xAVX

几个选项具有排他性，其中AVX是Intel高级矢量扩展指令集

cpu构架选择

-mtune=pentium3、pentium4

默认就是p4了...

过程间优化

禁止单文件过程间优化-no-ip

有些数学功能mkl未提供的情况下的解决办法

在链接库之后加-lm，可以再链接GCC的math库

Tag:个人编写的小程序，原始发布位置在matlab盲去卷积图像复原程序

07年根据老研学论坛上一个网友贴出的文献所写，只具有基本功能，其中paddedsize是dipum函数，如有转载，请注明出处。

——by realasking

代码

 matlab 
function image=bldconv(g)
%Blind-deconv example
%Created by realasking
%Welcome to visit my blog: hi.baidu.com/realasking
%finish time:5-22-2007
 
a1=0.1;
a2=0.01;
 
PQ=paddedsize(size(g));
G=fft2(g,PQ(1),PQ(2));
 
[y x]=size(G);
htemp=ones(3);
h0=freqz2(htemp,PQ(1),PQ(2));
R=Rcreat(y,x);
 
H=h0;
for k=1:10 %计算IMG和psf
iMG=(conj(H).*G)./((conj(H).*H)+a1*R);
H=conj(iMG).*G./((conj(iMG).*iMG)+a2*R);
end
 
IMG=mat2gray(real(ifft2(iMG)));
image=IMG(1:size(g,1),1:size(g,2));
imwrite(image,'testfile.tif')
imshow(image,[]);
hold on;
 
end
 
%===========================
%计算R矩阵的函数Rcreat
%===========================
 
function R=Rcreat(y,x)
%R矩阵生成子函数
%by Realasking
%为bdeconv.m编制
 
i=1:y
j=1:x
 
RI=zeros(y,x);
RJ=zeros(y,x);
R=zeros(y,x);
 
for k=1:y %向量化代码生成R的矩阵
RI(k,i)=-2*cos(2*pi*i./y);
end
 
for k=1:x
RJ(j,k)=-2*cos(2*pi*j'./x);
end
 
R=4+RI+RJ;
 
end

Tag:自己写的一个程序，原始发布在：matlab微分方程循环求解的一种方法，以前matlab功能还不完善的时候有用，现在已经没有什么用处了。

一个取巧的做法，以前自己实现的一个范例，如有转载，请注明出处。

——by realasking

代码

 matlab 
function A=test1
%微分方程求解，由realasking创建
%未经允许请勿转载
%如果经过允许转载时，请保留此信息
%2007年5月13日
q=[5;5;5;4;6;9];
t0=0;tf=20;
x0=[0 0 0 0]
m=size(q)
for i=1:size(q)  
fid=fopen('a.m','w')  
fprintf(fid,'function xdot=a(t,x)\n');  
fprintf(fid,'q=%d;\n',q(i,1));  
fprintf(fid,'xdot=zeros(4,1);\n');  
fprintf(fid,'xdot(1)=(1/45.4)*(1500*(x(2)-x(1))+22000*(x(4)-x(3))+192000*(q-x(3)));\n') ;
fprintf(fid,'xdot(2)=(1/317.5)*(1500*(x(1)-x(2))+(x(3)-x(4))*22000);\n') ;
fprintf(fid,'xdot(3)=x(1);\n') ;
fprintf(fid,'xdot(4)=x(2);\n') ;
fclose(fid) ;
[t,x]=ode23('a',t0,tf,x0)
A(:,:,i)=x;
delete a.m  
end 

Tag:自己写的一个程序，原始发布地址在：matlab处理大学物理实验数据程序之最小二乘法

以前自己根据公式写的，批改实验报告方便，也是以前在百度发布过的老文，如有转载，请注明出处。

——by realasking

代码

 matlab 
function [a,b,r,ua,ub]=line_fit(C,D)
%最小二乘法拟合y=ax+b直线，并计算a、b的A类不确定度
%程序编制：realasking
 
%个人博客：http://hi.baidu.com/realasking
%输入变量:C、D，两个列向量，分别是x，y
%输出变量:a,b,r,ua,ub,r是残差
tmp=(sum(C)^2)-size(C,1)*sum(C.^2);
b=(sum(C.*D)*sum(C)-sum(D)*sum(C.^2))/tmp;
a=(sum(C)*sum(D)-size(C,1)*sum(C.*D))/tmp;
r=D-(a*C+b);
sy=sqrt(sum(r.^2)/(size(C,1)-2));
sa=sy/sqrt(sum((C-mean(C)).^2));
sb=sy*sqrt(((mean(C)^2)/(sum((C-mean(C)).^2))+(1/size(C,1))));
t=[1 12.7;
    2 4.3;
        3 3.18;
      4 2.78;
    5 2.57;
        6 2.45;
      7 2.36;
    8 2.31;
        9 2.26;
      10 2.23;
    15 2.13;
        20 2.09;
      21 1.96]
if(size(C,1)<=10)
    tt=t(size(C,1),2);
    elseif(size(C,1)<=15)
        tt=t(10,2)+((t(11,2)-t(10,2))*((size(C,1)-10)/(15-10)));
  elseif(size(C,1)<=20)
      tt=t(11,2)+((t(12,2)-t(11,2))*((size(C,1)-15)/(20-15)));
else
    tt=t(13,2);
    end
    ua=tt*sa;
    ub=tt*sb;

参考文献

吴平，大学物理实验教程，机械工业出版社