结构力学求解器!!!

&& 查看话题
【求助】 求MAtlab求解程序!!!
小弟是化学专业的,在验证化学公式时由于数学基础不好,不知道怎么弄。能否用matlab编程算出下面的三个方程组,此外要算出置信度或者线性之类的东西,因为我要验证哪个公式符合化学规律,谢谢虫子们帮忙啊!!!
3.98*a+59.74*b+1/9*c+9d=67.81
21.76*a+46.39*b+3/7*c+7/3*d=81.19
51.25*a+28.84*b+c+d=93.49
75.31*a+14.43*b+7/3*c+3/7*d=86.71
93.25*a+6.68*b+9*c+1/9*d=97.45
3.98*a+59.74*b+1/9*c=67.81
21.76*a+46.39*b+3/7*c=81.19
51.25*a+28.84*b+c=93.49
75.31*a+14.43*b+7/3*c=86.71
93.25*a+6.68*b+9*c=97.45
1/9*a+9b=4.09
3/7*a+7/3*b=13.04
7/3*a+3/7*b= -3.03
9*a+1/9*b=-2.48
请一定要算出置信度,因为我要凭此看这三个方程组哪个有有效解。谢谢帮助,不胜感激!!!
=solve('等式1','等式2'...,'等式8'),你用这个能求取有效解,应该你试试 试试这个求解线性方程组
solve,linsolve
%矩阵的行之间用分号隔开,元素之间用逗号或空格
X=zeros(4,1);%建立一个4元列向量
X=linsolve(A,B)
diff(fun,var,n):对表达式fun中的变量var求n阶导数。
例如:F=sym('u(x,y)*v(x,y)'); %sym()用来定义一个符号表达式
diff(F); %matlab区分大小写
pretty(ans) %pretty():用习惯书写方式显示变量;ans是答案表达式 第一个方程,5方程4未知数,条件过剩。
第一个方程,5方程3未知数,条件过剩。
第一个方程,5方程2未知数,条件过剩。 问楼主一个问题,你这里的置信度是指什么:sweat: Originally posted by bluesine at
第一个方程,5方程4未知数,条件过剩。
第一个方程,5方程3未知数,条件过剩。
第一个方程,5方程2未知数,条件过剩。 可以求最小二乘解。 Originally posted by saladin983 at
可以求最小二乘解。 能告诉我怎么解吗,我这个是方程过剩的,所以求的是大体的结果,就是那个最小二乘法,能详细些吗?不胜感激! Originally posted by bluesine at
问楼主一个问题,你这里的置信度是指什么:sweat: 就是matlab上的那个“confidence bounds” Originally posted by swenson at
能告诉我怎么解吗,我这个是方程过剩的,所以求的是大体的结果,就是那个最小二乘法,能详细些吗?不胜感激! 是基本运算,给你找了一点中文说明。其实Matlab的说明文档写的很好,如果耐心看的话能学到很多的。另外可以参考一下lscov这个内建函数。
Matlab求解线性方程组
AX=B或XA=B
在MATLAB中,求解线性方程组时,主要采用前面章节介绍的除法运算符“/”和“\”。如:
X=A\B表示求矩阵方程AX=B的解;
X=B/A表示矩阵方程XA=B的解。
对方程组X=A\B,要求A和B用相同的行数,X和B有相同的列数,它的行数等于矩阵A的列数,方程X=B/A同理。
如果矩阵A不是方阵,其维数是m×n,则有:
m=n 恰定方程,求解精确解;
m>n 超定方程,寻求最小二乘解;
m<n 不定方程,寻求基本解,其中至多有m个非零元素。
针对不同的情况,MATLAB将采用不同的算法来求解。 Originally posted by saladin983 at
是基本运算,给你找了一点中文说明。其实Matlab的说明文档写的很好,如果耐心看的话能学到很多的。另外可以参考一下lscov这个内建函数。
Matlab求解线性方程组
AX=B或XA=B
在MATLAB中,求解线性方程组时 ... 我编了,可不知道哪错了,那东西看不太懂,你能给我编好吗?对你应该很简单,帮个忙吧!谢谢!! Originally posted by swenson at
我编了,可不知道哪错了,那东西看不太懂,你能给我编好吗?对你应该很简单,帮个忙吧!谢谢!! 把程序发上来看一眼吧,包含系数矩阵和右端向量的。 Originally posted by saladin983 at
把程序发上来看一眼吧,包含系数矩阵和右端向量的。
%f = fittype('A*=b');
% = fit(A,b,f)
A1=A'*A;b=A'*b;
x=inv(A1)*b
这种能算结果,但是我要看结果的置信度,那个怎么算啊,怎么知道结果的可靠性啊 Originally posted by swenson at
%f = fittype('A*=b');
% = fit(A,b,f)
A1=A'* ... 这样用也能算,没有错,不过通常我会用A\b来求最小二乘解,反斜杠这个运算更稳定更快捷。
这么小的样本做回归分析我觉得没法讨论置信度,当然,我对概率统计的东西本来就不熟。Matlab有个函数regress,统计工具箱里的,说明如下:
REGRESS Multiple linear regression using least squares.
& & b = REGRESS(y,X) returns the vector of regression coefficients, b,
& & in the linear model&&y = Xb, (X is an nxp matrix, y is the nx1
& & vector of observations).
& &&&= REGRESS(y,X,alpha) uses the input, ALPHA
& & to calculate 100(1 - ALPHA) confidence intervals for B and the
& & residual vector, R, in BINT and RINT respectively.&&The vector
& & STATS contains the R-square statistic along with the F and p
& & values for the regression.
& & The X matrix should include a column of ones so that the model
& & contains a constant term.&&The F and p values are computed under
& & the assumption that the model contains a constant term, and they
& & are not correct for models without a constant.&&The R-square
& & value is the ratio of the regression sum of squares to the
& & total sum of squares.
你的模型里没有常数项,所以只能算.我改了变量名,以对应函数说明。过程如下:
= regress(y,X,0.05);
得到的b是最小二乘解,而bint是置信度为95%(可以自行调整)的置信区间,看了输出的数据你就会知道这个置信区间太大,没什么意义,我想这是因为数据太少的缘故。如果非要比的话,可以计算相同置信度下的置信区间的大小,比如
norm(bint(:,1)-bint(:2))(计算矩阵bint第一列和第二列的差,然后求范数)。不过我觉得这么做意义不大。如何衡量你这个问题里线性拟合的有效性,这个问题我回答不了,需要做统计的人来说,或者,你提供一个类似的样例。 Originally posted by saladin983 at
这样用也能算,没有错,不过通常我会用A\b来求最小二乘解,反斜杠这个运算更稳定更快捷。
这么小的样本做回归分析我觉得没法讨论置信度,当然,我对概率统计的东西本来就不熟。Matlab有个函数regress,统计 ... 线性拟合的有效性我有一个样例,你看看,对你的耐心回答,我很感激!!
plot(x,y,'*')
f = fittype('1 +a + b*x- ((1+a+b*x)^2-4*b*x)^1/2');
= fit(x,y,f)
% 显示拟合曲线与数据点的结果
plot(c2,'m')&&% 绘制拟合曲线
hold on,plot(x,y,'*') %绘制原数据点 Originally posted by swenson at
线性拟合的有效性我有一个样例,你看看,对你的耐心回答,我很感激!!
y = [0..5768; 0.....07 ... 我机器上没有安装Curve fitting toolbox,晚些时候再看看。偶尔涉猎一下自己不熟悉的领域,复习一下基础知识,对我自己也有帮助。 Originally posted by saladin983 at
我机器上没有安装Curve fitting toolbox,晚些时候再看看。偶尔涉猎一下自己不熟悉的领域,复习一下基础知识,对我自己也有帮助。 恩,麻烦你了啊! 看了你给的这个例子,问题类型不太一样,这是个一元的非线性回归,而你的问题是多元线性回归。我看了一下Matlab曲线拟合工具箱的文档,没有看到任何关于多元回归的说明,所有的library model都是一元的。所以我认为这个工具箱在你的问题上没有用处。
网上稍微查了一下,regress函数是做多元线性回归最常用的。关于有效度的问题,我翻了这个文档(http://em./uploadfile/0244.ppt),现学了一下多元线性回归的内容,里面倒是有几种检验方法。但是你的模型里面有一个很大的差别,也就是前面说到的,没有常数项(注意:在regress说明档里面就提到了),这个跟一般的线性回归模型都不同,所以他们的一些校验参数似乎也用不上。
看起来这个有效度的参数只能自己定义了,我的设想是norm(Xb-y,2)^2/length(y)(变量定义同13楼),看起来还算是合理,但是我不知道这样做是否有什么理论基础。 Originally posted by saladin983 at
看了你给的这个例子,问题类型不太一样,这是个一元的非线性回归,而你的问题是多元线性回归。我看了一下Matlab曲线拟合工具箱的文档,没有看到任何关于多元回归的说明,所有的library model都是一元的。所以我认 ... 我研究研究,虽然好多定义对我很抽象,谢谢啦!! 刚看到Matlab帮助文档中关于曲线拟合的有效度有一段说明,参见Residual Analysis一节。照这个看来,似乎前面定义的量有一定的道理。 Originally posted by saladin983 at
刚看到Matlab帮助文档中关于曲线拟合的有效度有一段说明,参见Residual Analysis一节。照这个看来,似乎前面定义的量有一定的道理。 具体点好吗?求解!!!_百度知道
求解!!!
atl100.ll丢失,自己下了一个放进去还是不行,安装的墨菲显示它没被指定什么的!求解啊!
又出现这个了:模块atl100.dll可能与您正在运行的windows版本不兼容
我有更好的答案
那里也有这个问题的解决方案
w 不兼容就没办法了
您可能关注的推广
等待您来回答
下载知道APP
随时随地咨询
出门在外也不愁谁来告诉我沙发的合理高度是多少?求解!!!!!!!!!!!!!!_百度知道
谁来告诉我沙发的合理高度是多少?求解!!!!!!!!!!!!!!
提问者采纳
从地面算起,到沙发座面38-----48之间,通常采用42,单位cm。实际上就是小腿(坐下后大腿以下)的长度。所以这个高度是国人小腿的均高。祝您快快乐乐开心每一天!
提问者评价
谢谢您的回答哦!也祝您开心每天!
其他类似问题
为您推荐:
其他6条回答
沙发的总高度一般为850-950mm左右(包括靠背高度)。单独坐高通常为450mm左右。不同风格的沙发高度不一样,古典风格的沙发一般靠背较高,总高约900mm,有些欧式设计的酒吧沙发则有1000mm左右,而现代风格的沙发一般850mm左右高,当然还有更低的了。
你是为个人设计的吗?那样可以量身定做,这是最好的呢!也有固定的尺寸,但是也根据地域人群有分类!设法深度在80~100,坐深55~70 正常状态,希望能帮到您!我是学室内设计的,有问题可以交流一下!
沙发座的一般高度为400mm--450mm,常用高度为420mm,背的高度看款式来定。
一般为380——450,如果是欧式家具的话,会更高,有的能到500多
市场上沙发的尺寸,高度一般是85-90CM左右。
现代的沙发 400mm左右
欧式 差不多450左右
您可能关注的推广
沙发的相关知识
等待您来回答
下载知道APP
随时随地咨询
出门在外也不愁下了魔兽世界地图怎么安装到游戏里去? 求解》 急!!!!!!!_百度知道
下了魔兽世界地图怎么安装到游戏里去? 求解》 急!!!!!!!
提问者采纳
打开游戏安装目录,寻找map 文件夹,把下载好的图复制到这个文件夹里,就行了,然后进入游戏创建游戏在地图里就可以找到你的新地图了
提问者评价
谢了,兄弟》》
其他类似问题
为您推荐:
魔兽世界地图的相关知识
其他1条回答
你是想说魔兽称霸?把地图放在maps文件就行
等待您来回答
下载知道APP
随时随地咨询
出门在外也不愁L/C,在线求解!!! - 外贸单证 -
福步外贸论坛(FOB Business Forum) |中国第一外贸论坛
& L/C,在线求解!!!
UID 1617263
阅读权限 40
L/C,在线求解!!!
L/C客户需要的文件资料
Documents Required
1. MANUALLY SIGNED COMMERCIAL INVOICES IN 3 (SETS/FOLDS)
CERTIFIED BY CHAMBER OF COMMERCE
贸促会出的发票认证应该只有一份,我们可以自己打印三份发票盖章提供给客户吗?
+INVOICES SHOULD SHOW THE BREAK DOWN IN THE VALUE AS FOLLOWING:
FOB VALUE:
FREIGHT PREPAID :
INSURANCE PREMIUM:
TOTAL CIF COST :
2. FULL SET OF ORIGINAL CLEAN ON BOARD OCEAN BILLS OF LADING
MADE OUT TO THE ORDER OF AL RAJHI BANK MARKED FREIGHT PREPAID
AND NOTIFY APPLICANT INDICATING THIS DOCUMENTARY CREDIT NUMBER
AND NAME, ADDRESS, TELEPHONE NO OF THE CARRYING VESSELS AGENT AT
THE PORT OF DISCHARGE.
之前我听说提单上不能显示“CLEAN ON BOARD”,只能显示“SHIP ON BOARD”,到时候给银行审单的时候没问题的吗?
提单收件人要显示“AL RAJHI BANK”?
通知方那一栏要显示信用证号码?
在卸货港那一栏要显示船公司代理的名称,地址,电话?
3. CERTIFICATE OF ORIGIN ISSUED/ATTESTED BY CHAMBER OF COMMERCE
STATING THE NAME AND ADDRESS OF THE MANUFACTURER/OR PRODUCERS,
AND THAT GOODS EXPORTED ARE WHOLLY OF CHINA ORIGIN.
贸促会出的产地证写明制造商的名称跟地址,产品是完全原产。
4. PACKING LIST IN 3 FOLDS.
箱单要3份?
5. MARINE INSURANCE CERTIFICATE/POLICY IN DUPLICATE FOR 110PCT
OF INVOICE VALUE, ENDORSED TO THE ORDER OF AL RAJHI BANK WITH
CLAIMS PAYABLE IN SAUDI ARABIA IN THE CURRENCY OF THE CREDIT,
COVERING INSTITUTE CARGO CLAUSE A, INSTITUTE WAR CLAUSES
(CARGO), INSTITUTE STRIKE CLAUSES (CARGO),HIGHSEAS PIRACY RISKS,
BENEFICIARY WAREHOUSE TO APPLICANT WAREHOUSE, LAND
TRANSPORTATION, TRANSHIPMENT RISKS,RIOT,CIVIL COMMOTION, T.P.N.D
RISKS (IF APPLICABLE)
海上保险凭证是货代出的吗?110PT是什么意思?
福友们,帮我看看这个信用证上客户提到的资料是否容易做到?有没什么特别注意的地方?
[ 本帖最后由 Elsaxule 于
19:50 编辑 ]
(小韦&#9743;:)
单证代理QQ:
UID 1684046
积分 107215
帖子 12878
阅读权限 120
来自 产地证/商会认证/使馆认证
回复 #1 Elsaxule 的帖子
1. MANUALLY SIGNED COMMERCIAL INVOICES IN 3 (SETS/FOLDS)
CERTIFIED BY CHAMBER OF COMMERCE
三份商业发票手签并盖章,由贸促会认证!
( (毛毛))
UID 331320
积分 16530
阅读权限 80
来自 江苏常州
1. MANUALLY SIGNED COMMERCIAL INVOICES IN 3 (SETS/FOLDS)
CERTIFIED BY CHAMBER OF COMMERCE
贸促会出的发票认证应该只有一份,我们可以自己打印三份发票盖章提供给客户吗?
只要提供一份正本,其他二份可以是副本
(小韦&#9743;:)
单证代理QQ:
UID 1684046
积分 107215
帖子 12878
阅读权限 120
来自 产地证/商会认证/使馆认证
回复 #1 Elsaxule 的帖子
3. CERTIFICATE OF ORIGIN ISSUED/ATTESTED BY CHAMBER OF COMMERCE
STATING THE NAME AND ADDRESS OF THE MANUFACTURER/OR PRODUCERS,
AND THAT GOODS EXPORTED ARE WHOLLY OF CHINA ORIGIN.
产地证由商会办理,商会和贸促会是同一个机构,产地证在贸促会出证后在第五栏加盖一个商会章即可!并在产地证货描显示货物源于中国!THAT GOODS EXPORTED ARE WHOLLY OF CHINA ORIGIN.
( (毛毛))
UID 331320
积分 16530
阅读权限 80
来自 江苏常州
2. FULL SET OF ORIGINAL CLEAN ON BOARD OCEAN BILLS OF LADING
MADE OUT TO THE ORDER OF AL RAJHI BANK MARKED FREIGHT PREPAID
AND NOTIFY APPLICANT INDICATING THIS DOCUMENTARY CREDIT NUMBER
AND NAME, ADDRESS, TELEPHONE NO OF THE CARRYING VESSELS AGENT AT
THE PORT OF DISCHARGE.
之前我听说提单上不能显示“CLEAN ON BOARD”,只能显示“SHIP ON BOARD”,到时候给银行审单的时候没问题的吗?&&答:没问题。CLEAN指没有不良批注的提单。
提单收件人要显示“AL RAJHI BANK”?& && & 答:TO THE ORDER OF AL RAJHI BANK
通知方那一栏要显示信用证号码?& && & 答:写在货描处即可。
在卸货港那一栏要显示船公司代理的名称,地址,电话?& &&&答:显示在提单空白处即可
( (毛毛))
UID 331320
积分 16530
阅读权限 80
来自 江苏常州
5. MARINE INSURANCE CERTIFICATE/POLICY IN DUPLICATE FOR 110PCT
OF INVOICE VALUE, ENDORSED TO THE ORDER OF AL RAJHI BANK WITH
CLAIMS PAYABLE IN SAUDI ARABIA IN THE CURRENCY OF THE CREDIT,
COVERING INSTITUTE CARGO CLAUSE A, INSTITUTE WAR CLAUSES
(CARGO), INSTITUTE STRIKE CLAUSES (CARGO),HIGHSEAS PIRACY RISKS,
BENEFICIARY WAREHOUSE TO APPLICANT WAREHOUSE, LAND
TRANSPORTATION, TRANSHIPMENT RISKS,RIOT,CIVIL COMMOTION, T.P.N.D
RISKS (IF APPLICABLE)
海上保险凭证是货代出的吗?110PT是什么意思?
保单有保险公司出。一般货代可以代办。110PT 是指保险加成。发票金额上加10%。
(123叹口气)
UID 1075105
阅读权限 40
这种信用证相对我家的来说 实在是太简单了,很好做的
1.要做CCPIT发票认证,把正本发票及相关的资料递交给贸促会就能做了。一般三五个工作日就能好。另外要的三张给普通正本或是一正两副都可以。
记得要手签并且把发票金额分拆成FOB金额,保险和海运费,最后再写上总的CIF金额即可。
2.提单收货人做:TO THE ORDER OF AL RAJHI BANK&&一般情况下 SHIPPED ON BOARD= CLEAN ON BOARD
DOCUMENTARY CREDIT NUMBER -信用证号 ,这个显示在货描里就可以
NAME, ADDRESS, TELEPHONE NO OF THE CARRYING VESSELS AGENT AT THE PORT OF DISCHARGE. -提单确认时跟货代标注一下要显示目的港代理的详细信息 出单时就会有的。
3.这个是要做一般产地证的。THE NAME AND ADDRESS OF THE MANUFACTURER/OR PRODUCERS,
AND THAT GOODS EXPORTED ARE WHOLLY OF CHINA ORIGIN 这两句显示在货描里或是备注里都可以。 前者是要显示生产商或是制造商的名称和地址,后者是显示原产中国,这么照着抄就行了。
4.箱单三份
5.保单做 TO THE ORDER OF AL RAJHI BANK ,赔付地是 SAUDI ARABIA IN USD(做美元的话),
COVERING INSTITUTE CARGO CLAUSE A, INSTITUTE WAR CLAUSES (CARGO), INSTITUTE STRIKE CLAUSES (CARGO),HIGHSEAS PIRACY& & RISKS, BENEFICIARY WAREHOUSE TO APPLICANT WAREHOUSE, LAND TRANSPORTATION, TRANSHIPMENT RISKS,RIOT,CIVIL
COMMOTION, T.P.N.D RISKS-这段要在赔付条件的备注栏里写上,要不就是不符点。
保险金额是合同金额的1.1倍。即 保险金额=合同金额 x 110%
保单是自己公司出具的,有保险的系统即可。丘博保险,太平保险,太平洋保险 这些都可以的。
这种信用证实在是太简单了,没什么难的。
[ 本帖最后由 januarist 于
10:03 编辑 ]
(123叹口气)
UID 1075105
阅读权限 40
另外还要注意一下就是 到沙特的 选船公司的时候,不要用敌对国家,涉黑国家的。
订舱的时候都标注一下,省的到时麻烦。
(佛山货代)
佛山货代 乐从外贸仓库
UID 1460979
阅读权限 40
原帖由 januarist 于
09:59 发表
这种信用证相对我家的来说 实在是太简单了,很好做的
1.要做CCPIT发票认证,把正本发票及相关的资料递交给贸促会就能做了。一般三五个工作日就能好。另外要的三张给普通正本或是一正两副都可以。
记得要手签并 ... 分析的很详细呢!办理签证,不同大使馆出证时间不同的,注意下。
代理产地证,商会认证,使馆认证
UID 1999570
积分 86906
福步币 17 块
阅读权限 120
来自 原产地证CO,FORM A|FORM F
这些做起来都是比较容易的,放心好了!
UID 1617263
阅读权限 40
回复 #2 ckck 的帖子
贸促会会出三份?
UID 1617263
阅读权限 40
回复 #8 januarist 的帖子
选船公司还有要求的?
(专业货运保险 周先生)
提供专业服务,请加QQ:
UID 1570079
积分 35148
福步币 250 块
阅读权限 120
来自 24小时手机:
回复 #1 Elsaxule 的帖子
保险单被保险人:TO THE ORDER OF AL RAJHI BANK
按发票金额110%投保,
赔付地: SAUDI ARABIA IN THE CURRENCY OF THE CREDIT,
险别要求:COVERING INSTITUTE CARGO CLAUSE A, INSTITUTE WAR CLAUSES
(CARGO), INSTITUTE STRIKE CLAUSES (CARGO),HIGHSEAS PIRACY RISKS,
BENEFICIARY WAREHOUSE TO APPLICANT WAREHOUSE, LAND
TRANSPORTATION, TRANSHIPMENT RISKS,RIOT,CIVIL COMMOTION, T.P.N.D
RISKS (IF APPLICABLE)
(专业货运保险 周先生)
提供专业服务,请加QQ:
UID 1570079
积分 35148
福步币 250 块
阅读权限 120
来自 24小时手机:
保单是由保险公司出的,而不是货代公司出的,
你们可以直接找保险公司办理,或保险代理公司办理,比较专业些.
也可以让货代公司代办,货代公司对保险的专业性不高,收费也相对要高些.
UID 124804
积分 49709
福步币 48 块
阅读权限 120
1. MANUALLY SIGNED COMMERCIAL INVOICES IN 3 (SETS/FOLDS)
CERTIFIED BY CHAMBER OF COMMERCE
+INVOICES SHOULD SHOW THE BREAK DOWN IN THE VALUE AS FOLLOWING:
FOB VALUE:
FREIGHT PREPAID :
INSURANCE PREMIUM:
TOTAL CIF COST :
签署的发票3份,经商会认证,发票上必须分别列明FOB价值,预付运费的金额以及保险金额及CIF总价值。
提交一份认证的发票,其余2份用正本去交单就行。
当前时区 GMT+8, 现在时间是
Powered by D1scuz! && 2001-

我要回帖

更多关于 求解数学题 的文章

 

随机推荐