人大经济论坛 › 标签 › 方程组

标签: 方程组经管大学堂：名校名师名课

相关帖子	版块	作者	回复/查看	最后发表

联立方程组的3sls估计	Stata专版	hubifeng？ 2013-6-28	23 32595	Bullishheaven 2017-12-23 19:17:27
面板数据联立方程组需要做什么检验么，我用的3sls估计，多谢了	计量经济学与统计软件	446804599 2013-3-30	3 6966	浅蓝倒吊 2016-7-9 16:20:50
Stata13.0，新鲜试用中的一些小发现~	Stata专版	eblog 2013-8-18	51 23038	houyunhuang 2015-10-31 14:38:43
面板数据的联立方程组模型	EViews专版	Centigrade23 2013-1-23	9 10372	maodongjun 2015-9-17 14:44:39
哪位同学有SAS求解非线性方程组的案例啦，使用哪种迭代方法都可以	SAS专版	remeva 2013-7-16	1 1681	bakoll 2015-3-3 20:26:13
求助	MATLAB等数学软件专版	remeva 2013-7-15	0 1356	remeva 2013-7-15 19:48:32
mathcad未解出的方程组？大神帮我看看是为啥（悬赏50金币） - [悬赏 50 个论坛币]	计量经济学与统计软件	萱草椰子 2013-5-22	1 1464	萱草椰子 2013-7-11 18:16:28
基于联立方程组模型的辽宁省能源供需研究 - [!reward_solved!]	求助成功区	帅到被追杀 2013-7-11	3 941	veronica0620 2013-7-11 11:42:13
联立方程与IV	统计软件培训班VIP答疑区	peyzf 2013-6-29	1 1524	arlionn 2013-6-30 23:01:04
联立方程组模型中的迭代似不相关回归ISUR如何在STATA中实现？	Stata专版	skydeeye 2013-6-13	5 5322	h3327156 2013-6-26 23:13:58
金融数量分析：基于MATLAB编程一书的pdf&所有程序集	MATLAB等数学软件专版	shufelyc 2013-3-17	7 5545	七夜哥哥 2013-5-20 13:11:38
关于用CGARCH模型研究股票市场波动率的几个问题	EViews专版	wusi343922710 2013-5-13	1 3216	ermutuxia 2013-5-15 13:52:44
求助（急）：用1stopt解决非线性方程组的解（已有程序，需要较高版本的1stopt帮忙算出 - [悬赏 10 个论坛币]	MATLAB等数学软件专版	xiaokang312135 2013-5-4	0 3357	xiaokang312135 2013-5-4 19:55:07
MATLAB解非线性方程组，请帮忙看下代码哪里错了 - [!reward_solved!]	MATLAB等数学软件专版	萱草椰子 2013-4-27	6 1936	美髯客 2013-4-28 11:06:18
急急急！！在联立方程组中，内生变量与外生变量的交互项是内生变量吗	爱问频道	sophiaxbg 2013-3-5	2 2036	WinniePooh2008 2013-3-25 09:18:34
matlab 求解方程组 - [悬赏 100 个论坛币]	MATLAB等数学软件专版	anne_cheng 2013-3-4	4 1354	jacobian_T 2013-3-11 13:16:16
sureg（似不相关回归）：insufficient observations？ - [!reward_solved!]	Stata专版	TVB 2013-1-26	3 6366	蓝色 2013-1-26 13:37:43
MATLAB解方程组问题	爱问频道	游愚 2012-8-20	8 3186	游愚 2012-8-20 17:27:41

更多...

相关日志

