点扩散函数matlab程序

[设计]罚函数法MATLAB程序

一、进退法、0.618法、Powell法、罚函数法的Matlab程序设计罚函数法(通用)

function y=ff(x,k)

y=-17.86*0.42*x(1)/(0.8+0.42*x(1))*(1-exp(-

2*(0.8+0.42*x(1))/3))*exp(-1.6)*x(2)-22.

99*x(1)/(0.8+x(1))*(1-exp(-2*(0.8+x(1))/3))*x(3)+k*(x(2)-

(1.22*10^2*(9517.8*exp(-1

.6-2*0.42*x(1)/3)*x(2)+19035.6*exp(-

2*x(1)/3)*x(3)))/(1.22*10^2+9517.8*exp(-1.6-2

*0.42*x(1)/3)*x(2)+19035.6*exp(-2*x(1)/3)*x(3)))^2+k*(x(3)-exp(-0.8-2*x(1)/3)*x(3)

-exp(-2.4-2*0.42*x(1)/3)*x(2))^2; % 主函数，参数包括未知数的个数n，惩罚因子q，惩罚因子增长系数k，初值x0，以及允许的误差r function G=FHS(x0,q,k,n,r,h,a) l=1;

while

这是做卫星通信的课程设计时涉及的卫星星下点的作图程序,matlab编写。

我们对卫星星下点的代码做了点修改，把地球自转也考虑进去了。

卫星星下点的运行轨迹如下：

卫星纬度050100150200

卫星经度250300350400

代码如下：（红色的为修改部分）

clc

close all

clear all

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

义 % % % % % % % % % % % % % % % % % % % % % %

PI = 3.1415926;

W_EARTH = 7.29e-5;% 地球自转速度(rad)

GGC = 3.986e5;% 地心引力常数

N_T =10;% 实验观测周期数，可延长观察时间

% 轨道六根数

a = 20000;% 轨道长半轴（Km）

e = 0.3;% 轨道偏心率

i = 30 * PI/180;% 轨道倾角（rad）

w = 70 * PI/180;% 轨道近地点幅角（rad）

RAAN = 10 * PI/180;% 升交点赤经

T = 2 * PI * sqrt((a^3) / GGC);

% % % % % 变量定

这是做卫星通信的课程设计时涉及的卫星星下点的作图程序,matlab编写。

% 轨道真近心角f(rad)，

这是做卫星通信的课程设计时涉及的卫星星下点的作图程序,matlab编写。

我们对卫星星下点的代码做了点修改，把地球自转也考虑进去了。

卫星星下点的运行轨迹如下：

卫星纬度050100150200

卫星经度250300350400

代码如下：（红色的为修改部分）

clc

close all

clear all

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

义 % % % % % % % % % % % % % % % % % % % % % %

PI = 3.1415926;

W_EARTH = 7.29e-5;% 地球自转速度(rad)

GGC = 3.986e5;% 地心引力常数

N_T =10;% 实验观测周期数，可延长观察时间

% 轨道六根数

a = 20000;% 轨道长半轴（Km）

e = 0.3;% 轨道偏心率

i = 30 * PI/180;% 轨道倾角（rad）

w = 70 * PI/180;% 轨道近地点幅角（rad）

RAAN = 10 * PI/180;% 升交点赤经

T = 2 * PI * sqrt((a^3) / GGC);

% % % % % 变量定

这是做卫星通信的课程设计时涉及的卫星星下点的作图程序,matlab编写。

% 轨道真近心角f(rad)，

小波去噪举例

MATLAB中用wnoise函数测试去噪算法

sqrt_snr=3; init=231434;

[x,xn]=wnoise(3,11,sqrt_snr,init); % 加噪，信噪比为3 subplot(3,2,1),plot(x)

title('original test function') subplot(3,2,2),plot(xn) title('noised function') lev=5;

xd=wden(x,'heursure','s','one',lev,'sym8');%利用小波对一维信号进行降噪, XD为降噪后

的%信号，CXD，LXD为XD的小波分解结构 % 's' or 'h'决定阈值的使用方式，SCAL决定阈值是%否随噪声变化：'one' 不调整， 'sln'对第一层系%数的层噪声分别进行估计和调整； 'mln'对各层%系数的层噪声分别进行估计和调整；

subplot(3,2,3),plot(xd)

title('One de-noised function')

xd=wden(x,'heursure','s','sln',lev,'sym8'); subplot(3,2,4),pl

电力系统计算潮流 14节点经典例子

ZL=20;

WC=0.01;

