应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 其他 » Matlab习题求解
41/1返回列表
查看: 820  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

杨幸遇

金虫 (小有名气)

[求助] Matlab习题求解 已有1人参与

小弟不才，Matlab的程序代码出现问题，求大神给出代码模板或者此题程序代码。期末试题啊，求给力！感激涕零！
已知水在不同温度下的密度如下表所示,请用一元三点拉格朗日插值法计算水在12,33和66℃时的密度。
T（℃）           0          20          40          60          80         100
ρ（kg/m3）999.9        998.2        992.2        983.2        971.8        958.4
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
翅膀的痕迹
1楼 2014-05-26 22:11:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
杨幸遇: 金币+15, ★★★★★最佳答案 2014-05-28 08:55:07
CODE:
function Lag_interp
clear all;clc
Ti=[12 33 66];
T=[ 0          20          40          60          80         100];
p=[999.9        998.2        992.2        983.2        971.8        958.4];

for i=1:length(Ti);
    p_answer(i)=polyinterp(T,p,Ti(i));
end
Result=[Ti' p_answer']
   


function v=polyinterp(x,y,u)
n=length(x);
v=zeros(size(u));
for k=1:n;
    w=ones(size(u));
    for j=[1:k-1 k+1:n];
        w=(u-x(j))./(x(k)-x(j)).*w;
    end
v=v+w*y(k);
end

计算结果：

Result =

   12.0000  999.5037
   33.0000  994.6881
   66.0000  980.0144
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
2楼2014-05-27 18:33:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

杨幸遇

金虫 (小有名气)
引用回帖:
2楼: Originally posted by 月只蓝 at 2014-05-27 18:33:19
function Lag_interp
clear all;clc
Ti=;
T=;
p=;

for i=1:length(Ti);
    p_answer(i)=polyinterp(T,p,Ti(i));
end
Result=
   


function v=polyinterp(x,y,u)
n=length(x);
v=zeros(size(u)) ...

万分感谢！金币送上！
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

杨幸遇
翅膀的痕迹
3楼2014-05-28 08:54:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

杨幸遇

金虫 (小有名气)
送红花一朵
引用回帖:
3楼: Originally posted by 杨幸遇 at 2014-05-28 08:54:54
万分感谢！金币送上！...

大神 我可否再问你一道题啊，实在是对MATLAB编程无感，只为完成作业。
赞一下(1人) 回复此楼
翅膀的痕迹
4楼2014-05-28 16:07:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 杨幸遇 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