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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 材料区 » 材料综合 » 塑料与橡胶 » 二维拟合程序改为三维的，程序看不懂，不知道怎么改，求各路大神指教批评
11/1返回列表
查看: 187  |  回复: 0
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

fengniao2011

新虫 (初入文坛)

[求助] 二维拟合程序改为三维的，程序看不懂，不知道怎么改，求各路大神指教批评

clc;clear;
%曲线拟合
x=[0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0];
y=[0, 4, 5, 14, 15, 14.5, 14, 12, 10, 5, 4];
scatter(x,y,'filled');
len_x = length(x);
max_x = max(x);
min_x = min(x);
num = 100;
delta = (max_x-min_x)/num;
x_f =[];
f =[];
max_delta = (max_x-min_x)*3/10;
for i=0:num
    x_val = min_x + i*delta;
    x_f = [x_f,x_val];
    A = zeros(2,2);
    B = [];
    for j=1:len_x
        s = abs(x(j)-x_val)/max_delta;
        if s<=0.5
            w = 2/3-4*s^2+4*s^3;
        elseif s<=1
            w = 4/3-4*s+4*s^2-4*s^3/3;
        else
            w = 0;
        end
        A = A + w*[1;x(j)]*[1,x(j)];
        B = [B,w*[1;x(j)]];
    end
    f =[f,[1,x_val]*inv(A)*B*y'];
end
hold on
plot(x_f,f,'r');
回复此楼

» 猜你喜欢
1楼 2014-08-20 16:04:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子，希望有用哦~

科研从小木虫开始，人人为我，我为人人
相关版块跳转 我要订阅楼主 fengniao2011 的主题更新
11/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