matlab源程序怎么保存

%main_OFDM.m

%一个相对完整的OFDM通信系统的仿真设计，包括编码，调制，IFFT， %上下变频，高斯信道建模，FFT，PAPR抑制，各种同步，解调和解码等模 %块，并统括系统性能的仿真验证了系统设计的可靠性。

clear all close all clc

%++++++++++++++++++++++++++全局变量++++++++++++++++++++++++++++++ % seq_num 表示当前帧是第几帧 % count_dds_up 上变频处的控制字的累加

% count_dds_down 下变频处的控制字的累加（整整） % count_dds_down_tmp 下变频处的控制字的累加（小数） % dingshi 定时同步的定位

% m_syn 记录定时同步中的自相关平台 global seq_num global count_dds_up global count_dds_down

global count_dds_down_tmp global dingshi global m_syn

%+

%main_OFDM.m

%一个相对完整的OFDM通信系统的仿真设计，包括编码，调制，IFFT， %上下变频，高斯信道建模，FFT，PAPR抑制，各种同步，解调和解码等模 %块，并统括系统性能的仿真验证了系统设计的可靠性。

clear all close all clc

%++++++++++++++++++++++++++全局变量++++++++++++++++++++++++++++++ % seq_num 表示当前帧是第几帧 % count_dds_up 上变频处的控制字的累加

% count_dds_down 下变频处的控制字的累加（整整） % count_dds_down_tmp 下变频处的控制字的累加（小数） % dingshi 定时同步的定位

% m_syn 记录定时同步中的自相关平台 global seq_num global count_dds_up global count_dds_down

global count_dds_down_tmp global dingshi global m_syn

%+

书中Matlab源程序 第1 章 绪论

【例1-1】有一名学生，期末有5门功课要考试，可用的复习时间有18小时。假定这五门课程分别是数学、英语、计算机基础、画法几何和专业概论。如果不复习直接参加考试，这五门功课预期的考试成绩分别为65分、60分、70分、60分和65分。复习以1小时为一单元，每增加1小时复习时间，各门功课考试成绩就有可能提高，每复习1小时各门功课考试成绩提高的分数分别为3分、4分、5分、4分和6分。问如何安排各门功课的复习时间可使平均成绩不低于80分，并且数学和英语成绩分别不低于70分和75分。

解：设分配在数学、英语、计算机基础、画法几何和专业概论这五门功课的复习时间分别为x1,x2,x3,x4,x5，则可列出如下的目标函数和限制条件为：

f(x)?x1?x2?x3?x4?x5 x1?x2?x3?x4?x5?18

(3x1?4x2?5x3?4x4?6x5?320)/5?80

3x1?65?70

本例具体程序如下：

%li_1_1 f=[1 1 1 1 1];