ZP=[1,2,0.01335,0.04211,0.0,0.0;

1,3,0.0,0.20912,0.0,0.0;

1,4,0.0,0.55618,0.0,0.0;

1,10,0.05811,0.17632,0.0,0.0340;

1,11,0.06701,0.17103,0.0,0.0128;

2,10,0.05695,0.17388,0.0,0.0346;

2,12,0.0,0.25202,0.0,0.0;

2,14,0.05403,0.22304,0.0,0.0492;

3,4,0.0,0.11001,0.0,0.0;

3,13,0.0,0.17615,0.0,0.0;

4,5,0.03181,0.08450,0.0,0.0;

4,9,0.12711,0.27038,0.0,0.0;

5,6,0.08205,0.19207,0.0,0.0;

6,12,0.09498,0.19890,0.0,0.0;

7,8,0.22092,0.19988,0.0,0.0;

7,12,0.12291,0.25581,0.0,0.0;

8,9,0.17093,0.34802,0.0,0.0;

8,12,0.06615,0.130

Matlab写的函数优化遗传算法程序

function [BestPop,Trace]=fmaxga(FUN,LB,UB,eranum,popsize,pcross,pmutation,options)
% [BestPop,Trace]=fmaxga(FUN,LB,UB,eranum,popsize,pcross,pmutation)
% Finds a maximum of a function of several variables.
% fmaxga solves problems of the form:
% max F(X) subject to: LB <= X <= UB
% BestPop--------最优的群体即为最优的染色体群
% Trace-----------最佳染色体所对应的目标函数值
% FUN------------目标函数
% LB--------------自变量下限
% UB--------------自变量上限
% eranum----------种群的代数,取100--1000(默认1000)
% popsize---------每一代种群的规模；此可取50--100(默认50)
% pcross----

函数逼近与曲线(面)拟合的MATLAB程序

7.1 曲线拟合、误差及其MATLAB程序

例7.1.1 已知函数y f(x) 5x3 14x 7sin2(2 x)和一组数据(xi,yi)列入表7–1中，比较最大误差，平均误差，均方根误差和误差平方和.

表7–1 例7.1.1的一组数据(

x,y)

解 由给定的函数和数据，在MATLAB工作窗口输入

>> x=[-2.5,-1.7,-1.1,-0.8,0,0.1,0.5,3.6]; n=length(x);

y=[-43.50 5.69 11.34 14.16 0 1.02 -6.37 185.84];

f=5.*x.^3-14.*x+7.*(sin(2*pi*x)).^2; fy=abs(f-y);

fy2=fy.^2; [x',y',f',fy',fy2'], Ew=max(fy),

E1=sum(fy)/n, E2=sqrt((sum(fy2))/n), E=sum(fy2)

运行后屏幕显示如下

x y f fy fy2

-2.5000 -43.5000 -43.1250 0.3750 0.1406

-1.

合成孔径雷达成像系统点目标仿真

源程序：

clc

close all

C=3e8; %光速 Fc=1e9; %载波频率 lambda=C/Fc; %波长

%成像区域 Xmin=0; Xmax=50; Yc=10000; Y0=500;

%SAR基本参数

V=100; %雷达平台速度 H=0; %雷达平台高度 R0=sqrt(Yc^2+H^2);

D=4; %天线孔径长度 Lsar=lambda*R0/D; %合成孔径长度 Tsar=Lsar/V; %合成孔径时间

Ka=-2*V^2/lambda/R0;%线性调频率 Ba=abs(Ka*Tsar);

PRF=2*Ba; %脉冲重复频率 PRT=1/PRF;

ds=PRT; %脉冲重复周期

Nslow=ceil((Xmax-Xmin+Lsar)/V/ds);%脉冲数 Nslow=2^nextpow2(Nslow); %量化为2的指数

sn=linspace((Xmin-Lsar/2)/V,(Xmax+Lsar/2)/V,Nslow); %创建时间

信源函数

randerr 产生比特误差样本

randint 产生均匀分布的随机整数矩阵 randsrc 根据给定的数字表产生随机矩阵 wgn 产生高斯白噪声 信号分析函数

biterr 计算比特误差数和比特误差率 eyediagram 绘制眼图 scatterplot 绘制分布图

symerr 计算符号误差数和符号误差率 信源编码

compand mu律/A律 压缩/扩张

dpcmdeco DPCM（差分脉冲编码调制）解码dpcmenco DPCM编码 dpcmopt 优化DPCM参数

lloyds Lloyd法则优化量化器参数 quantiz 给出量化后的级和输出值 误差控制编码

bchpoly 给出二进制BCH码的性能参数和产生多项式

convenc 产生卷积码

cyclgen 产生循环码的奇偶校验阵和生成矩阵

cyclpoly 产生循环码的生成多项式 decode 分组码解码器 encode 分组码编码器

gen2par 将奇偶校验阵和生成矩阵互相转换

gfweight 计算线性分组码的最小距离

hammgen 产生汉明码的奇偶校验阵和生成矩阵

rsdecof 对Reed-Solomon编码的ASCII文件解码

rsencof

MATLAB函数大全（转）

Matlab有没有求矩阵行数/列数/维数的函数？

ndims(A)返回A的维数

size(A)返回A各个维的最大元素个数 length(A)返回max(size(A))

[m,n]=size(A)如果A是二维数组，返回行数和列数

nnz(A)返回A中非0元素的个数

MATLAB的取整函数:fix(x), floor(x) :,ceil(x) , round(x) (1)fix(x) : 截尾取整.

>> fix( [3.12 -3.12])

ans =

3 -3

(2)floor(x):不超过x 的最大整数.(高斯取整)

>> floor( [3.12 -3.12])

ans =

3 -4

(3)ceil(x) : 大于x 的最小整数

>> ceil( [3.12 -3.12])

ans =

4 -3

(4)四舍五入取整

>> round(3.12 -3.12)

ans =

0

>> round([3.12 -3.12])

ans =

3 -3 >>

如何用matlab生成随机数函数 rand