分享核天体物理: accumulation 2015-7-4 14:42; 在本论文中，我们在相对论平均场近似框架下利用推广的手征强子模型（即FST模型）和非线性Walecka模型系统讨论了中子星内的奇异相（包括K 介子凝聚，超子和奇异夸克物质）及其对中子星整体性质的影响。在相对论平均场理论的有效模型下，重子（包括核子和超子）和K 介子之间通过交换δ,ω,和ρ介子来传递相互作用。假设中子星内部为无限大的均匀的各向同性的核物质，利用平均场的方法取介子场为它们在均匀核物质基态中的平均值，矢量介子场只存在时间分量，而ρ介子场只存在同位旋第三分量。从有效模型的拉氏密度出发，利用标准的变分方法得到重子场、轻子场以及介子场的运动方程，结合中子星的系统条件（粒子化学势的ß平衡条件、系统整体电中性条件以及重子数守恒条件）组成了完备的非线性方程组进行自洽求解，从而得到中子星的成分含量、压强以及能量密度等性质。把得到的中子星状态方程代入Tolman-Oppenheimer-Volkoff (TOV) 方程组进行计算，我们可以得到中子星的质量-半径关系。我们系统讨论了K 介子凝聚和超子在中子星内的生成以及这些奇异相对中子星整体性质的影响。计算结果表明，K 介子凝聚和超子明显地软化中子星物质的状态方程，从而降低中子星的最大质量。由于中微子在前中子星和中子星的产生和演化过程起至关重要的作用，因此，我们详细分析了前中子星物质内的中微子束缚对前中子星内K 介子凝聚和前中子星的质量－半径关系的影响。通过讨论可以得到结论，中微子束缚把1K 介子凝聚的临界密度增加到更高的密度；同时中微子束缚使得前中子星的状态方程变硬，从而增大前中子星的最大质量。我们还计算了中子星内奇异夸克物质的生成及其对中子星性质的影响。首先，我们利用简单的MIT 袋模型描述中子星内的夸克相，结果表明，奇异夸克相的存在排除了中子星内出现1K 介子凝聚的可能性，同时奇异夸克相降低中子星的最大质量。当考虑夸克物质中的色超导效应，我们利用色味锁定(CFL)夸克模型，结果表明：当选取较小的袋常数，中子星内部能够生成CFL夸克物质，在中子星核心存在着强子-CFL夸克混合相；与未成对的夸克物质相比，色超导效应能够软化系统的状态方程，从而降低中子星的最大质量。同时，我们还在CFL夸克模型中引入了微扰QCD修正参数®s 来考虑夸克间通过传递胶子进行的相互作用，计算表明，前中子星的最大质量和相应的半径对QCD修正参数并– I –不敏感。我们还计算了超子之间（Y-Y）的相互作用强度对（前）中子星奇异相以及整体性质的影响。随着天文观测技术的发展，对中子星质量和引力红移等物理量的观测能够提供重要的限制条件来验证中子星理论模型的准确性，因此，我们仔细分析了中子星的天文观测数据，选取中子星的观测质量和引力红移做为中子星理论模型的限制条件。通过比较理论计算结果与中子星观测的限制条件，从而研究判断有效模型的适用性以及奇异相能否在中子星内部生成，结果表明，在严格的观测质量和红移限制条件下,K 介子凝聚能够在中子星内部出现；同时质量和红移限制条件不利于超子和奇异夸克物质在中子星内部的生成。; 个人分类: 物理学|0 个评论