A=[1 1 1 1 1; -3 -4 -5 -4 -6; -3 0 0 0 0; 0 -4 0 0 0; 3 0

clear all;

clc;

tic;

X=zeros(1000,4);Y=zeros(1000,4);U=zeros(1000,1);

M=3;

a=[28 18 74 74 70 72 60 36 12 18 14 90 78 24 54 62 98 36 38 32]; b=[42 50 34 6 18 98 50 40 4 20 78 36 20 52 6 60 14 58 88 54]; Aeq=zeros(20,80);

for j=1:20

for k=1:20

if j==k

Aeq(j,k)=1;

else

Aeq(j,k)=0;

end

end

for k=21:40

if k==j+20

Aeq(j,k)=1;

else

Aeq(j,k)=0;

end

end

for k=41:60

if k==j+40

Aeq(j,k)=1;

else

Aeq(j,k)=0;

end

end

for k=61:80

if k==j+60

Aeq(j,k)=1;

else

Aeq(j,k)=0;

end

end

end

Aeq;

A=zeros(4,80);

for i=1:4

for k=1:80

if i==1

for k=1:20

A(i,k)=1;

end

end

if i==2

for k=21:40

A(i,k)

%*************************************************************************************

% This function pertains to the addition of AWGN with mean zero and % parameter 'variance' to an input signal. %

% AUTHOR: Wenbin Luo % DATE : 04/12/01 % % SYNOPSIS: y = awgn(x,var) % x ---> input signal % var ---> variance

% y ---> y = x + AWGN

%***********************************************************************************

function y = awgn(x,var) w = randn(1,lengt

电工专业必看

%简单潮流计算的牛顿拉夫逊程序，相关的原始数据数据数据输入格式如下： %B1是支路参数矩阵，第一列和第二列是节点编号。节点编号由小到大编写 %对于含有变压器的支路，第一列为低压侧节点编号，第二列为高压侧节点 %编号，将变压器的串联阻抗置于低压侧处理。

%第三列为支路的串列阻抗参数。

%第四列为支路的对地导纳参数。

%第五列为含变压器支路的变压器的变比

%第六列为变压器是否含有变压器的参数，其中“1”为含有变压器，

%“0”为不含有变压器。

%B2为节点参数矩阵，其中第一列为节点注入发电功率参数；第二列为节点负荷功率参数；第三列为节点电压参数；第六列为节点类型参数，其中“1”为平衡节点，“2”为PQ节点，“3”为PV节点参数。

%X为节点号和对地参数矩阵。其中第一列为节点编号，第二列为节点对地 %参数。

n=input('请输入节点数:n=');

n1=input('请输入支路数:n1=');

isb=input('请输入平衡节点号:isb=');

pr=input('请输入误差精度:pr=');

B1=input('请输入支路参数:B1=');

B2=input('请输入节点参数:B2=');

X=input('节点号和对地参数:X=');

Y=zeros(n

%*************************************************************************************

% This function pertains to the addition of AWGN with mean zero and % parameter 'variance' to an input signal. %

% AUTHOR: Wenbin Luo % DATE : 04/12/01 % % SYNOPSIS: y = awgn(x,var) % x ---> input signal % var ---> variance

% y ---> y = x + AWGN

%***********************************************************************************

function y = awgn(x,var) w = randn(1,lengt

1、Newdon迭代法求解非线性方程

function [x k t]=NewdonToEquation(f,df,x0,eps) %牛顿迭代法解线性方程

%[x k t]=NewdonToEquation(f,df,x0,eps) %x:近似解 %k:迭代次数 %t:运算时间

%f:原函数，定义为内联函数 ?:函数的倒数，定义为内联函数 %x0:初始值 %eps:误差限 %

%应用举例：

%f=inline('x^3+4*x^2-10'); ?=inline('3*x^2+8*x');

%x=NewdonToEquation(f,df,1,0.5e-6) %[x k]=NewdonToEquation(f,df,1,0.5e-6) %[x k t]=NewdonToEquation(f,df,1,0.5e-6)

%函数的最后一个参数也可以不写。默认情况下，eps=0.5e-6 %[x k t]=NewdonToEquation(f,df,1)

if nargin==3 eps=\end tic; k=0; while 1

x=\ k=\

if abs(x-x0) < eps || k >30

流水线车间调度问题matlab源程序

流水线型车间作业调度问题遗传算法Matlab源码
流水线型车间作业调度问题可以描述如下：n个任务在流水线上进行m个阶段的加工，每一阶段至少有一台机器且至少有一个阶段存在多台机器，并且同一阶段上各机器的处理性能相同，在每一阶段各任务均要完成一道工序，各任务的每道工序可以在相应阶段上的任意一台机器上加工，已知任务各道工序的处理时间，要求确定所有任务的排序以及每一阶段上机器的分配情况，使得调度指标(一般求Makespan)最小。下面的源码是求解流水线型车间作业调度问题的遗传算法通用MATLAB源码，属于GreenSim团队原创作品，转载请注明。



function [Zp,Y1p,Y2p,Y3p,Xp,LC1,LC2]=JSPGA(M,N,Pm,T,P)
%--------------------------------------------------------------------------
% JSPGA.m
% 流水线型车间作业调度遗传算法
% GreenSim团队原创作品，转载请注明
% Email:greensim@
% GreenSim团队主页：/greensim
% 欢迎访问GreenSim——算

计算题

3.6 已知三类训练样本为

?1：[ -1，-1 ]T， ?2：[ 0，0 ]T， ?3：[ 1，1 ]

试用多类感知器算法求解判别函数。

解：采用多类情况3的方式分类，将训练样本写成增广向量形式，有

X1= [ -1，-1，1 ]T， X2= [ 0，0，1 ]T， X3= [ 1，1，1 ]T

任取初始权向量为：

W1(1) = W2(1) = W3(1)= [ 0，0，0 ]T 取校正增量c = 1。 迭代过程如下：

第一次迭代，k = 1，以X1= [ -1，-1，1 ]T作为训练样本，计算得 d1(1) = W1T(1) X1= 0 d2(1) = W2T(1) X1= 0 d3(1) = W3T(1) X1= 0

X1??1，但d1(1)>d2(1)且d1(1)>d3(1)不成立，故修改3个劝向量，即 W1(2) = W1(1) + X1= [ -1，-1，1 ]