通原软件scilab基础实验林承绩

合集下载

scilab使用辅导

scilab使用辅导
信息技术教研组 袁中果
Hale Waihona Puke 国人民大学附属中学练习P22
信息技术教研组 袁中果
中国人民大学附属中学
复杂判断结构
If 表达式 1 语句序列1 elseif表达式 1 语句序列2 else 语句序列3 …… end
信息技术教研组 袁中果
中国人民大学附属中学
练习
给出一个百分制成绩,要求输出成 绩等级‘A’、 ‘B’、 ‘C’、 ‘D’、 ‘E’。90分以上为‘A’ ,80~89分为 ‘B’,70~79分为‘C’,60~69分为 ‘D’,60分以下为‘E’。
信息技术教研组 袁中果
中国人民大学附属中学
三种基本结构
循环结构
While型 Until型
A
A
P1 Y 成立 N成立
P2
不成立
图2.7
信息技术教研组 袁中果
中国人民大学附属中学
选择结构
If 表达式 语句序列1 else
语句序列2
If 表达式 语句序列1 end
end
信息技术教研组 袁中果
中国人民大学附属中学
信息技术教研组 袁中果
中国人民大学附属中学
P7练习A 4
3x-2y+z=-4 5x+2y-3z=12 x+3y+5z=-1 x+y+z=12 3x=3y-z=16 x-y-z=-2
信息技术教研组 袁中果
中国人民大学附属中学
程序设计
变量 在scilab中变量的定义是很自由的,几乎大部 分的字串都可以做变量。例如x=3;a2=14; 但是 变量开头不可以是数字。同時也不可以將Scilab 內定常量更改它的数值。 常量 Scilab的內定常常多半以‘%’開頭,例如圓周 率π,它以‘%pi’來表示。自然常數e:%e。虛 數i:%i。如對這些常量重新定义會出現錯誤信 息。

开源科学计算软件Scilab及其教学应用

开源科学计算软件Scilab及其教学应用

开源科学计算软件Scilab及其教学应用作者:王凤蕊王文宏来源:《中小学信息技术教育》2006年第11期新的《高中数学课程标准》已将“数学建模”课设置为高中数学的一个专题课程。

同时,高中物理、化学等课程的教与学活动也经常需要对各种问题进行数学建模和计算。

因此,科学计算软件的使用能力应成为广大教师和学生的基本技能之一。

开源科学计算软件Scilab不仅能解决各种各样的计算问题,而且能将计算过程和结果可视化,同时还能模拟一些事物的变化过程。

比如,在解析几何的学习中,可以借助Scilab动态可视化生成三维空间,加深学生的理解。

与商业科学计算软件Matlab相比,在功能上Scilab基本可以替代Matlab,且Scilab最诱人之处在于,它是一款开放源码的软件,使用正版软件完全免费。

本文首先对Scilab的功能、特点、下载、安装与编程方法进行了简要介绍,最后结合两个实例探讨了其在中学教学中的应用。

一、Scilab简介Scilab是法国国家信息与自动化研究院(INRIA) 开发的“开放源码”科学计算自由软件,它主要有数值计算、仿真与模拟、计算结果可视化等功能。

Scilab数据类型丰富,可以方便实现各种矩阵运算,并允许用户在线建立各种自定义函数。

此外,Scilab还具有图形显示功能,可实现各种常规形式的计算结果可视化。

Scilab是一种解释性语言,能运行于Windows、Linux以及Unix等操作系统环境下。

作为开放源码的软件,Scilab的源代码、用户手册及二进制的可执行文件都是免费的。

用户不仅可以在Scilab的许可证条件下自由使用该软件,还可以根据需要修改源代码,使之更加符合自身需要。

此外,Scilab还包括一些应用于不同科学计算领域的工具箱,如科学计算、数学建模、信号处理、网络分析、决策优化、线性与非线性控制等。

二、Scilab的安装及运行Scilab的源代码和执行码程序可以从Scilab英文网站/或中文站点/下载,Scilab的安装程序、说明文档、应用案例以及一些最新消息都能从上述网站得到。

用SCIlab软件处理物理化学实验数据

用SCIlab软件处理物理化学实验数据