分享线性方程组—三角方程组: accumulation 2015-4-9 11:20; 1.最简单的是对角形方程与三角形方程； 2.把一般性方程化为三角形方程； 3.输入结果：fin.txt；输出结果：fout.txt； 4.上三角方程组与下三角方程组； module tri_eq !----------------------------------------module coment ! Version : V1.0 ! Coded by : syz ! Date : 2010-4-8 !----------------------------------------------------- ! Description : 用于解上、下三角形线性方程组的回带方法模块 ! !----------------------------------------------------- ! Contains : ! 1. ! 2. !----------------------------------------------------- contains subroutine uptri(A,b,x,N) !---------------------------------subroutine comment ! Version : V1.0 ! Coded by : syz ! Date : 2010-4-8 !----------------------------------------------------- ! Purpose : 上三角方程组的回带方法 ! Ax=b !----------------------------------------------------- ! Input parameters : ! 1. A(N,N)系数矩阵 ! 2. b(N)右向量 ! 3. N方程维数 ! Output parameters : ! 1. x 方程的根 ! 2. ! Common parameters : ! !---------------------------------------------------- implicit real*8 (a-z) integer ::i,j,k,N real*8::A(N,N),b(N),x(N) x(N)=b(N)/A(N,N) !回带部分 do i=n-1,1,-1 x(i)=b(i) do j=i+1,N x(i)=x(i)-a(i,j)*x(j) end do x(i)=x(i)/A(i,i) end do end subroutine uptri subroutine downtri(A,b,x,N) !---------------------------------subroutine comment ! Version : V1.0 ! Coded by : syz ! Date : 2010-4-9 !----------------------------------------------------- ! Purpose : 下三角方程组的回带方法 ! Ax=b !----------------------------------------------------- ! Input parameters : ! 1. A(N,N)系数矩阵 ! 2. b(N)右向量 ! 3. N方程维数 ! Output parameters : ! 1. x 方程的根 ! 2. ! Common parameters : ! !---------------------------------------------------- implicit real*8 (a-z) integer ::i,j,N real*8 ::A(N,N),b(N),x(N) x(1)=b(1)/a(1,1) do k=2,N x(k)=b(k) do i=1,k-1 x(k)=x(k)-a(k,i)*x(i) end do x(k)=x(k)/a(k,k) end do end subroutine downtri end module tri_eq !############################################################## module driver !----------------------------------------module coment ! Version : V1.0 ! Coded by : syz ! Date : !----------------------------------------------------- ! Description : 驱动函数模块 ! !----------------------------------------------------- ! Parameters : ! 1. ! 2. ! Contains : ! 1. dir_main 驱动函数入口 ! 2. dri_up 当读到关键字uptri 时启动该函数 ! 3. dri_down 当读到关键字 downtri时启动该函数 !----------------------------------------------------- ! Post Script : ! 1. ! 2. !----------------------------------------------------- contains subroutine dri_main() !---------------------------------subroutine comment ! Version : V1.0 ! Coded by : syz ! Date : 2010-4-9 !----------------------------------------------------- ! Purpose : 驱动程序入口函数 ! !----------------------------------------------------- ! Input files : ! 1. fin.txt 准备数据 ! 2. ! Output files : ! 1. fout.txt 输出结果文件 ! 2. ! !---------------------------------------------------- implicit real*8 (a-z) integer ::ioerr character (20)::upordown open (unit=11,file='fin.txt',iostat=ioerr) open (unit=12,file='fout.txt') do read (11,*)upordown ! 读输入文件 ! 当读到关键字uptri时启动上三角矩阵计算 ! 当读到关键字downtri时候启动下三角矩阵计算 if ( upordown(1:5) == 'uptri' ) then call dri_up() else if (upordown(1:)=='downtri') then call dri_down() end if if (ioerr/=0) exit !读到文件结束时，退出读文件 end do end subroutine dri_main subroutine dri_up() !---------------------------------subroutine comment ! Version : V1.0 ! Coded by : syz ! Date : 2010-4-9 !----------------------------------------------------- ! Purpose : 启动上三角阵的计算 ! !----------------------------------------------------- use tri_eq implicit real*8 (a-z) integer,parameter ::N=4 integer ::i,j real*8 ::A(N,N),b(N),x(N) read (11,*)((A(i,j),j=1,N),i=1,N) !读入B向量 read (11,*) b call uptri(A,b,x,N) write (12,101)x 101 format (T5,'上三角形方程组的解',/,T4,'x=',4(/F12.8)) end subroutine dri_up subroutine dri_down() !---------------------------------subroutine comment ! Version : V1.0 ! Coded by : syz ! Date : 2010-4-9 !----------------------------------------------------- ! Purpose : 启动下三角阵的计算 ! !----------------------------------------------------- use tri_eq implicit real*8 (a-z) integer,parameter ::N=4 integer ::i,j real*8 ::A(N,N),b(N),x(N) read (11,*)((A(i,j),j=1,N),i=1,N) !读入B向量 read (11,*) b call downtri(A,b,x,N) write (12,101)x 101 format (/,T5,'下三角形方程组的解',/,T4,'x=',4(/F12.8)) end subroutine dri_down end module driver !----------------------------------------------------- program main !--------------------------------------program comment ! Version : V1.0 ! Coded by : syz ! Date : 2010-4-9 !----------------------------------------------------- ! Purpose : 计算上、下三角形方程组 ! !----------------------------------------------------- ! In put data files : ! 1. fin.txt 输入方程系数 ! 2. ! Output data files : ! 1. fout.txt 计算结果 ! 2. !----------------------------------------------------- ! Post Script : ! 1. 需要准备输入数据 ! ! 2. 由主函数启动驱动程序进行计算 !----------------------------------------------------- !use tri_eq use driver !调用驱动函数 call dri_main() end program main; 个人分类: 裂变模型|0 个评论