P= P+ ,; o( 。 2 )
优博 气压平 均值
察温度存 于向
t=7 . 7 . 7 . 8. 8 . 8. 8. 8 . 8 . 9 . l [ 4 66 8 01 21 45 5 7 92 06 4 2 9 2
如果 没有 指定 颜色 ,ca si b内部 会有 一个颜 色 列表 l
图2 a hp与 1 / T的关 系
1 . 首先进行数据输入 和预处理( 具体命令 可参见 Si b ca l
的相 关 文 商 ) 。
t2 . ; o 88 = 5 L= 0 o3 ;
境温度 , 茎起点读数 儡
策优化 、 线性 、 非线性控制等各个方面 , 同时具有 良好的数据可视化功能 ,它可以将科学计算的数
据 用二 维或 三维 图形表 现 出来 。
图、 求斜率等 。传统方法以手工计算为基础 , 主要 在于人工作图处理中存在误差 ,且误差大小因人
而异等不足【 。近几年来相应的物理化学实验教 材 中相继介绍计算机处理相应实验数据的计算机 程序 ,实验数据 的计算机化处理方面的工作有所 报道 ,改进 这类数据处理方法能更大程度地提高 实验教学及科研工作效率。但是大多数软件为商 业软件 , 使用版权受到限制。 本文引入一个开源数 据处理程序 s Ia cL b 处理物理化学实验 数据 , 它 包含各种矩 阵运算 ,能很方便地对实验数据进行 计算并拟合 出图形 ,避免人为的计算误差 和作图 不精确 ,在实验教学 中可以起到不可替代 的辅助 作用 , 可大大提高实验的效率与结果的精确度。
第2卷 第1 6 2期
21年 1 01 2月
乐 山 师范 学 院 学 报
J u a f eh n T a h r C l g o r lo s a e c es ol e n L e

第一讲 Scilab基本操作

第一讲  Scilab基本操作

-->x=0:0.2:1; -->y=x.*x
// 生成当x=0,0.2,0.4,0.6,0.8,1时, 0.16 0.36 0.64
y=x2函数值序列
y =
0 0.04 1.
>> 表达式 y=x.*x中的“.*”表示对应元素相乘; 另外(./)表示对应元素相除,(.^)表示对应元素求幂.
-->x=0:0.2:1; -->y=x.^2 y =
>> 在Command窗口下进行数学运算,只需将表达式直接输入提示号(--> )之后,并按【Enter】 键即可。Scilab会将计算结果直接存入变量ans,代表Scilab运算后的答案 (Answer),并在 Command窗口下显示其数值。 >> Scilab将所有数值均保存为double形式,故不需要变量声明,Scilab会自动进行变量的创建和 回收。 >> Scilab的基本运算符:加(+)、 减(-)、 乘(*)、 除(/)、指数(^); 上述运算符不仅运用于通常实数之间计算,也可运用于矩阵的运算。
Scilab 预定义的一些常数: i: 基本虚数单位 pi: 圆周率(π=3.1415926...) e: 欧拉常数(е=2.7182818...) inf: 表示无限大, 例如1/0 nan或NaN: 表示非数值(Not a number), 例如0/0
注: 在使用预定义常数之前加%,如下图。
【例2-2】简单矩阵的输入和运算。
0
0.04
0.16
0.36
0.64
1.
Scilab常用的基本数学函数
abs(x): 纯量的绝对值或向量的模 sqrt(x):开平方 fix(x): 求实数x的整数部分 modulo(x,y):求x除以y余数 round(x):四舍五入至最近整数 floor(x):不超过x的最大整数 sign(x): 符号函数 当x<0时,sign(x)=-1; 当x=0时,sign(x)=0; 当x>0时,sign(x)=1。 exp(x):自然指数 log(x):以e为底的对数,即自然对数 log2(x):以2为底的对数 log10(x):以10为底的对数 sin(x): 正弦函数 cos(x):余弦函数 tan(x):正切函数 asin(x):反正弦函数 acos(x):反余弦函数 atan(x):反正切函数

科学计算自由软件SCILAB基础教程

科学计算自由软件SCILAB基础教程

第1章 SCILAB简介 1.1 引言1.2 SCILAB软件的构成1.3 安装SCILAB的系统需求1.4 SCILAB主窗口介绍图1.1 SCILAB 3.0的主窗口图1.2 SCILAB 3.0的文件菜单项图1.3 SCILAB 3.0的编辑菜单项图1.4 SCILAB的选择菜单项图1.5 SCILAB的控制菜单项图1.6 SCILAB的编辑器窗口图1.7 SCILAB的应用窗口图1.8 SCILAB的帮助窗口1.5 SCILAB中的基本操作与预定对象 1.5.1 SCILAB中的文件操作1.5.2 界面层次的控制操作1.5.3 SCILAB主窗口中的快捷键操作1.5.4 SCILAB中预先定义的对象图1.9 who指令后的SCILAB主窗口图1.10 whos指令后的SCILAB主窗口1.6 谈谈如何学习SCILAB图1.11 SCILAB的演示窗口图1.12 SCILAB的选择图形演示窗口图1.13 SCILAB的选择图形演示窗口图1.14 SCILAB的一个图形信息窗口图1.15 SCILAB的一个图形演示窗口第2章 数据类型 2.1 引言2.2 特定符号与常数2.2.1 特定符号2.2.2 特定常数2.3 标量的数值计算2.4 数值型向量与矩阵的定义及基本运算 2.4.1 数值型向量与矩阵的定义2.4.2 数值型向量与矩阵的运算2.5 与数值型矩阵有关的若干常用函数 2.5.1 常用矩阵的生成函数2.5.2 size函数和matrix函数2.5.3 从已知矩阵提取部分元素来构成同阶新矩阵的若干函数2.5.4 与方阵的行列式求值、求逆、线性代数方程组的求解、求矩阵特征值与特征向量等有关的函数2.5.5 与矩阵(数组)或向量有关的数据统计函数2.6 向量与矩阵(数组)元素的引用2.7 整型数组2.8 布尔型向量与矩阵的定义及基本运算2.9 字符串型数据的定义与运算 2.9.1 字符串的定义2.9.2 字符串的运算2.9.3 与处理字符串有关的一些常用函数2.10 多项式类型 2.10.1 多项式的定义2.10.2 多项式运算2.10.3 有关多项式的几个常用函数2.11 表类型2.11.1 表类型的定义2.11.2 表类型数据的引用,域的插入与删除2.11.3 tlist和mlist类型第3章 SCILAB中的程序设计、脚本文件与函数 3.1 引言3.2 顺序结构程序设计3.2.1 赋值语句3.2.2 输入输出语句3.3 选择结构程序设计 3.3.1 if语句图3.1 不完整if语句的框图图3.2 完整if语句的框图3.3.2 select_case语句图3.3 select 语句框图(else缺省情形)3.4 循环结构程序设计 3.4.1 for语句图3.4 for语句的框图3.4.2 while 语句图3.5 while 型循环框图3.4.3 循环语句的嵌套(多重循环语句)3.4.4 continue 语句和break 语句图3.6 循环体内含有continue语句的框图图3.7 循环体内含有break语句的框图3.5 脚本文件与函数3.5.1 脚本文件3.5.2 函数3.5.3 局部变量与全局变量、函数的嵌套定义与递归调用图3.8 用递归调用求 n!示意图(n =4)3.5.4 有关程序调试的几个常用指令3.6 函数的应用3.6.1 函数名作形式参数——二分法求非线性方程的根图3.9 用二分法求方程p(x)=0的根的框图3.6.2 函数的递归调用——求两正整数的最大公因子图3.10 用递归调用求n!示意图3.6.3 多重循环的应用——线性方程组的顺序消元法图3.11 顺序消元法的框图第4章 计算结果可视化 4.1 引言图4.1 图形窗口图4.2 图形窗口的文件菜单图4.3 图形窗口的工具菜单图4.4 图形窗口的编辑菜单4.2 二维图形的绘制4.2.1 plot指令图4.5 利用指令plot(y,'X', 'Y', 'y=f(x)')画图。

《讲Scilab基本操作》课件

《讲Scilab基本操作》课件
《讲Scilab基本操作》PPT 课件
通过这个PPT课件,我们将学习Scilab的基本操作,让你在科学计算方面快速 上手。准备好探索Scilab的世界吧!
一、Scilab的介绍
Scilab的概述
了解Scilab是一个开源的数 值计算和数据分析软件, 适用于工程学、科学研究 和教育领域。
Scilab的特点
探讨Scilab在不同领域的应用 前景,如工程学、科学研究 和数据分析。
Scilab的学习建议和 总结
提供学习Scilab的建议和总结, 帮助大家更好地掌握Scilab的 基本操作。
Scilab的数据分析与拟 合
掌握如何使用Scilab进行数 据分析和曲线拟合,以从 数据中获得有用信息。
Scilab与其他编程语言 的集成
了解如何将Scilab与其他编 程语言集成,从ቤተ መጻሕፍቲ ባይዱ扩展其 功能和应用领域。
五、Scilab常见问题解答
1 Scilab的常见错误及解决方法
指导解决常见的Scilab错误,如语法问题、函数调用和脚本调试。
2 Scilab的优化与调试技巧
分享Scilab的性能优化和调试技巧,确保代码的正确性和效率。
3 Scilab的相关资源推荐
为学习Scilab提供相关资源推荐,如文档、教程和论坛。
结束语
Scilab的未来发展
展望Scilab的未来发展方向, 包括新功能、性能改进和社 区贡献。
Scilab的应用前景
掌握Scilab的特性,如多维 数组、函数库、图形显示 和内置编辑器等。
Scilab的优势
了解为什么Scilab是一种强 大的工具,可以用于解决 复杂的数学和工程问题。
二、Scilab的安装与环境配置

科学计算自由软件SCILAB在考试分析中的应用

科学计算自由软件SCILAB在考试分析中的应用
总 项 分 布 图
43 不 同学校 同一科 目的成 绩指标 的 比较 . 教 育研究 者经 常需要 对 同一科 目不 同考试 群体 成绩进 行深 入 比较 , 就需 要对考 试 成绩 的计算 考 这
试 成绩 的各项 指标 。表 1 贵州省 贵 阳市两所 学校 是 20 0 9年 中考 语文 成绩 的两项 统计 指标 量 :
作者简 介: 曲霄2 ̄ 9 0 ) 女 ,山西应县人 , 1 18 一, ( 硕士 , 讲师 , 研究方 向:网络通信 与信息 安全 。
2 1 年 01
f 霄 红 : 学 计 算 自由软 件 S L B在 考 试 分 析 中 的应 用 抖 1 科 CIA
・5 1
全距= 最大 值一 小值 最 伞距 简 明易 懂 , 算 简单 ,但它 容 易受 到 两极 计 端数 据值 的影 响 ,并且 只能反 映 数据 组 的两 极端 值 的 差 异 程 度 ,不 能 很 好 地 反 映 中 间 数 据 的差 异 情 况 , 映不够 灵 敏 。 因此 , 只能 粗 略地反 应数 据 反 它 的差 异程 度 ,只有 与 其他 统 计 量结 合 起来 使 用 ,才 能 比较 全 面地 反映 出数 据分 布 的差 异情 况 。
关关 系 等 。
图 1 频 数分 布直 方 图
42 同一科 目的不 同学校 的 成绩 差异 比较 _ 图 2是 根 据 贵 州 省 贵 阳 市 2 0 0 9年 中考 各 区
3 数 据 分 布 的 图表 表 示
数据 分 布 是指 反 映观 察 值 及 其呈 现 规 律 的 表 、
表 1 统 计 指 图 3是根据上表 中指标 绘制的正态分布 曲线 , 从
5 结 束 语
由 于 科 学 计 算 自 由软 件 S IAB可 以 免 费 下 CL

物理实验技术中的实验数据处理与分析工具推荐

物理实验技术中的实验数据处理与分析工具推荐

物理实验技术中的实验数据处理与分析工具推荐引言:在物理实验中,获得准确和可靠的实验数据并进行合理的处理和分析是至关重要的。

为了提高数据处理的效率和准确性,科学家和研究人员们积极探索并开发各种实验数据处理和分析工具。

本文将为大家推荐几个常用的物理实验数据处理和分析工具,并简要介绍其主要特点和应用。

一、MATLAB(Matrix Laboratory)MATLAB是一款功能强大的科学计算软件,特别适用于处理实验数据。

它提供了丰富的数据分析和处理函数,如数值计算、数据可视化、信号处理、图像处理等。

对于物理学研究者而言,MATLAB可以用来处理实验数据,进行统计分析和回归分析,帮助研究者发现规律并评估实验结果的有效性。

二、PythonPython是一种高级编程语言,具有优雅、简洁和易学的特点。

它广泛应用于科学计算和数据分析领域。

在物理实验中,Python可以利用众多强大的第三方库,如NumPy、SciPy和Pandas等,进行数据处理与分析。

相比于其他工具,Python 具有较高的自由度和可扩展性,能够满足不同实验需求的要求。

三、OriginOrigin是一款非常受欢迎的科学绘图和数据分析软件。

它提供了丰富的数据处理和分析工具,如傅里叶变换、线性回归、非线性拟合等,并且有着直观的图形用户界面。

Origin既适用于初学者进行简单的数据处理,也适用于专业研究者进行高级的数据分析和可视化。

四、ExcelExcel是一款常见的电子表格软件,广泛应用于数据处理和分析。

对于简单的实验数据处理,Excel提供了一系列基本的函数和工具,如排序、过滤、求和等。

此外,Excel还支持绘制简单的图表,如折线图、柱状图等。

虽然Excel功能相对有限,但它易于学习和操作,尤其适合初学者或对数据处理要求较低的实验。

五、LabVIEW(Laboratory Virtual Instrument Engineering Workbench)LabVIEW是一款用于测量和控制系统的开发环境,也常用于实验数据处理和分析。

系统分析与设计 (2)PPT课件

系统分析与设计 (2)PPT课件

需求分析
系統分析 與設計
系統實施

系統分析與設計─理論與實務應用(三版)
圖1-1b 七階段資訊系統 開發步驟
使用者 需求分析
軟硬體 需求分析
系統分析
系統設計
編碼
測試
操作與 維護

系統分析與設計─理論與實務應用(三版)
資訊系統開發的原則
人類為滿足需求而建立了各種系統。例如,萬 里長城是防禦系統、電影是娛樂產品,而企業 之電腦軟體則是管理資訊系統。

系統分析與設計─理論與實務應用(三版)
Course Outline
上課時請往前座入坐。 上課時將手機置於震動 (以聽不到鈴聲為原則)。 上課時請勿干擾臨座同學聽課或本人上課,制止不聽
者,扣分並請離教室。 期中、期末測驗各佔30%成績,其中,期中測驗以筆
試方式,期末以簡報方式報告展示系統。 平時成績計算方式:
資訊系統開發要考慮人、方法、科技及企業經營 等四個重要的構面,其關係有如鑽石結構,各構 面需要均衡發展: 人
方法
科技 企業經營

系統分析與設計─理論與實務應用(三版)
資訊系統開發的特性與挑戰
資訊系統開發的特性如下: 軟體是人造的邏輯系統。設計者對於電腦程 式中的邏輯必須完全瞭解。 軟體開發是一種專業的工作。軟體開發者之 專業門檻較高,且養成時間也較長。 軟體開發是一種智力活動。人的智力潛能無 窮、個別差異大,同時也充滿著錯誤和偏 誤。
1. 平時考 (第八週、第十二週實施)筆試、需求分析報 告、初步系統分析報告30% 2. 出席10% 3. 平時考未按時參加者、輪到未按時報告者,以0分計 4. 一次點名未到扣2分,遲到或早退減半扣

系統分析與設計─理論與實務應用(三版)

scilab-使用教程

scilab-使用教程

与MATLAB的互操作性
01
SCILAB可以与MATLAB进行互 操作,允许用户在SCILAB中调 用MATLAB函数或编辑器中的 脚本。
02
互操作性可以通过MATLAB Engine API for Scilab实现, 该API提供了与MATLAB的接 口,以便在SCILAB中执行 MATLAB代码。
图像处理
SCILAB提供了图像处理和分析的功能,可以用 于图像处理和计算机视觉等领域。
ABCD
工程仿真
SCILAB可以用于工程设计和仿真的数值计算和 分析,如流体动力学、结构力学等。
系统控制
SCILAB可以用于控制系统设计和分析,如控制 系统仿真和优化等。
02 SCILAB基础操作
启动与退出SCILAB
易用性
SCILAB的用户界面友好,易于学 习和使用,支持多种编程语言接 口,方便用户进行编程和开发。
开放性
SCILAB是开源软件,用户可以自 由获取和使用源代码,并根据自 己的需求进行定制和开发。
SCILAB的应用领域
科学计算
SCILAB广阵元素
使用单引号或双引号将元素索引括起来,例如 `A(1,1)`或`A"[1,1]"`。
矩阵运算
支持常见的矩阵运算,如加法、减法、乘法 等。
03 SCILAB编程基础
控制流语句
顺序执行
SCILAB中的代码按照从上到下的顺序 执行,无需使用特定的控制流语句。
循环语句
使用`for`循环重复执行特定代码块, 例如`for i=1:10, statements, endfor`。
初值条件
求解常微分方程时需要指定初值条件,可以使用"="指定初值。例如, y(t0)=y0将指定在t0时刻的初值。

开源软件Scilab在自控原理教学中的仿真应用

开源软件Scilab在自控原理教学中的仿真应用
口, 用 户 可 以 自己独 立编 写 自己需要 的 T o o l k i t , 开发 新 应用很 方便 , 值得推广 和使用 。
二、 基于 S c i l a b的编 程
i p i v , &u n ) ; }
图 1 C代 码 范例
( 二) 在S c i l a b里编 译 、 连 接 代 码
{ c o n s t i n t u n= 1 ; C 2 F( d l a s w p ) ( n , A, i n , &u n , mi p i v ,
到理解和掌握课程 的相关理论 很难 。 有 些教师尝试在课 程 中引 入 Ma t l a b软件 ,但 是 M a l f a b是商业工 程计算 软
如下所述 : ( 一) 编写基 于 c语 言的专用库
生成动态库写代码 是第一步 , 编译连接 代码是第 二
步, 要 想在 S e i l a b 里使 用这些库 ( 工具箱 ) 还需 要进行第
三 步就是 加载共 享库。如果选用 的系统平 台是 Wi n d o w s
比较 多 、 比较 抽象 , 因此 多数 同学在 学 习过程 中都会 感
e x t e r n i n t C 2 F ( d l a s w p ) ( i n t 六n , d o u b l e ★A, i n t ★I d a ,
i n t ★k l , i n t ★k 2 , i n t ★i p i v , i n t ★s t e p ) ; i n t d l a s w p l( d o u b l e ★A, i n t ★I n , i n t ★n , i n t ★m i p i v , i n t i p i v )

Scilab程序设计教程_v0.02

Scilab程序设计教程_v0.02

Scilab程序设计教程内容摘要本书是为适应教学的需要而编写。

全书以Scilab4.1版为基础,介绍Scilab程序设计的方法与应用,主要内容包括:Scilab操作基础、Scilab矩阵及其运算、Scilab程序设计、Scilab 文件操作、Scilab绘图、Scilab数据分析与多项式计算、Scicos图形化仿真环境以及Scilab 常用工具箱。

本书符合教学规律和特点,内容丰富、取材新颖、注重实用、便于教学。

本书可作为高校理工科各专业大学生、研究生学习的教材,也可供广大科技工作者阅读使用。

编者列表:序前言目录1. Scilab基础 (4)1.1. Scilab概述 (4)1.2. Scilab的运行环境与安装 (4)1.3. Scilab集成环境 (6)1.4. Scilab帮助系统 (8)2. Scilab矩阵及其运算 (8)2.1.变量和数据操作 (8)2.2. Scilab矩阵 (11)2.3. Scilab运算 (15)2.4.字符串 (16)2.5.稀疏矩阵 (20)3. Scilab程序设计 (17)3.1. SCI文件 (17)3.2.程序控制结构 (17)3.3.函数文件 (18)3.4.程序举例 (18)3.5.程序调试 (18)4. Scilab文件操作 (19)4.1.文件的打开与关闭 (19)4.2.文件的读写操作 (19)4.3.数据文件定位 (19)5. Scilab绘图 (19)5.1.二维数据曲线图 (19)5.2.其他二维图形 (19)5.3.三维图形 (19)5.4.图形修饰处理 (19)5.5.图像处理及动画制作 (19)6. Scilab数据分析与多项式计算 (20)6.1.数据统计处理 (20)6.2.数据插值 (20)6.3.曲线拟合 (20)6.4.离散傅立叶变换 (20)6.5.多项式计算 (20)7. Scilab解方程与函数极值............................................................................. 错误!未定义书签。

开源软件scilab在大学物理教学中的应用研究

开源软件scilab在大学物理教学中的应用研究

开源软件scilab在大学物理教学中的应用研究
陈永冰
【期刊名称】《吉林省教育学院学报:中旬》
【年(卷),期】2013(029)009
【摘要】本文列举大学物理中的典型实例,介绍开源软件scilab在大学物理教学中的具体应用。

在大学物理课程中适当地加入scilab图像和动画,推动基础课程教学方法的进步,可以使物理教学变得生动并更有科学性。

【总页数】2页(P68-69)
【作者】陈永冰
【作者单位】贵州师范学院物理与电子科学学院,贵州贵阳550018
【正文语种】中文
【中图分类】G641
【相关文献】
1.开源软件Scilab在自控原理教学中的仿真应用 [J], 董南萍;廖文江;周冠玲
2.开源软件Scilab在自控原理教学中的仿真应用 [J], 董南萍;廖文江;周冠玲;
3.开源软件Scilab在自控原理教学中的仿真应用 [J], 董南萍; 廖文江; 周冠玲
4.虚拟仿真技术在大学物理课程教学中的应用研究与展望 [J], 杨亚玲;马驰
5.微课和翻转课堂在大学物理实验教学中的应用研究 [J], 高银浩;梁明超;曹永华因版权原因,仅展示原文概要,查看原文内容请购买。

应用Visual Basic 6.0开发木材标本管理系统

应用Visual Basic 6.0开发木材标本管理系统

应用Visual Basic 6.0开发木材标本管理系统
徐斌;许成林
【期刊名称】《木材工业》
【年(卷),期】2005(019)003
【摘要】应用Visual basic 6.0语言开发了安徽农业大学木材标本馆标本管理系统,系统设有检索、输入、统计、打印、维护等功能模块.本系统适当调整后也适用于林业系统的树木学、植物、动物、昆虫等标本的管理.
【总页数】3页(P37-38,43)
【作者】徐斌;许成林
【作者单位】安徽农业大学森林利用学院,合肥,230036;安徽农业大学森林利用学院,合肥,230036
【正文语种】中文
【中图分类】S781;TP312VB
【相关文献】
1.利用Visual Basic 6.0开发人员管理系统 [J], 潘波滢
2.基于Visual Basic 6.0及Access 2003数据库技术开发资料室图书管理系统 [J], 宋晓英;李丹
3.Visual Basic 6.0在高校教务管理系统中的应用 [J], 黄华;付铮;陈轩
4.Visual Basic6.0开发儿童影楼管理系统初探 [J], 武帅
5.基于Visual Basic
6.0开发的《倒闸操作危险点预控管理系统》 [J], 李智
因版权原因,仅展示原文概要,查看原文内容请购买。

通原实验报告labview

通原实验报告labview

THE COMMUNICATION FUNDAMENTALSEXPERIMENTNAME: 单赟吉STUDENT ID:11211105CONTENTSLEB2 digital communication (3)1.Experiment Principle (3)2.Experiment Purpose (5)3.Experimental Content (5)4.Experimental Block Diagram (6)PART1:The UART communication system (7)PART2:An Image over an LTI Channel and Dealing with Non-Ideality (8)5. Experimental Results (9)LEB5Altitude modualtion (11)1.Experiment Principle (11)2.Experiment Purpose (13)3.Experimental Content (13)4.Experimental Block Diagram (14)PART1:AM transmitter (14)PART2:A Single-/Multi-tone AM Transmitter (14)PART3:A *.wav AM transmitter (15)PART4:AM receiver (16)5. Experimental Results (17)LAB 2DIGTIAL COMMUNICATION1.Experiment Principle1.1 UART protocolAs aforementioned, UART is widely used for digital chip-to-chip communications on embedded systems. While there are many configurations , a common UART packet structure is 10 bits- 1 START bit, 8 DATA bits (one byte), and 1 STOP bit. This structure is shown in figure 1 below, and is commonly referred to as 8n1 (8 data bits, no parity, 1 stop bit). Sometimes an additional bit is included for parity checking (a simple form of error detection), but that is not covered in this lab.The data byte is sent LSB (Least Significant Bit) first, which means we effectively send the byte "backwards." For instance, if we wanted to send the byte 00101001, we would transmit the following 10-bit sequence: 010******* (one start bit (0), the 8 data bits LSB->MSB (10010100), and then one stop bit (1)).These bits are encoded as voltages: a HIGH voltage (i.e. 5 volts) to represent a 1,and a LOW voltage (i.e. 0 volts) to represent a 0. The START bit is always a LOW state, and the STOP bit is always a HIGH state. Each bit state is held for someperiod of time - the symbol period - to generate a ’square-wave’ type of signal. The amount of time for which we hold a LOW or HIGH bit determines the baud rate (speed) of the link. In most systems, the transmitter and receiver agree upon thisbaud (or bit-clock) rate in advance. The resulting packet signal is sent along a wire from the transmitter to the receiver.UART packets can be sent from the transmitter at any point in time , and there canbe any amount of time in between packets, hence the asynchronous nature of this communication link. In order to achieve synchronization with an incoming packet,the communication wire idles in the HIGH (1) state in between packets. Since the START bit is always a LOW, we know a packet has begun when this transition occurs. After we synchronize to the start of a packet, we use the known baud rate to estimate the center of each data bit, and sample the voltage of the signal at this point. Infigure 2 below, we illustrate the start of a packet and the data-bit sampling interval. After the receiver decodes the entire data packet, we reverse the result (to get theoriginal MSB->LSB) byte, and we’re done! We’ve received one UART packet. The stop bit simply returns the communications wire to the original IDLE (HIGH) state,and the receiver begins waiting for the next START bit which signals the beginning of the next packet.1.2 UART decoderThe transmitter encodes 0’s as0Volts and 1’s as 5Volts. Each character of the message is represented using ASCII encoding. ASCII is a one-byte(8 bit) representation for 256 of the most of the commonly used characters in theEnglish language. As a result, each UART packet you will receive consists of a singlecharacter from this message. The characters of the message are sent in order, so theyshould be kept in the order that they are received. After decoding the entire stream, you will have received a question.1.3Estimating an LTI channelT he additional ’byte’ (array element)which is pre-pended to the picture array. It contains decimal 85, which isequivalent to binary 01010101. Your entire first UART packet, including the startand stop bit, will then be 010******* (remember the data order goes LSB->MSB).We’ll call this packet the synchronization & training packet. Thefollowing additional pieces of information:1. Before the start of the first UART packet, the transmitter sends ’11111111’.This is simply to ensure that we begin with the communication line IDLEing HIGH.2. Following these first eight 1’s, we send the aforementioned "synchronization" packet.3. The LTI channel impulse response is no more than 8 taps long (h[i] = 08 i>7), and the channel is noise-free.4. The first tap of the channel impulse is guaranteed to be non-zero (h[0] 6= 0)5. Before the UART transmission starts, the state of the communication line isunknown ( x[i] unknown 8 i< 0 )6. After this first packet, the communication continues as usual: it sends the pixels of the image.1.4Correcting for the channelWe will consider a simple example to demonstrate how you are going to ’correct’ for the channel: Consider the case where you have a channel which has a response of length 2, and we know the first tap is non-zero. The output y[n] at any given time is thenNow, let’s assume we successfully built a channel estimation block to determine h[0] and h[1]. Let’s also assume that we know x[n-1], since that was the previous symbol we should have already decoded (think about this for a moment - it’s a recursive argument). Since we know everything besides x[n], we just re-arrange the above equation to get an estimate of the current sample:2.Experiment PurposeIn Lab 2, our goal is to know some aspects of digital signaling, including simple channel codes and the effect of sending these signals over non-idealcommunication channels.3.Experimental Content3.1 Part1:The UART communication systemIn this lab, the main mission is the decoder. We can put it into 3 steps. Step 1: Build the Receiver State Machine. Step 2: Decode the incoming sample stream. Step 3: Read the question!3.2 Part2:An Image over an LTI Channel, and Dealing withNon-IdealityThe goal for this part of the lab is to know somethingabout non-ideal channels and thetheory of deconvolution. First, Verify that your receiver still works. Then Estimating an LTI channel. Finally, Correcting for the channel.4.Experimental Block DiagramPART1:The UART communication systema.The whole block diagramThe ASCII transmitter and the converter have already finished, thus we only build the UART receiver system to complete the communication system.b.the UART receiverThe diagram given is too simple,I must do a lot of work to finish it. But there is the same receiver in the Lab 2 which has finished. So I just link every modules. Finally, the whole system of the receiver is as blow. I will show you how it works in the following paper.c. Decision regenerationBefore into the state machine, the signal must be identified as high electrical level or low. We need a decision regeneration module to choose the state of the code.d.The addition of 1D arrayIn the READ state we must add the array after decimatinge. The RegisterWe also need a register to save the data after the decoding and output all at last to get the resultf. The IDLE stateThis link form is waiting to see a transition from HIGH->LOW. On transition, go to state READ.PART2:An Image over an LTI Channel, and Dealing withNon-IdealityThe UART receiver has already built up, the mission is to complete the channel correction.The whole block diagram of the channel correction is shown as this figure.5. Experimental ResultsPART1: the receiver can receive the code “How are you?”by recovering the ASCII form the transmitter.LAB5Amplitude Modulation1.Experiment Principle1.1 Amplitude ModulationAmplitude Modulation (AM) is an analog modulation scheme where the amplitude (A) of a fixed-frequency carrier signal is continuously modified to represent data in a message. The carrier signal is generally a high frequency sine wave used to “carry” the information on the envelope of the message. The result is a double-sideband signal, centered on the carrier frequency, with twice the bandwidth of the original signal. The mathematical expression of modulated wave is:Here,A0represents DC component of the modulation signal m(t), f(t)represents AC component of the modulation signal, the carrier wave we use has unitamplitude,angular frequency is a fixed valuewc, the initial phase of the carrier signal c(t) is θ0.From(1.1),we can conclude that Amplitude Modulation is the product of modulation signal m(t) and carrier signal c(t). In order to realize that the amplitude of carrier signal coordinated by f(t) linearly, m(t) should consist of DC component to make surem(t)≥0, i.e.,Therefore, the envelop ofSam(t) is completely above time axis, as shown in figure 1.Figure.1 Time domain of an AM signalFor we can obtain its spectrum byFourier transform directly as shown in Figure.2:There are three fundamental AM methods known as DSB (Double Side-Band),SSB(Single Side-Band),VSB( Vestigial Side-Band).1.2 Amplitude DemodulationSince the core idea of the software radio is try to attach the analog-to-digital converter to an antenna and get the software as close to the antenna as feasible, the process of various modulation and demodulation towards the source signal is carried out in digital domain.PC terminal sends a source signal to the general software radio platform USRP - 2920 via gigabit Ethernet interface, and then use two high speed D/A converter to up convert source signal , finally manipulate the orthogonalization I/Q modulation after low pass filtering. Among them, the I/Q represents the in-phase component (I) and quadrature component (Q), and the resulting I/Q sampling points after pa eventually sent to the transmitting antenna a gigabit Ethernet interface is sent to the PC, the whole process of signal processing in LabVIEWFor a single tone signal mathematical expression of amplitude modulated signal can be written as:Among them, 0A is DC component of modulated signal m(t), βam represents the modulation index. We can obtain in-phase and quadrature components after the orthogonal decomposition:In-phase componentQuadrature componentThe square root of the sum of in-phase and quadrature components is:where, m(n) denotes modulated signal which has be digitalized , ϕ0 denotes the sampling frequency of I/Q signal .We can get the modulated signal m(n)by subtract the DC component A0 from (1.5). This method has stronger ability to resist frequency mismatch, in other words, a certain frequency deviation is tolerant between the local carrier and carrier of modulated signal. When there is a frequency or phase offset between the local carrier and carrier of modulated signal due to transfer signals or other reasons, the in-phase component and quadrature component can be represented as:Where,represents frequency offset and phase offset respectively, which can be constant or random variable;whats more, ωLO and φLO is frequency and initial phase angle value of the local carrier. We can also obtain formula(1.13)by extracting the quadratic sum of them. Therefore, it’s not obligatory that frequency and phase of signal’s carrier is synchronization with the local crystal carrier strictly when we adopt orthogonal demodulation algorithm to demodulate AM signals.2.Experiment PurposeFamiliar with the principle of the signal modulation, the transmittion of single Tone signal, the system of the wav player and the receiver.3.Experimental Content3.1 Amplitude Modulation SimulationWe need tobuild a VI which implements the longer of the two equations shown above for Amplitude Modulation. Inspect the front panel and block diagram that has already been created for you. When this VI is completed, We will be able to select the amplitude and frequency of both the carrier and data signals as well as see the time and frequency domain representation of the signals.3.2A Single-/Multi-tone AM Transmitterwe will generate a real single tone signal, modulate it by its amplitude and transmit it on VHF band after up-conversion based on USRP 2920.3.3A *.wav AM TransmitterOpen AM_Tx_WAV.vi which consists of the basic transmitter side code. The basic structureof the code is in place, I must fill in the missing components.3.4AM Receiver4.Experimental Block DiagramPART1:AM TransmitterJust build the module and link each one follow the instructor.The block diagram is just likeThe system realized the function of AM modulation by the equationPART2: A Single-/Multi-tone AMTransmitterwe will generate a real single tone signal, modulate it by its amplitude and transmit it on VHF band after up-conversion based on USRP 2920. Unfortunately, we cannot use the AM modulated signal in chapter 2.1 as a subVI of AM transmitter directly. A motherboard of USRP product family provides the following subsystems: clock generation and synchronization, FPGA, ADCs, DACs, host processor interface, andpower regulation."TX parameters:" This is a control which sets the parameters of the USRP.(a) "TX IQ Rate" is the number of samples per second that are transmitted by the USRP. Make sure this is set to the default value of 200k .(b) "TX Carrier Frequency" is the center frequency at which the transmitter transmits. (c) "TX Gain" gives you a control over the power at which to transmit. Typically for close range communication, the default value of 0 dB works fine.(d) "TX Antenna" is the antenna being used for transmission. Make sure this is set to "TX1".USRP programming partincluding the open USRP channel, parameter configuration, begin to collect, obtain baseband waveform data on the frequency of the waveform data, read into the niUSRP Write Tx Data VI, the While loop to read waveform modulated signal. Finally stop and close the USRP, the release of resources.Here we use ssb:()()()t m t m k t s DSB max 1∙+=PART3:A *.wav AM TransmitterUSRP programming including the open USRP channel, parameter configuration, begin to collect, obtain baseband waveform data on the frequency of the waveform data, read into the niUSRP Write Tx Data VI, the While loop to read waveform modulated signal. Finally stop and close the USRP, the release of resources.PART4:AM receiverThe whole block diagram is:USRP programming section as shown, including open USRP receive path, parameter configuration, start collecting, continuous access to the baseband waveform data which after down-converted and read out of these into the right side of the baseband IQ shift register, While loop initialized the left side of the shift register . Finally, stop and turn off the USRP, release resources.The baseband waveform of the envelope detection processing is in the middle, first, the complex basebandwaveform is converted to polarcoordinates, then find the amplitude of the envelope, and then stored in the waveform of the shift register to the right, the small box to the left of the shift register initialized.5.Experimental Results5.1Amplitude Modulation Simulation5.2A *.wav AM TransmitterAfter setting the data of the frequency and input the USRP IP address, we can receive a single tone signal by USRP.5.3AM receiverwe can find a AM radio station by setting the RX carrier Fre by USRP。

北邮scilab通信原理软件实验报告

北邮scilab通信原理软件实验报告

信息与通信工程学院通信原理软件实验报告实验二时域仿真精度分析一、实验目的1. 了解时域取样对仿真精度的影响2. 学会提高仿真精度的方法二、实验原理一般来说,任意信号s(t)是定义在时间区间(-无穷,+无穷)上的连续函数,但所有计算机的CPU 都只能按指令周期离散运行,同时计算机也不能处理这样一个时间段。

为此将把s(t)按区间[-T/2 ,+T/2 ]截短为按时间间隔dert T均匀取样,得到的取样点数为N=T/dert T.仿真时用这个样值集合来表示信号s(t)。

Dert T反映了仿真系统对信号波形的分辨率,越小则仿真的精确度越高。

据通信原理所学,信号被取样以后,对应的频谱是频率的周期函数,其重复周期是1/t;。

如果信号的最高频率为那么必须有才能保证不发生频域混叠失真,这是奈奎斯特抽样定理。

设则称为仿真系统的系统带宽。

如果在仿真程序中设定的采样间隔是,那么不能用此仿真程序来研究带宽大于这的信号或系统。

换句话说,就是当系统带宽一定的情况下,信号的采样频率最小不得小于2*Bs,如此便可以保证信号的不失真,在此基础上时域采样频率越高,其时域波形对原信号的还原度也越高,信号波形越平滑。

也就是说,要保证信号的通信成功,必须要满足奈奎斯特抽样定理,如果需要观察时域波形的某些特性,那么采样点数越多,可得到越真实的时域信号。

三、实验步骤1.将正弦波发生器模块、示波器模块、时钟模块按下图连接:四、实验结果对比分析时钟设置0.01,得到的结果如下:时钟设置0.3,以后得到的结果如下:五、思考题(1)观察分析两图的区别,解释其原因。

答:因为信号周期是1,而第一个图的采样周期是0.01,所以一个周期能采样100个点,仿真出来的波形能较精确地显示成完整波形,而第二个图采样周期是0.3,所以一个周期只有三个采样点,故信号失真了。

(2)将示波器的控制时钟的period的参数改为0.5,观察仿真结果,分析其原因。

结果如下:可见,此时根本没有信号显示了。

用Scilab增强Excel的数据处理能力

用Scilab增强Excel的数据处理能力

用Scilab增强Excel的数据处理能力
李保源
【期刊名称】《福建电脑》
【年(卷),期】2008(024)008
【摘要】借助R/Scilab(D)COM Sever,可将Excel输入数据和存储数据非常方便的优点和Scilab强大的数值计算和计算结果可视化能力相结合,从而大大增强Excel的数据处理能力.在Excel中调用Scilab,用VBA编程极其方便,只需改动供Scilab调用的命令字符串,就可实现原Excel无法办到的作图和计算功能.
【总页数】2页(P185,161)
【作者】李保源
【作者单位】南通大学理学院,江苏,南通,226007
【正文语种】中文
【中图分类】TP3
【相关文献】
1.插值法在测井数据处理中的应用——基于Excel的数据处理 [J], 孔玉霞;廖明先;易娟子
2.巧用Excel提升工资管理数据处理能力 [J], 郜向茹
3.加强数据处理训练,增强学生综合实验能力 [J], 剡文杰
4.浅谈定量包装秤的检定及Excel数据处理应用 [J], 任志斌
5.关于药学类专业开设实验数据处理方法课程的必要性与可行性的探讨——基于Microsoft Excel软件的数据处理方法的活用 [J], 白文明;陈建平;牡丹;武旭芳;布仁
因版权原因,仅展示原文概要,查看原文内容请购买。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
8
北京邮电大学信息与通信工程学院
综合实验报告
ylabel ( "y(t)(V)" ) mtlb_axis( [0,3, -3,+3] ); else plot ( f ,abs(Y) ) xlabel ("f(kHz)" ) ylabel ( "|Y(f)|" ) mtlb_axis( [-22,22, 0,40] ) end //end if end //end i
六、分析讨论 根据以上图形进行分析,总结如下几点: (1)AM、DSB-SC、SSB 调制信号的包络就是基带信号。在 AM中不存在 过调幅的现象。通过时域可以看出,如果把包络检测出来,再通过直流检波 器就可以还原出基带信号。 (2) 原始信号 基带频谱 AM 信号 双边带抑制载波调幅信号的频谱加上 离散的大载波分量 DSB-SC 抑制了载波分量 SSB 只包上含边带频率。
实验二 一、实验目的 1、 了解 FM 信号的 scilab 产生的原理和调制的实现方法,特别是积分计算 如何在计算机种通过累加运算完成。 2、 了解FM信号波形及频谱的特点,观察原信号与包络之间关系。 3、 掌握 FM 信号的通信数学模型,通过仿真出来的频率谱与卡松公式计算 出来的理论上的频谱宽度做比较。 4、 进一步深入了解 scilab软件的编程语法与工具箱。 二、实验内容 假 设 基 带 信 号 为 m( t ) s i n ( 2 0 0t 0 ) 2 cos (1t0 0 0 ) 4 s itn ( 5 , 00载波频 ) 率为 3 40kHz,仿真产生 FM 信号,观察波形与频谱,并与卡松公式做对照。FM 的 频率偏移常数是 5kHz/V。 三、仿真模型
北京邮电大学信息与通信工程学院
综合实验报告
2010—2011 学年第二学期 《通信原理软件》实验报告
专业班级 姓 学 名 号
信通 14 班 林承绩 09210405 信通院实验中心 2011 年 11 月 29 日
开课系室 报告日期
பைடு நூலகம் 北京邮电大学信息与通信工程学院
综合实验报告
摘 要
本文对基带信号 mt=sin2000 π t+2cos(1000 π t) 进行了载波 频率为 20kHz , AM 、 DSB-SB 、 SSB 的仿真,并且对基带信号 mt=sin2000πt+2cos 1000πt+4sin(500πt+π/ 3)进行载波频率为 40kHz 的 FM 信号(FM 的频率偏移常数是 5kHz/V),还仿真测量 了占空比为 25%、50%、75%以及 100%的单双极性归零码波形及 其功率谱,分析不同占空比对仿真结果的影响,以及仿真了测量滚降 系统系数为α =0.25 的根升余弦滚降系统的发送功率谱密度及眼图。 最后进行了对基带传输的综合仿真。 关键词:基带传输、调频调幅、单双极性(不)归零码
AM 信 号 的 表 达 式 为 s(t ) Ac [1 m(t )]cos(2 fct ) , 其 傅 里 叶 频 谱 为 A S ( f ) c [M ( f f c ) ( f f c ) M ( f f c ) ( f f c )] ,由波形图知其信号的包络为 2 正值,频谱分布为在 f f c , f fc , f fc f m , f fc f m , f fc f m ,
9
北京邮电大学信息与通信工程学院
综合实验报告
(3)由调制信号的波形和频谱可见,AM信号比较容易从包络中恢复出原基 带信号,在工程上经常用包络检波器进行解调。两外两种从包络中恢复出基 带信号比较复杂,在工程上经常用相干解调方法进行解调。 七.所遇问题及解决: 好久没有用matlab了,也许是刚开始scilab暂时有点迁移不过来。 原理上问题不大,编程挺顺利。但是参考例题的语句太罗嗦了,我进行了优 化。 八、实验总结 通过本次实验,我对基带信号的 AM,DSB-SC,SSB调制的原理和 SCILAB 的实现有了更加清晰的了解。通过自己亲手编写程序,调试,观察波形,很 好的巩固了以前通信原理课上学习的原理。我学过matlab,所以对 SCILAB 的编程基本语法和编程结构比较熟悉。只要原理清晰,就能较好较省时的完 成实验。
f fc f m 处。即双边带抑制载波调幅信号的频谱。
SSB 信号:
SSB 信号的表达式为: sSSB (t ) Ac m(t ) cos(2 f ct ) Ac m(t )sin(2 f ct ) ,其傅里叶频 谱为上边带 S ( f ) Ac M ( f fc )u( f f c ) Ac M ( f f c )u ( f f c ) ,上边带频谱分布在 f f c f m , f fc f m 处。 四、仿真设计 (1)仿真思想 1)由于对信号的截短造成与理想的差别,对离散序列进行傅里叶变化时会有 一个较小的虚部,避免产生误差故需要取其实部; 2) SSB 调制中希尔伯特变换可在频域进行, 然后再把频域进行傅里叶反变换 得到时域信号的希尔伯特变换。 通过工具箱进行画图并分析不同调制方法的正确性与差异性。 五、实验程序及结果 //mt=sin(2000πt) +2cos(1000πt) //载波频率为 20kHz clear all exec t2f.sci ; exec f2t.sci ; N=2^12; //采样点数 fs=80; //khz,采样速率 Bs=fs/2; //系统带宽 T=N/ fs; //截短时间 t=-T/2+[0:N-1]/fs; //时域采样点 f=-Bs+[0:N-1]/T; //频域采样点 f1=1;f2=0.5; A1=1; A2=2; phi1=0;phi2=0; //待观测正弦波的频率、
基带信号 m(t)对载波进行线性调频,调制信号的表达
10
北京邮电大学信息与通信工程学院
综合实验报告
式为: 在离散时间下,FM 信号的表达式为:
其中 si、mk 是 FM 信号和基带信号的采样点值。通过离散点的求和来求连 续信号的积分,这要求离散点足够多,也即采样间隔足够小,使得求和能较 好的逼近积分。 三、编程原理与思想: t 调频信号为 s(t ) Ac cos[2 fc (t )] ,其中 (t ) 2 K f m( )d 因此,只要先积分求 得 (t ) ,进而求得 s(t),即是调频信号。经过傅里叶变换,得到频谱图。 四、代码及截图 //mt=sin(2000πt)+2cos(1000πt)+4sin(500πt+π/3)载波频率为 20kHz //仿真产生 FM 信号,FM 的频率偏移常数是 5kHz/V clear all exec t2f.sci ; exec f2t.sci ; N=2^12; //采样点数 fs=160; //khz,采样速率 dt=1/fs; Bs=fs/2; //系统带宽 T=N/ fs; //截短时间 t=-T/2+[0:N-1]/fs; //时域采样点 f=-Bs+[0:N-1]/T; //频域采样点 f1=1;f2=0.5;f3=0.25; A1=1;A2=2;A3=4; phi1=0;phi2=0;phi3=%pi/3;// 待观测正弦波的频率、幅度和初相 Kf=5; //kHz/V ; fc=40; //kHz //以上是初始化参数设置 m=A1*sin((2*%pi)*f1*t+phi1)+A2*cos((2*%pi)*f2*t+phi2)+A3*sin((2* %pi)*f3*t+phi3); M=t2f(m,fs); //傅里叶变换 mm=real(f2t(M,fs)) ; //傅里叶反变换 phi=(2*%pi) *Kf*cumsum(m)*dt ; //cumsum 为积分 dt? fm=cos((2*%pi)*fc*t+phi) ; //FM 信号 FM=t2f (fm,fs) ; //以上是仿真计算部分,以下是绘图部分 xset("window",1) subplot(2,2,1)
北京邮电大学信息与通信工程学院
综合实验报告
实验一 一、实验目的 1、初步了解SCILAB软件的使用,能写出简单的小程序,看懂代码。 2、理解连续信号在计算机中的表示形式,理解时域与频域的采样定理,理解 傅里叶变换在信号分析中的重要作用。 3、了解 AM、DSB-SC、SSB 信号的具体的数学模型,并了解用 SCILAB 产 生的原理和调制的实现方法。 4、观察AM、DSB-SC、SSB信号波形及振幅频谱的特点。 5、 掌握 AM、 DSB-SC、 SSB 信号调制的通信模型和原理, 学会使用 SCILAB 软件对信号进行仿真的原理和设计方法。 二、实验内容 假设基带信号为m(t)=sin(2000*pi*t)+2cos(1000*pi*t),载波频率为 20KHZ,仿真 AM、DSB-SC、SSB信号,观察其波形和频谱。 三、仿真模型 四、AM、DSB-SC、SSB 信号调制原理: 这三种调制,都是用频率为fm的基带信号对频率fc为的载波信号进行调制, 得到搬移到频率为fc的调制信号,在具体的实现上有所不同。 m(t) AM A Vcmcoswct
f fc f m 处。即双边带抑制载波调幅信号的频谱加上离散的大载波分量。
DSB-SC 信号: DSB-SC 信 号 的 表 达 式 为 : s(t ) m(t ) Acco s(2 fct ) , 其 傅 里 叶 频 谱 为 :
6
北京邮电大学信息与通信工程学院
综合实验报告
Ac [M ( f f c ) M ( f f c )] ,由波形图知,时域中此信号有π相位的突变, 2 信号包络与 m(t) 不同。频谱分布为在 f fc f m , f fc f m , f fc f m , S( f )
7

北京邮电大学信息与通信工程学院
相关文档
最新文档