分享计算物理学——线性代数方程组的直接求解: accumulation 2015-3-9 17:02; 1.追赶法解三对角线性方程组（1）对系数矩阵A进行克劳特分解A=LU0 （2）追赶法顺利实施的条件（3）追赶法算法的伪代码 2.列主元消元法（1）列主元消去法（2）全主元消去法（3）列主元消元法算法; 个人分类: 原子核物理|0 个评论

分享计算物理学——线性代数方程组的直接求解: accumulation 2015-3-9 16:35; 1.初等消去阵：（1）初等消去阵或高斯变换；（2）初等消去阵的性质；（3）高斯消元与回代过程；（4）高斯消元法的步骤（伪代码）；（5）高斯消元：最终三角方程组；（6）增广矩阵与系数矩阵；（7）高斯消元法的计算复杂度； 2.矩阵的LU分解（1）初等变换矩阵；（2）矩阵A分解为一个下三角和上三角矩阵的乘积；（3）L及其逆矩阵；（4）矩阵的直接LU分解；（5）杜利特尔分解与克劳特分解；（6）非奇异矩阵不一定存在LU分解；（7）矩阵的直接LU分解的用途；（8）对角占优矩阵；（9）对角占优矩阵在高斯消元算法中稳定; 个人分类: 原子核物理|0 个评论

分享计算物理学——线性代数方程组的直接求解: accumulation 2015-3-9 15:46; 1.引言：系数矩阵、解向量与右端向量； 2.线性方程组求解的重要性：（1）线性方程组求解与矩阵的特征值问题；（2）物理学与经济学实例；（3）微分方程的差分化；（4）数值求解的问题； 3.线性代数的基本概念回顾：（1）线性相关与线性无关；（2）非奇异矩阵的4种定义形式：逆矩阵、行列式、秩序、非零向量；（3）方程组Ax=b解的存在性与唯一性；（4）奇异矩阵与非奇异矩阵的例子； 4.数值求解线性方程组：（1）克莱姆法则；（2）计算复杂度； 5.几种特殊的矩阵：（1）对角矩阵；（2）三对角矩阵；（3）上下三角形矩阵；（4）对称矩阵；（5）对称正定矩阵；（6）对称半正定矩阵；（7）正交矩阵；（8）几种特殊矩阵的性质； 6.三角形方程组的求解：（1）上三角方程组的回代算法；（2）下三角方程组的前代算法； 7.线性方程组的三种基本等价变换：（1）非奇异线性变换的简化；（2）简化为：对角矩阵和三角形矩阵；（3）PAx=Pb——PA=U；（4）上三角化与对角化；（5）特征值问题：相似变换；（6）A相似于B：P(-1)*A*P； 8.线性方程组的求解的两大类方法：（1）直接解法—高斯消元；（2）迭代解法—构造向量序列；（3）迭代解法的分类；; 个人分类: 原子核物理|0 个评论